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Our current knowledge of the quark-gluon plasma in thermodynamical equilibrium is reviewed. 



Abstract 

> 

^r) \ The phase diagram of strongly interacting matter is discussed, with emphasis on the quark-hadron 

phase transition and the color-superconducting phases of quark matter. Lattice QCD results on 
the order of the phase transition, the thermodynamical functions, the heavy quark free energy, 
[ mesonic spectral functions, and recent results for nonzero quark chemical potential are presented. 

Analytic attempts to compute the thermodynamical properties of strongly interacting matter, such 
as perturbation theory, quasiparticle models, "hard-thermal-loop" (HTL)-resummed perturbation 
theory, the Polyakov-loop model, as well as linear sigma models are discussed. Finally, color- 
superconducting quark matter is considered in the limit of weak coupling. The gap equation and 
the excitation spectrum are derived. The solution of the gap equation, gap parameters in various 
color-superconducting phases, and critical temperatures for the transition to normal-conducting 
quark matter are presented. A summary of gluon and photon properties in color superconductors 
is given. 



1 Introduction and Summary 

Quantum chromodynamics (QCD) is the fundamental theory of the strong interaction. QCD is an 
asymptotically free theory [T], i.e., interactions between quarks and gluons become weaker as the mutual 
distance decreases or as the exchanged momentum increases. Consequently, at very large temperatures 
and/or densities, the interactions which confine quarks and gluons inside hadrons should become suf- 
ficiently weak to release them |2j. The phase where quarks and gluons are deconfined is termed the 
quark-gluon plasma (QGP). Lattice QCD calculations have established the existence of such a phase 
of strongly interacting matter at temperatures larger than ~ 150 MeV and zero net-baryon density. 
Depending on the number of quark flavors and the masses of the quarks, the transition between ordi- 
nary hadronic matter and the QGP could be a thermodynamic phase transition of first order, of second 
order, or simply a crossover transition. 

The QGP was certainly present in the evolution of the early universe. The early universe was 
very hot, but close to net-baryon free. In the opposite limit of small temperature and large baryon 
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density, the QGP may exist even nowadays in the interior of compact stellar objects such as neutron 
stars. The main effort in present-day nuclear physics is to create the QGP under controlled conditions 
in the laboratory via collisions of heavy nuclei at ultrarelativistic energies [3] • The temperatures and 
net-baryon densities reached in nuclear collisions depend on the bombarding energy. They interpolate 
between the extreme conditions of the early universe on one side and compact stellar objects on the 
other. 

If at all, the QGP is only transiently created in a nuclear collision; free quarks and gluons will 
not be released. Therefore, deciding whether a QGP was formed or not is not easy. Detectors in the 
laboratory can only measure hadrons, leptons, or photons. The bulk of the particles emerging from a 
nuclear collision are hadrons with transverse momenta of order ~ 1 GeV. They carry information about 
the final stage of the collision after hadronization of the QGP. The formation of the latter can only 
indirectly influence this final stage, for instance by modifying the collective dynamics of the system 
through a softening of the equation of state in the hadronization transition [4 . Very few hadrons are 
emitted with transverse momenta of the order of several GeV. They arise from the fragmentation of 
jets and may carry information also about the earlier stages of the collision. Of particular interest is 
the situation where the jet has to traverse hot and dense matter possibly formed in the collision and 
is quenched by multiple rescattering in the medium jH]. From this "jet-quenching" process one may 
indirectly learn about the properties of the hot and dense environment. Finally, leptons and photons 
only interact electromagnetically. Once formed in the early stage of the collision, they leave the system 
unimpededly and carry information about this stage to the detector H . The difficulty is to disentangle 
the thermal radiation from a hot, equilibrated QGP [7j from the initial production of leptons and 
photons in the very first, highly energetic partonic collisions and from the thermal radiation of hot 
hadronic matter j^j. 

In order to decide whether a QGP was formed, one has to have detailed knowledge about its 
properties. Otherwise it is impossible to find an unambiguous signature for QGP formation in nuclear 
collisions. In this review, I present an overview of the thermodynamical properties of the QGP. Sectional 
contains a detailed discussion of the phase diagram of strongly interacting matter. The present status of 
knowledge is shown schematically in Fig. ^ Depending on the temperature, T, and the quark chemical 
potential, fx, strongly interacting matter may occur in three distinct phases: the hadronic phase, the 
QGP, and color-superconducting quark matter. The ground state of (infinite) nuclear matter is at 
(T, /i) = (0,308) MeV. There is a line of first-order phase transitions emerging from this point and 
terminating in a critical endpoint at a temperature of order ~ 10 MeV. At this point, the transition 
is of second order. This phase transition is the nuclear liquid-gas transition 9j. To the left of the line 
nuclear matter is in the gaseous phase, and to the right in the liquid phase. Above the critical endpoint, 
there is no distinction between these two phases. 

For temperatures below ~ 160 MeV and quark chemical potentials below ~ 350 MeV (corresponding 
to net-baryon densities which are a few times the ground state density of nuclear matter), strongly 
interacting matter is in the hadronic phase. Quite similar to the liquid-gas transition, there is a line 
of first-order phase transitions which separates the hadronic phase from the QGP and terminates in a 
critical endpoint where the transition is of second order. This endpoint is approximately at (T, /i) ~ 
(160,240) MeV, cf. Sec. 13.61 For smaller quark chemical potentials (smaller net-baryon densities), the 
transition becomes crossover, and there is no real distinction between hadronic matter and the QGP. 
As will be discussed in detail in Sec. the position of the critical endpoint depends on the value of the 
quark masses. Finally, at large quark chemical potential (large baryon density) and small temperature, 
quark matter becomes a color superconductor. There can be multitude of color-superconducting phases, 
depending on the symmetries of the order parameter for condensation of quark Cooper pairs. The 
discussion in Sec. 121 is qualitative and is meant to give an overview of the phase structure of strongly 
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Figure 1: The phase diagram of strongly interacting matter (schematic). 



interacting matter at nonzero temperature and quark chemical potential. The discussion in the following 
sections is both more quantitative as well as technical and focusses on the properties of the QGP and 
color-superconducting quark matter. 

The early universe evolved close to the temperature axis in the phase diagram of strongly interacting 
matter. Matter in the core of compact stellar objects, like neutron stars, is close to the quark chemical 
potential axis, at values of fi around 400 - 500 MeV. Nuclear collisions at bombarding energies around 
Ehab ~ 1 AGeV explore a region of temperatures and quark chemical potentials around (T, fi) ~ 
(70, 250) MeV. Collisions at current RHIC energies of a/s = 200 AGeV are expected to excite matter 
in a region around and above (T, /i) ~ (170, 10) MeV. Collision energies in between these two extremes 
cover the intermediate region and, in particular, may probe the critical endpoint. 

Section El presents a review of lattice QCD. After a brief introduction to the basic principles, results 
on the order of the QCD phase transition, the equation of state of strongly interacting matter, the heavy 
quark free energy, and mesonic spectral functions are collected. For technical reasons, most lattice QCD 
calculations have been done at zero quark chemical potential. An extension to nonzero values of fi is 
difficult and has been started only fairly recently. First results will also be discussed. 

Lattice QCD is in principle an exact method to solve QCD. If one had sufficiently large computer 
power, one could not only decrease the lattice spacing and increase the size of the system to come 
sufficiently close to the continuum and thermodynamic limit, one could also sample over a sufficiently 
large number of configurations to make the statistical errors arbitrarily small. However, one still has 
to interpret the results in physical terms. In this respect, analytic approaches to solve QCD have a 
certain advantage over lattice QCD. In an analytic approach, one has complete control over the physical 
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assumptions entering the calculation. Section |U gives an overview of what is known about the QGP 
from analytic calculations. 

The most simple approach from a conceptual (albeit not technical) point of view is to determine 
the thermodynamical properties of the QGP by a perturbative computation of the QCD partition 
function in terms of a power series in the strong coupling constant, g. This can be done up to terms 
of order 0(g 6 \ng). At order 0(g 6 ), the perturbative series breaks down ^21 E], and the remaining 
nonperturbative contribution has to be determined, for instance, from a lattice QCD calculation. Those 
terms of the perturbative series, which are analytically computable, are rapidly decreasing in magnitude 
at high temperatures where the strong coupling constant is small, g <C 1. This gives rise to the hope 
that only the first few terms of the perturbative series are actually necessary to obtain physically 
reasonable values for the QCD partition function. For temperatures of order ~ 150 MeV, however, 
g ~ 1 and the perturbative series is obviously not convergent. Therefore, one has tried other ways 
to compute the partition function, either by expanding around a nonperturbative ground state or by 
resumming certain classes of diagrams to improve the convergence properties of the perturbative series. 
In both approaches, quarks and gluons are considered as quasiparticles with a dispersion relation which 
is modified as compared to the one in the vacuum. Still another approach is to construct an effective 
theory for QCD which can be solved exactly or at least within a many-body approximation scheme. 
All these approaches will be reviewed in Sec. 0] 

Section contains an introduction to color superconductivity at large quark chemical potentials. 
In this case, analytic calculations are well under control, because corrections can be systematically 
computed in terms of powers of g. After a derivation of the gap equation for the color-superconducting 
gap function, the excitation spectrum in a color superconductor is presented. The solution of the 
gap equation is discussed and the magnitude of the gap parameter is determined. As in ordinary 
superconductors, quark Cooper pairs break up, if the thermal motion becomes too large. The critical 
temperature for the transition between the normal- and the superconducting phase is computed. Finally, 
the properties of gluons and photons in color superconductors are discussed. Section El concludes this 
review with a brief summary of the material and an outlook towards future directions of research in 
this area. 

A lot of the material in this review can also be found in other places. The standard review for 
properties of hot and dense, strongly interacting matter is the seminal work of Gross, Pisarski, and 
Yaffe [T2j. The contents as well as more recent developments have found their way into textbooks 
^OlE]- For early reviews focussing on the properties of the QGP, see Ref. |14j . An introduction to 
lattice QCD and recent results can be found in Ref. [TSJ. The present status of lattice QCD is also 
reviewed in Ref. JH]. Resummation techniques which attempt to compute the QCD partition function 
analytically are discussed in great detail in Ref. |TT] . 

The present review tries to give a balanced overview of all subfields concerned with the equilibrium 
properties of the QGP. Therefore, the presentation is not as detailed as in a more specialized review. 
On the other hand, I tried to explain the basics in somewhat greater detail than usually found in the 
literature. My hope is that in this way, this review will become useful for early-stage researchers working 
in both theory as well as experiment, and for all researchers who would like to get an overview of the 
theoretical activity related to equilibrium properties of the QGP. 

The only somewhat more specialized and thus more technical part is the section on color super- 
conductivity. This field has seen a lot of activity only fairly recently, but there are already a couple 
of excellent reviews jTH]- These, however, focus mainly on the basics of the phenomenon of color su- 
perconductivity and its phenomeno logical implications. In contrast, Sec. 03 contains a very detailed 
discussion of how to compute properties of the quasiparticle excitations in a color superconductor in 
the weak-coupling limit. By clarifying some of the technical details, I hope to remove the obstacles that 
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may prevent other researchers to enter this rapidly evolving and rather exciting new field of strongly 
interacting matter physics. 

Due to space restrictions this review is forced to omit many things that could (and possibly, should) 
also have been said about strongly interacting matter at high temperature and/or density. Fortunately, 
most of these have already been covered in review articles. These are, for instance, non-equilibrium 
properties of the QGP ^H] or the physics of instantons in non-Abelian gauge theories [2U|. Another 
important topic which is not mentioned in this work, but for which excellent reviews exist [21] . are 
the experimental signatures for the QGP. Recent developments in the field of color superconductivity 
are mainly focussed on deriving effective theories for quarks around the Fermi surface. These greatly 
simplify calculations and allow to systematically study effects of nonzero quark masses, for details, see 
Ref. [22]- Finally, the list of references is, necessarily, far from complete. I would like to apologize to 
all authors whose work should have been (but was not) mentioned. 

The units are h = c = kg = I. I work in Euclidean space-time at nonzero temperature T, i.e., 
space-time integrals are f x = $ T dr f v d 3 x, where V is the 3- volume of the system. Energy-momentum 
integrals are J K = TJ2 n f d 3 k/(27r) 3 ; n labels the Matsubara frequencies oj^ = 2nnT for bosons and 
uj { n = (2n + 1)ttT for fermions, n = 0, ±1, ±2, .... I denote 4-vectors with capital letters, but unless 
mentioned otherwise, retain a notation familiar from Minkowski space: X M = (t,x), where t = it, and 

= (k , k), where ko = iu n , with the metric tensor g^ u = diag(+, — , — , — ). 3-vectors k have length 
k = |k| and direction k = k/k. 



2 The QCD Phase Diagram 

2.1 Basics 

In order to understand the phase structure of strongly interacting matter, one has to know its equation 
of state. In the grand canonical ensemble, the equation of state is determined by the grand partition 
function ^ [T3| 

Z(T,V,fi) = Jv$Thj>VA£ exp [£(£ + // J\0] , (1) 

where \i is the quark chemical potential associated with (net) quark number conservation. The QCD 
Lagrangian is given by 

C = $ (yfDp - m) $ - - FST + • (2) 



For N c colors and Nf flavors, if) is the 4iV c iVf-diHiensional spinor of quark fields, ip = V^7o is the Dirac 
conjugate spinor, 7^ are the Dirac matrices and m is the quark mass matrix. The covariant derivative 
is defined as = <9 M — igA a ^T ai with the strong coupling constant g, the gluon fields A%, and the 
generators T a of the local [SU(N C ) C ] symmetry. (Throughout this paper, I indicate local, i.e. gauged, 
symmetries by square brackets.) The latter are usually taken as T a = A a /2, where A a are the Gell-Mann 
matrices. The gluonic field strength tensor is defined as 

FT = d»K - d v A» + gf abc A% A" c , (3) 

where f a bc are the structure constants of [SU(N C ) C ]. The term £ ga u g e in Eq. (J2J will not be specified 
further. It comprises gauge fixing terms and the contribution from Faddev-Popov ghosts. The number 
density operator associated with the conserved (net) quark number is Af = xjj 70 ip. 

For any finite volume V and nonzero temperature T, the partition function is defined for a compact 
Euclidean space-time volume Vxl/T. For the sake of simplicity (but without loss of generality), assume 
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that the spatial volume V is a box, V = L 3 , with L being the length of the box in one spatial dimension. 
All fields are then usually taken to be periodic in spatial directions, <f>(r, 0, y, z) = <f>(r, L,y, z), where 
stands generically for if), if), and A%. Bosonic fields, such as gluons, are periodic also in temporal 
direction, A%(0, x) = A%(l/T, x), while fermionic fields, such as quarks, are antiperiodic in temporal 
direction, -?/>(0,x) = — if)(l/T, x). 

From the grand partition function, one can derive other thermodynamic quantities, for instance the 
pressure, 

p(T,u)=T^ -> ^lnZ(y^oo). (4) 

In the thermodynamic limit, In Z is an extensive quantity (~ V) and the dependence of the pressure 
on V drops out. 

Phase transitions are determined by studying the derivatives of the pressure with respect to T and u 
for a given point (T, fi) in the phase diagram of the independent thermodynamic variables temperature 
and chemical potential. For a phase transition of first order, the first derivatives 
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are discontinuous while the pressure p is continuous at the point (T, /x). Here, s is the entropy density 
and n the (net) quark number density. For a phase transition of second order, the second derivatives 
are discontinuous, while p and its first derivatives are continuous. In this way, phase transitions of 
arbitrarily high order can be defined. One speaks of a crossover transition, if thermodynamic properties 
change rapidly within a narrow range of values T and /i, but the pressure and all its derivatives remain 
continuous. Usually, the points (T, /i) where a phase transition occurs are continuously connected and 
thus determine a line of phase transitions in the phase diagram. These lines usually start on either the 
T or the \x axis. They may terminate for nonzero values of T and \i. Two examples for this behavior, 
the liquid-gas transition in nuclear matter and the quark-hadron transition, have already been seen in 
the phase diagram of strongly interacting matter, Fig. ^ and will be discussed in more detail in the 
following. 



2.2 The Liquid- Gas Phase Transition 

The liquid-gas transition in nuclear matter is a consequence of the fact that nuclear matter assumes 
its ground state at a nonvanishing baryon density ub,o — 0.17 fm~ 3 at zero temperature T = 0. The 
underlying microscopic mechanism for this phenomenon is a competition between attractive and repul- 
sive forces among nucleons, with the attraction winning at this particular value of the baryon density. 
(This is good, because otherwise there would be no stable atomic nuclei, precluding the existence of 
our universe as we know it.) In infinite, isospin-symmetric nuclear matter, nucleons in the ground 
state are bound by —16 MeV (if one neglects the Coulomb repulsion), i.e., the energy per baryon is 
(E/Nb)o = (e/rifljo = — 16 MeV ~ 924 MeV, where e is the energy density, ub the baryon density, 
and mjv — 939 MeV is the rest mass of the nucleon. Nuclear matter is mechanically stable in the 
ground state, such that the pressure vanishes, p = 0. From the fundamental relation of thermodynam- 
ics, e = Ts + fin — p, one then concludes that the baryon chemical potential in the ground state is 
identical to the energy per baryon, hb,o = (e/^_B)o — 924 MeV. Since a baryon consists of three quarks, 
Ub = n/3 and [lb = 3/i. Hence, the ground state of nuclear matter is the point (T, /j) — (0,308) MeV 
in the nuclear matter phase diagram. 

Obviously, it costs energy to compress nuclear matter to baryon densities % > ub,q- Such an 
increase in energy leads to an increase in pressure. At zero temperature, this can be immediately 
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seen from the identity p = n 2 B <\( k E/Nb)/ ( ^Ub- Since the pressure is a monotonous function of the 
thermodynamic variables, and since it vanishes in the ground state, there are only two possibilities for 
the behavior of p for densities n B < ub,o'- either the pressure remains zero as the density decreases, 
p = 0, or the pressure further decreases such that p < 0. The latter possibility implies that the system 
becomes mechanically unstable. This can be prevented by fragmenting nuclear matter into droplets. 
These droplets are mechanically stable, i.e., the density inside each droplet is equal to n BQ and the 
pressure vanishes. The total baryon density in the system can be further reduced by simply decreasing 
the droplet density. The pressure in such a system remains zero down to arbitrarily small densities, 
because compression just results in a decrease of the space between droplets. Thus, p = from % = 
to tib = n,B,o, and then p > for ris > %,o- 

At small, but nonzero temperatures, this picture remains valid, with the additional possibility to 
thermally evaporate single nucleons from the surface of the droplets. At small temperatures and densities 
below the ground state density, one thus has a mixed phase of nucleons and droplets of nuclear matter. 
This is reminiscent of water which, at room temperature and normal pressure, consists of a mixed 
phase of water molecules and water droplets. Changing the density one can alter the relative fraction of 
molecules and droplets. Beyond the density where droplets fill the entire volume one enters the liquid 
phase, while below the density where the last droplet fragments into molecules one enters the gas phase. 
This behavior is typical for a first-order phase transition. In this case, this is the so-called liquid-gas 
transition in water. 

Nuclear matter shows a similar behavior, featuring a "gaseous" phase of nucleons at small chemical 
potentials (densities) and a "liquid" phase of nuclear matter at large chemical potentials (densities), 
cf. Fig. ^ At small temperatures the transition between the two phases is of first order. Thus, in 
the (T, /i) phase diagram there is a line of first-order phase transitions extending from the nuclear 
ground state (0, 308) MeV up towards larger values of T and smaller values of [i. As for water, this line 
terminates at a critical point where the transition becomes of second order. The critical temperature is 
of the order of 10 MeV. As for water, one cannot distinguish between the gaseous and the liquid phase 
for temperatures above this critical temperature. The existence of the liquid-gas phase transition has 
been confirmed in heavy-ion collision experiments at BEVALAC and GSI energies (i^Lab ~ 1 AGeV), 
although the precise value for the critical temperature and the critical exponents remain a matter of 
debate [Oj. 

The liquid-gas transition is also a feature of phenomenological models for the nuclear interaction, 
for instance, the Walecka model I n the following section another phase transition in strongly 

interacting matter is discussed, which very much resembles the liquid-gas transition in that it (most 
likely) is of first order for small temperatures and terminates in a critical point where the transition 
becomes of second order. This transition is the so-called quark-hadron transition. 

2.3 The Quark-Hadron Phase Transition 
2.3.1 Qualitative arguments 

For a non-interacting, translationally invariant system a convenient basis of states are the single-particle 
momentum eigenstates. Due to the Pauli principle, the density in a fermionic system can only be 
increased by successively filling states with higher momentum. The highest filled state defines the 
Fermi surface of the system, and the corresponding momentum is the Fermi momentum, kp. For non- 
interacting systems at zero temperature, the single-particle density n is given in terms of the Fermi 
momentum as 

n = J-kl, (6) 



7 



where d counts the internal degrees of freedom of the fermion (like spin, color, flavor, etc.). Thus, at 
large densities the Fermi momentum becomes large. 

In a cold, dense fermionic system particles can only scatter elastically if their momenta lie on the 
Fermi surface, as states below the Fermi surface are not accessible due to the Pauli principle (the so 
called "Pauli-blocking" phenomenon), and states above the Fermi surface are not accessible due to 
energy conservation. If the Fermi momentum exceeds the QCD scale parameter Aqcd ~ 200 MeV, 
scattering events between nucleons start to probe distances of the order 1 fm or less, i.e., the nucleonic 
substructure of quarks and gluons becomes visible. The Fermi momentum in the ground state of 
nuclear matter can be inferred from Eq. (JHJ) to be kp t o — 250 MeV. This is already of the order of Aqcd- 
Nevertheless, a description of nuclear matter in terms of nucleonic degrees of freedom is certainly feasible 
around the ground state. At which densities does a description in terms of quark and gluon degrees of 
freedom become more appropriate? The "volume" occupied by a single nucleon can be estimated from 
its charge radius to be ~ 2fm 3 . On the other hand, the specific volume of the system in the ground 
state is n^ ~ 6 fm 3 . In this sense, nuclear matter in the ground state is dilute. However, increasing the 
density to about 3n B , the system becomes densely packed with nucleons. At larger densities, they will 
even start to overlap. Therefore, around densities of a few times nuclear matter ground state density, 
one expects that a description of the system in terms of quarks and gluons is more appropriate. 

Similar arguments also apply to a system at nonzero temperature, even when the net-baryon number 
density is small. At nonzero temperature, nuclear matter consists not only of nucleons but also of other, 
thermally excited hadrons. For a non-interacting system in thermodynamical equilibrium and neglecting 
quantum statistics, the hadron number densities are proportional to n« ~ mf TK 2 {mi/T) e^ /T , where 
i labels the hadron species, m{ is their mass, fii is their chemical potential, and K2(x) is a modified 
Bessel function of the second kind. For nonzero temperature and small net-baryon number density, the 
lightest hadrons, the pions, are most abundant. At nonzero temperature and small baryon chemical 
potential, the typical momentum scale for scattering events between hadrons is set by the temperature 
T. If the temperature is on the order of or larger than Aq CD , scattering between hadrons starts to probe 
their quark-gluon substructure. Moreover, since the particle density increases with the temperature, 
the hadronic wave functions will start to overlap for large temperatures. Consequently, above a certain 
temperature one expects a description of nuclear matter in terms of quark and gluon degrees of freedom 
to be more appropriate. 

The picture which emerges from these considerations is the following: for quark chemical potentials 
\i which are on the order of 350 MeV or smaller, and for temperatures T < Aqcd ~ 200 MeV, nuclear 
matter is a gas of hadrons. (At very small temperatures T < 10 MeV, there is a gaseous and a liquid 
nucleonic phase, cf. Sec. 12.21 ) On the other hand, for T,ji^> Aqcd, nuclear matter consists of quarks 
and gluons. The natural question which emerges is, whether the "hadron phase" and the "quark-gluon 
phase" (the QGP) are separated by a phase transition in the thermodynamic sense. The rigorous way 
to decide this question is by identifying an order parameter which is nonzero in one phase and zero in 
the other. This will be discussed in more detail in the following. 



2.3.2 Pure gauge theory 



Let us first study the pure [SU(N C )} gauge theory, i.e., QCD without dynamical quarks (sometimes also 
termed the Nf = case). In this theory, there is a phase transition between a low-temperature and a 
high-temperature phase, cf. Sec. 13.21 The order parameter for this transition is the expectation value 
(L(x)) of the Polyakov loop (or Wilson line) 
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(7) 
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where P stands for path-ordering. The expectation value of an operator O in the grand canonical 
ensemble is defined as 

(O) = — Jv^V^jVA^O exp J x (£ + A*A0 . (8) 

The expectation value (L(x)) vanishes in the low-temperature phase. If one interpretes this expectation 
value as ~ exp(— Fq/T), where Fq is the free energy of an infinitely heavy quark [21], then (L(x)) = 
implies that the free energy is infinite, corresponding to confinement of color charges. In the high- 
temperature phase, (L(x)) ^ 0, which implies that the free energy of an infinitely heavy quark is finite. 
This indicates the liberation of colored degrees of freedom, i.e., deconfinement. The expectation value 
of the Polyakov loop is therefore the order parameter for the deconfinement transition. 

For an [SU(N C )\ gauge theory the action has a global Z(N C ) symmetry: the action does not change 
when multiplying all time-like links at a given spatial position x by an element z = exp(i2irn/N c ) of the 
center Z(N C ) of the gauge group [SU (N c )]. In the high-temperature phase, the nonzero expectation value 
of the Polyakov loop breaks this symmetry spontaneously. In the low-temperature phase, {L(x)) = 0, 
and this symmetry is restored. For two colors, N c = 2, the effective theory in the critical region around 
the phase transition is given by a Z(2) spin model, i.e., it is in the same universality class as the Ising 
model [21]. This model has a second-order phase transition. For N c = 3, the effective theory is that of a 
Z(3) spin model j2H]j i-e., it is in the universality class of the 3-state Potts model which has a first-order 
phase transition The transition temperature was computed to be T c ~ 270 MeV [13 EE], see also 
Sec.O 

2.3.3 Dynamical quarks 

In the presence of dynamical quarks, Nf > 0, the picture becomes somewhat more complicated. The 
fermionic term in the QCD Lagrangian (J2J) breaks the Z(N C ) symmetry explicitly, and thus there 
is strictly speaking no order parameter for deconfinement. Nevertheless, the QCD transition in the 
presence of massless dynamical quarks has an order parameter, which is related to the chiral symmetry 
of QCD. While the QCD Lagrangian (j2J) is chirally symmetric when m — 0, the ground state of QCD is 
not, i.e., chiral symmetry is spontaneously broken. It is instructive to review these arguments in more 
detail. 

In the chiral limit, where the quark mass matrix is zero, m — 0, the QCD Lagrangian (J2J) is invariant 
under global chiral U(Nf) r x U(Nf)i rotations of the quark fields. To see this, decompose the quark 
spinors into right- and left-handed spinors, 

1 ±75 

ifj = lfj r + lfj e , VV,f = V T Jl V > ^r,£ = ~ , (9) 

where V r ^ are chirality projectors. Then perform a U{N/) r ^ transformation on the right /left-handed 
quark spinors, 

/ N f-' \ 

VV U r ,efp r ,e , U r/ = exp i ^ a" £ T a eU(N f ) r/} (10) 

V a=0 / 

where a a rt are the parameters and T a the generators of U(Nf) r £. The Lagrangian (J2J) remains invariant 
under this transformation, C(if> r , if)g) = C(U r ip ri XJg ipi). The chiral group U (Nf) r x U (Nf)t is isomorphic 
to the group U{Nf)y x U(Nf)^ of unitary vector and axial transformations, where V = r + £, A = r — £, 
i.e., ay = (a r + ae)/2, a a = {oe r ~ oe)/2. Any unitary group is the direct product of a special 
unitary group and a complex phase, U(N f ) = SU(N f ) x 17(1). Thus, U(N f ) r x U(N f ) e ^ SU(N f ) r x 
SU(N f ) e x U(l) r x U{l)i ^ SU(N f ) r x SU(N f ) e x U{l)y x U{1) A . The vector subgroup U{l) v of 
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this symmetry group corresponds to quark number conservation. As physical states trivially conserve 
quark number, this subgroup does not affect the chiral dynamics and can be omitted from the further 
symmetry consideration. This leaves an SU(Nf) r x SU(Nf)g x U(1)a symmetry. The axial U(1)a 
symmetry is broken explicitly by instantons (the so-called U(1)a anomaly of QCD) [27|, leaving an 
SU(Nf) r x SU(Nf)i symmetry which determines the chiral dynamics. Since instantons are screened 
in a hot and/or dense medium the U(1)a symmetry may become effectively restored in matter. 
Then, the chiral symmetry is again enlarged to SU(Nf) r x SU(Nf)i x U{1)a- 

Nonzero quark masses break the chiral symmetry of the QCD Lagrangian explicitly. The quark 
mass term in Eq. (j2J) is 

i> % rriij ip 3 = i> % T rriij + ^\ rriij ip J r , (11) 

where flavor indices i,j = l,...,Nf are explicitly written, a sum over repeated indices is implied, 
and the properties V r ,£jo = 7o^,r, V r Vi = VeV r = of the chirality projectors were used. Now 
suppose all quark masses were equal, = mSy. Performing chiral SU(Nf) r x SU(Nf)t x U(1)a 
rotations of the quark fields, one observes that the mass term flllj) preserves an SU(Nf)y symmetry. 
All axial symmetries are explicitly broken. If less than Nf quark masses are equal, say M < Nf, the 
preserved vector symmetry is SU(M)y. In nature, where m u ~ <C m s <C m c <C <C m t , one 
only has the well-known (approximate) SU(2)y isospin symmetry. Consequently, exotic hadrons with 
strange, charm, bottom, or top degrees of freedom are not degenerate in mass with their non-strange 
counterparts. 

The mass term ip l rriij i n ^ ne QCD Lagrangian is of the same form as the term H - S in spin models, 
which couples the spin S to an external magnetic field H. Obviously, the operator %p l ipi corresponds to 
the spin S, while the quark mass matrix assumes the role of the external magnetic field H. Thus, 
the expectation value %pi) is the analogue of the expectation value of the spin, the magnetization 
M = (S). While the mass term explicitly breaks the chiral symmetry SU{Nf) r x SU(Nf)t x U(1)a of 
the QCD Lagrangian to an (approximate) SU(2)y symmetry, the external magnetic field in spin models 
explicitly breaks the rotational symmetry 0(3) to 0(2). 

The analogy between QCD and spin models, however, extends further than this. Even in the 
absence of external magnetic fields, in spin models with ferromagnetic interactions rotational symmetry 
is spontaneously broken due to a nonvanishing magnetization M ^ in the ferromagnetic phase. Ana- 
logously, in the QCD vacuum, chiral symmetry is spontaneously broken by a nonvanishing expectation 
value (jj) 1 V> J )vac. 7^ 0. Let us introduce the so-called chiral condensate $ y ' and its complex conjugate, 
$ ij , via 

A nonvanishing expectation value ffiifri) ^ is then equivalent to $ lJ + $ y ^ 0. Just like the mass 
term in the QCD Lagrangian, a nonvanishing chiral condensate breaks the chiral symmetry. In the chiral 
limit, rriij = 0, nothing distinguishes one quark flavor from another and, consequently, = 0o ^ ■ ( in 
principle, there is another possibility how a chiral condensate could break the chiral symmetry, for a more 
detailed discussion see below and Ref. 29-) This chiral condensate breaks the chiral U(Nf) r x U(Nf)i 
symmetry spontaneously to U(Nf)y. To see this, note that the chiral condensate is still invariant under 
vector transformations, $ — > U r $ U\ = if U r — Ui = Uy, but not under axial transformations, 
$ -»• U r $ U\ ^ $, if U r = U\ = U a- 

According to Goldstone's theorem, the breaking of the global chiral symmetry leads to the occurrence 
of massless modes, the so-called Goldstone bosons. The number of Goldstone bosons is identical to the 
number of broken generators. In the QCD vacuum, where the U(1)a anomaly is present, the breaking 
pattern is SU(Nf) r x SU(Nf)g —> SU(Nf)y, i.e., in this case there are Nj — 1 broken generators, 
corresponding to the generators of the broken axial symmetry SU(Nf)A- For Nf = 1, there is no 
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global chiral symmetry that could be broken, and thus no Goldstone boson. For Nf = 2, the Goldstone 
bosons are the three pions, the lightest hadronic species. In nature, the pions are not completely 
massless, because the chiral symmetry is explicitly broken by the small, but nonzero quark mass term 
in the QCD Lagrangian. This turns the Goldstone bosons into so-called pseudo-Goldstone bosons. For 
Nf = 3, the pseudo-Goldstone bosons correspond to the pseudoscalar meson octet, comprising pions, 
kaons, and the eta meson. Since chiral symmetry is more strongly broken by the larger strange quark 
mass, the pseudo-Goldstone bosons carrying strangeness are heavier than the pions. For Nf > 4, the 
explicit symmetry breaking by the heavy exotic quark flavors is so strong that the would-be Goldstone 
bosons are actually heavier than the ordinary (i.e. non- Goldstone) non-strange bosons. 

In spin models, rotational symmetry is restored above some critical temperature and the magnetiza- 
tion vanishes. The magnetization is the order parameter for this so-called ferromagnet-diamagnet phase 
transition. By analogy, one expects a similar mechanism to occur in QCD, i.e., to vanish above some 
critical temperature. The symmetry of the ground state is then restored to the original chiral symmetry, 
i.e., SU(Nf) r x SU(N f ) e , if the U(1) A anomaly is still present, or SU(N f ) r x SU(N f ) e xU(l) A , if instan- 
tons are sufficiently screened at the transition in order to effectively restore the U(1)a symmetry. Lattice 
QCD calculations show that this expectation is indeed fulfilled: there is a phase transition between the 
low-temperature phase where chiral symmetry is broken and the high-temperature phase where it is 
restored, for details see Sec. 0J Just like the magnetization in spin models, the chiral condensate $ y ' is 
the order parameter for this so-called chiral phase transition. 

In the case of massless quarks, = 0, based on universality arguments one can analyse the order 
of the chiral transition in the framework of a linear sigma model for the order parameter field $* J [29J. 
This linear sigma model is an effective theory, i.e, all terms allowed by the original chiral symmetry 
must in principle appear, 

C eS = Tr (d &d°$) - v 2 Tr (v$ f ■ V$) - V eS ($) , (13) 
where the effective potential 

V eS (<$>) = m 2 Tr + Ai [Tr 2 + A 2 Tr ($t$) 2 - c (det$ + det$ f ) + . . . (14) 



determines the ground state of the theory. In Eq. ()13j) it was assumed that the first term is canonically 
normalized. However, since Lorentz symmetry is explicitly broken in a medium at nonzero temperature, 
the coefficient v 2 in Eq. (|13|) may in general be different from one. In Eq. ()14j) . . . . denote higher- 
dimensional operators which are irrelevant for the discussion of the order of the phase transition. For c ^ 
0, the chiral symmetry of C cS is SU(N f ) r x SU(N f ) e , while for c = 0, it is SU(N f ) r x SU(N f ) e x U{l) A . 
Thus, the U(1) A anomaly is present for c ^ 0, and absent for c = 0. While these chiral symmetries 
are manifest in the Lagrangian (|13|). the ground state of the theory respects them only for c = and 
m 2 > 0. For c = and m 2 < the chiral symmetry is spontaneously broken by a nonvanishing vacuum 
expectation value for the order parameter. Consequently, if the linear sigma model is to describe the 
chiral transition in QCD, one has to ensure that m 2 < for c = 0. There are still two possibilities how 
the order parameter can break the symmetry. As shown in Ref . [2H] , if A 2 > the ground state is given 
by $^ c . = 0o^j while for A 2 < the ground state is given by $^ ac = 4>q5 %1 8^ 1 . (The choice of the 
1-direction in right- and left-handed flavor space is arbitrary.) Nature realizes the first possibility. For 
the case c 7^ 0, no general arguments can be made; whether the ground state of the theory breaks chiral 
symmetry spontaneously depends on the particular values for the coupling constants c, Ai, A 2 and the 
number of flavors Nf. 

For Nf = 1, there is no difference between the two quartic invariants in Eq. (|14|) . and one may 
set Ai + A 2 = A. In the presence of the U(1) A anomaly, c 7^ 0, there is also no chiral symmetry, and 
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the transition is crossover, due to the linear term ~ c, which tilts the effective potential such that the 
(thermal) ground state ($)t 7^ 0. If the U(1)a anomaly is absent, c = 0, the effective theory for the 
order parameter falls in the same universality class as that of the 0(2) Heisenberg magnet, and thus 
the transition is of second order. 

For Nf = 2 and in the presence of the U(1)a anomaly, the chiral symmetry is SU(2) r x SU(2)£, 
which is isomorphic to 0(4). The effective theory for the order parameter is in the universality class 
of the 0(4) Heisenberg magnet. Consequently, the transition is of second order [2H]. If the U(1)a 
symmetry is effectively restored at the phase transition temperature, the symmetry group is larger, 
SU{2) r x SU(2)e x U(1)a, which is isomorphic to 0(4) x 0(2), and the transition is of first order. 
Lattice QCD calculations determine the transition temperature to be T c ~ 172 MeV [T5| I16j. 

For Nf = 3, the chiral transition is of first order, both when the U(1)a symmetry is explicitly broken 
by instantons or when it is effectively restored at the transition. In the first case, the effective theory 
features a cubic invariant in the order parameter field (the term ~ det$ + det$^), which drives the 
chiral transition first order In the second case, the transition is fluctuation-induced of first order 
[29j. This also holds for Nf > 4, irrespective of whether the U(1)a symmetry is explicitly broken or 
not. For Nf = 3, lattice QCD calculations find the transition temperature to be T c ~ 155 MeV |16j . 
cf. Sec. 13.21 Note that nonvanishing quark masses can also be accounted for by adding a term 



to the right-hand side of Eq. (|13|). As discussed above, this term is the analogue of the term H • S in 
spin models. Consequently, the external "magnetic field" is proportional to the quark mass matrix 



2.3.4 The quark-mass diagram 

Nonvanishing quark masses lead to the term (|15|) in the effective theory for the order parameter field. 
This term is linear in $, such that the effective potential is tilted. This may render a first or second 
order phase transition a crossover transition (similar to the case Nf = 1 with U(1)a anomaly discussed 
in Sec. 12.3.31 where a tilt in the potential is induced by the linear term ~ c). For instance, the second- 
order transition for QCD with Nf = 2 massless flavors is rendered crossover by a nonvanishing quark 
mass. The first-order phase transition for QCD with Nf = 3 massless flavors can also become crossover, 
if the quark masses are sufficiently large. In the real world, the up and down quark are approximately 
of the same mass, while the strange quark is much heavier. It is customary to put m q = m u ~ 
and identify first-order regions, second-order lines, and crossover regions in an (m q ,m s ) diagram, see 
Fig. |21 To simplify the following discussion, only the case where the U(1)a anomaly is present will be 
considered. 

The origin in Fig. ^corresponds to the massless 3-flavor case, and the transition is of first order. The 
upper left corner corresponds to the massless 2-flavor case, since the strange quark is infinitely heavy. 
Here, the transition is of second order. The lower right corner is the case of one massless flavor. The 
transition is crossover. The upper right corner, where all quark flavors are infinitely heavy, corresponds 
to the pure gauge theory. At this point the transition is of first order. 

The first-order regions around the origin and the upper right corner extend somewhat into the 
(m q ,m s ) plane and are bounded by critical lines where the transition is of second order. Along these 
critical lines, the second-order phase transitions are in the universality class of the Ising model, Z(2). In 
between the critical lines, the transition is crossover. There is also a second-order phase transition line 
(with a phase transition in the 0(4) universality class) extending downwards from the upper left corner 
along the m s axis. There is a tricritical point where this line meets the second-order phase transition 
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Figure 2: The quark- mass diagram (from Ref. [TK]). 



line bordering the first-order region around the origin [3U|. It is an interesting question whether the 
real world, where m q ~ 5 MeV -C m s ~ 100 MeV, is still inside the first-order phase transition region or 
already in the crossover region. There are ongoing lattice QCD studies to decide this question, which 
at present favor the latter possibility [TB] , see also Sec. 13.21 for more details. 



2.3.5 Nonzero quark chemical potential 

So far, the quark-hadron phase transition was studied at /z = 0. Let us finally discuss the case of nonzero 
quark chemical potential. For many years, lattice QCD studies at nonzero chemical potential were ham- 
pered by numerical problems related to the so-called sign problem of the fermion determinant. Only 
recently there have been attempts to compute the order of the phase transition, as well as thermodynam- 
ical properties, at nonzero quark chemical potential; for details, see Sec. 13.61 So far, these calculations 
have been done on fairly small lattices with rather heavy quarks. Consequently, they show a crossover 
transition at jj, — 0. This crossover transition extends to the point (T, /i) cr = (160 ± 3.5, 242 ± 12) MeV, 
see also Fig. [TUJ 

This point is a critical point where the transition is of second order. It is in the unversality class 
of the Ising model, i.e., Z(2). For smaller temperatures and larger chemical potentials, the transition 
becomes of first order. The critical point will move towards the T axis when the quark masses are 
decreased. From the discussion in Sec. 12.3.41 one cannot exclude the possibility that, for realistic quark 
masses, the first-order phase transition line emerges directly from the T axis. 

Finally, the question arises whether the line of first-order phase transitions extends all the way down 
to T = 0, and if so, at which point it hits the /i axis. Renormalization group arguments [31] suggest 
that the behavior at zero temperature is very similar to the one at nonzero temperature: the transition 
is of first order for Nf > 3 as well as for Nf = 2 in the absence of the U(1)a anomaly, while it could 
be of second order for Nf = 2 in the presence of the U{1)a anomaly. In the latter case, however, it 
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would remain of second order along the whole phase transition line, and only the universality class 
would change from 0(4) critical behavior at nonzero temperature to Gaussian critical behavior at zero 
temperature. Since lattice QCD calculations [32] indicate that the transition becomes of first order for 
temperatures below the critical point, and since it is hard to imagine that the transition switches back 
to second order as the temperature decreases further, this possibility can most probably be ruled out. 
Note, however, that if quark-gluon matter is in a color-superconducting phase to the right of the QCD 
phase transition line, other possibilities emerge, for details, see Sec. 12.41 

In any case, model calculations [33] within a Nambu-Jona-Lasinio (NJL) model 34 support the 
picture that the transition remains of first order below the critical point all the way down to the \i axis. 
The value of /i, where the first-order phase transition line meets the \x axis, depends sensitively on the 
parameters used in these model calculations. Its actual value should not be taken too seriously, because 
the NJL model with quark degrees of freedom does not have a phase where matter consists of nucleons. 
Instead, the transition to quark matter happens at a quark chemical potential which is of the order of 
the ground state quark chemical potential, /x — 308 MeV. Since we know that nucleonic matter exists, 
this behavior is clearly unphysical. 

The critical point has recently received a fair amount of attention [35J. For a second-order phase 
transition in the Z(2) universality class, there must be one massless degree of freedom. The fact that 
this degree of freedom is massless causes critical fluctuations. These fluctuations were suggested to be 
an experimental signature for the critical point in heavy-ion collisions. Which physical particle does 
the massless degree of freedom correspond to? For realistic quark masses, the pions are not massless 
in the vacuum, and it is unlikely that they become massless at the critical point (the pions usually 
get more massive when the temperature is increased). Moreover, since isospin is still a good symmetry 
at the critical point, all pions would simultaneously become massless. Then, one would have three 
massless modes instead of just one. Consequently, the pions cannot assume the role of the massless 
mode. Since the critical point exists even when considering only Nf = 2 quark flavors and since, for 
Nf = 2, there is only one other degree of freedom in the effective theory besides the pions, it must be 
this degree of freedom that becomes massless: the scalar a meson [331 EE3 ■ in the vacuum, this meson 
has a mass of about 400 — 1200 MeV [3H]. If it becomes massless, it can be copiously produced. When 
these a mesons decouple from the collision region, they assume their vacuum masses and rapidly decay 
into pions. Besides critical fluctuations, another signature for the critical point would thus be the late 
emission of a large amount of pions in a heavy-ion collision. 

2.4 The Color- Superconducting Phases of QCD 
2.4.1 Proof of existence of color superconductivity 

There are other phases in the phase diagram of nuclear matter, which have recently received much 
attention in the literature, the so-called color-superconducting phases in sufficiently cold and dense 
quark matter [TH]. Color superconductivity occurs, because there is an attractive interaction between 
two quarks at the Fermi surface [37| EE] ■ Then, by Cooper's theorem, these quarks form Cooper pairs 
which condense in the new ground state of the system. 

At least at asymptotically large quark chemical potentials, the existence of an attractive interaction 
between quarks at the Fermi surface, and thus the existence of color superconductivity, can be rigorously 
proven. Due to asymptotic freedom 1], when /x ^> Aq C d, the strong coupling constant of QCD, 
evaluated at the scale //, becomes small, g(fx) <C 1, such that the dominant interaction between quarks 
is given by single-gluon exchange. The scattering amplitude for single-gluon exchange in an [SU(N C ) C ] 
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gauge theory is proportional to 



[Tajki (T a )lj 
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4iV c 



($jk $il + 8ik 8jl) , 



(16) 



where i,j are the fundamental colors of the two quarks in the incoming channel, and k, I their respective 
colors in the outgoing channel. Under the exchange of the color indices of either the incoming or 
the outgoing quarks the first term is antisymmetric, while the second term is symmetric. In group 
theoretical language, for [,S77(3) C ] Eq. (fTo]) represents the coupling of two fundamental color triplets to 
an (antisymmetric) color antitriplet and a (symmetric) color sextet, 



The minus sign in front of the antitriplet contribution in Eq. (|16|) signifies the fact that this channel 
is attractive, while the sextet channel is repulsive. Therefore, one expects that quark Cooper pairs 
condense in the color-antitriplet channel. 

This argument holds rigorously at asymptotically large densities. The highest densities of nuclear 
matter that can be achieved in the laboratory through heavy-ion collisions, or that occur in nature 
in the interior of neutron stars, are of the order of ten times nuclear matter ground state density. At 
these densities, the quark chemical potential is still fairly small, /J ~ 0.5 GeV. For phenomenology it is 
therefore important to answer the question whether color superconductivity also exists at the (compar- 
atively moderate) densities occurring in nature. There is no rigorous way to answer this question, as 
an extrapolation of the above asymptotic argument becomes unreliable when g(fi) ~ 1. Nevertheless, 
calculations in the framework of the NJL model show that color superconductivity does seem to 
occur also at moderate densities. In this case, the attractive interaction could be mediated by instanton 
(instead of single-gluon) exchange. 

Color superconductivity is a much more complicated phenomenon than ordinary superconductivity. 
From a very qualitative point of view, in comparison to electrons, quarks carry additional quantum 
numbers such as color and flavor. The wave function of a Cooper pair has to be antisymmetric under 
the exchange of the two fermions forming the pair. Consequently, the possible color and flavor repre- 
sentations of the two-fermion state have to be chosen in a way which respects this antisymmetry. This 
requirement helps to classify all possible color-superconducting condensates ^OJ E] • This classification 
will be presented in the following. 

2.4.2 Classification of color- superconducting phases 

The baryon density in heavy-ion collisions or in neutron stars is at most of the order of ten times 
the nuclear ground state density. Therefore, the quark chemical potential is unlikely to assume values 
beyond \i ~ 1 GeV. At zero temperature, one can only have fermions if the Fermi energy /x exceeds 
their mass. Thus, only the three lightest quark flavors, up, down, and strange, play any role; charm, 
bottom, and top quarks are too heavy. For small temperatures, these heavy flavors can be thermally 
excited, but their abundance is exponentially suppressed ~ exp(— rrif/T), cf. Sec. I2.3.H Therefore, they 
will be excluded from the following consideration. For fi ~ 1 GeV, up and down quarks can be treated 
as truly ultrarelativistic particles, as m q /fi ~ 1CT 3 . To first approximation, also the strange quark will 
be considered to be massless. Corrections due to the strange quark mass can be treated perturbatively, 
the correction factor being of order m s /{i ~ 10 _1 [22J. 

For ultrarelativistic particles, spin S and angular momentum L are not separately good quantum 
numbers, only the total spin J = L + S is. Therefore, possible Cooper pair wave functions should be 
classified according to their total spin J. Let us first focus on the spin-zero channel, J = 0. The J = 



[3] c x [3] c = [3} c a + [6} c s . 



(17) 
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representation of the SU(2)j spin group is totally antisymmetric. Therefore, the remaining color and 
flavor part of the Cooper pair wave function has to be symmetric under the simultaneous exchange 
of color and flavor indices in order to fulfill the requirement of overall antisymmetry. If one assumes 
that quarks pair in the antisymmetric color-antitriplet channel, one has no choice but to also chose 
an antisymmetric flavor representation. For Nf = 1 flavor, this is impossible, as there is no flavor 
symmetry group with an antisymmetric representation. 

One therefore has to consider at least Nf = 2 quark flavors (for instance, up and down), where the 
most simple representation is the antisymmetric flavor singlet [1]{ representation of the SU(2)y flavor 
group. Therefore, the most simple J = quark Cooper pair condensate has the form 

®{? = e ijk ef°® k . (18) 

Here, i, j = 1, . . . , N c are the color indices of the quarks forming the Cooper pair, while f,g= 1, . . . , Nf 
are the corresponding flavor indices. The two totally antisymmetric tensors on the right-hand side 
ensure that the condensate belongs to the [3]„ in color, as well as the representation in flavor space. 
The color-superconducting phase represented by the condensate (|T%|) is commonly called the "2SC" 
phase (for "2-flavor color SuperConductor"). 

Condensation of quark Cooper pairs occurs if the quantity $fc on the right-hand side of Eq. (fTHj) is 
nonzero, ^ 0. Thus, the quark Cooper pair condensate carries a fundamental color index k. This 
indicates that the local [SU(3) C ] color symmetry is spontaneously broken by the quark Cooper pair 
condensate, similarly to the spontaneous breaking of the global chiral symmetry SU(Nf) r x SU(Nf)i 
by the chiral condensate (J 12)) in the QCD vacuum discussed in Sec. 12.3.31 In this sense, <3>fc is the order 
parameter for color superconductivity. It is zero in the phase of unbroken [S77(3) c ] symmetry, and 
nonzero in the broken phase where condensation of quark Cooper pairs occurs. 

Of course, a local symmetry can never be truly broken spontaneously [12] • However, after fixing 
the gauge, spontaneous breaking does occur, just like in ordinary superconductivity or in the standard 
model of electroweak interactions. In ordinary superconductivity, the condensation of electron Cooper 
pairs breaks the [C/(l) em ] gauge symmetry of electromagnetism, while in the standard model, the Higgs 
field breaks [SU{2)(\ x [{7(l)y] to [{7(l) em ]- Analogously, the quark Cooper pair condensate (fTHj) breaks 
the [S77(3) c ] color gauge symmetry. 

Since quarks carry baryon and electric charge, the Cooper pair condensate (|18|) in principle also 
breaks the global U(l)y of baryon number conservation and the local [{7(1 ) em ] of electromagnetism. 
In the discussion of chiral symmetry breaking, these symmetries were never broken because the chiral 
condensate consists of a quark and an antiquark which carry opposite baryon and electric charge. The 
chiral condensate is thus a singlet under U(l)y and [{7(1 ) em ] and consequently preserves these symme- 
tries. This is different for a color-superconducting condensate which consists of two quarks. It turns out, 
however, that there exists a "rotated" baryon number {7(1) symmetry and a "rotated" electromagnetic 
[{7(1)] symmetry, which are formed from the original baryon number and electromagnetic symmetries 
and the eighth generator of [5*{7(3) c ] [13]. This is similar to electroweak symmetry breaking, where the 
[{7(l) em ] symmetry of electromagnetism is a combination of the [{7(l)y] hypercharge symmetry and the 
third generator of [SU(2)i]. Thus, the U(l)y and [{7(l) em ] symmetries are not really broken in the 2SC 
phase, but "rotated" . The rotation angle is the analogue of the Weinberg angle in the standard model 
of electroweak interactions; for more details, see Sec. 15.41 

Spontaneous symmetry breaking in gauge theories does not lead to Goldstone bosons. Rather, 
what would have been a Goldstone mode will be "eaten" by a gauge boson which in turn becomes 
massive and thus acquires an additional longitudinal degree of freedom. There are as many massive 
gauge bosons as there would have been Goldstone modes due to spontaneous symmetry breaking. 
In ordinary superconductors, the electromagnetic [{7(l) em ] symmetry is broken, which has only one 
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generator. Consequently, there is a single Goldstone mode which is "eaten" by the single gauge boson 
present in this case, the photon. The photon acquires a so-called Meissner mass. What happens 
physically is that magnetic fields are damped on length scales of the order of the inverse Meissner mass, 
which in turn leads to the Meissner effect, the expulsion of magnetic flux from the superconductor. In 
the standard model of electroweak interactions, [SU{2)t\ x [?7(l)y] is broken to [C/(l) em ], i.e., there are 
three Goldstone modes which in turn lead to the massive gauge bosons of the weak interaction, W ± 
and Z. The photon is massless. This is required, since it is the gauge boson of the residual [6 r (l) em ] 
symmetry of electromagnetism. Analogously, in a color superconductor one expects some of the gluons 
to become massive. Exactly how many gluons acquire a mass depends on the pattern of symmetry 
breaking. For the condensate ()18j) . one can clarify this via the following argument. 

By a global color rotation, one can always orient the order parameter to point in the 3-direction 
in color space (more precisely, the antz-3-direction, as condensation occurs in the color arziz-triplet 
channel) , 

= (19) 

Physically, this means that if we call color 1 red, color 2 green, and color 3 blue, red up (or down) quarks 
and green down (or up) quarks condense to form an anti-blue Cooper pair condensate. Blue up and 
down quarks do not participate in condensation. The condensate ([!§]) does not break the color [S77(3) c ] 
gauge symmetry completely. The residual symmetry is a local [SU(2) C ] symmetry in the space of the 
first two colors (in our conventions, red and green). Including electromagnetism, the symmetry breaking 
pattern for the condensate (fT9Jl is therefore [577 (3) c ] x [t/(l) cm ] — ► [SU(2) C ] x [^(1)]- Consequently, 
there are 8 — 3 = 5 broken generators, which lead to five massive gluons. The remaining three gluons 
must remain massless as they correspond to the gauge bosons of the residual local [SU(2) C ] symmetry. 
This is also borne out by an explicit calculation of the gluon Meissner masses in the 2SC phase |4*4^ I45j . 
The gauge boson of the local [C/(l)] symmetry (the "rotated" photon) is also massless. For more details, 
see Sec. 15.41 

For Nf = 3 flavors, condensation of quark Cooper pairs becomes considerably more interesting. 
First, to preserve the antisymmetry of the Cooper pair wave function the two quarks have to be in 
the [3]l representation of the global SU(3)v flavor symmetry. Consequently, the quark Cooper pair 
condensate has the form 

*(f = (20) 

The difference to Eq. (jl8|) is that, to ensure antisymmetry in flavor space, one is required to use 
the totally antisymmetric tensor of rank 3, e^ 9h , rather than its rank-2 counterpart. Consequently, 
an additional flavor index h appears in the order parameter, A nonvanishing order parameter 
automatically implies that not only the local [S77(3) c ] color, but also the global SU(3)v flavor symmetry 
is broken. The situation is not unlike the one encountered in superfluid helium-3 [46J. Superfluid 
helium-3 forms Cooper pairs with spin S = 1 and angular momentum L = 1. (Both spin and angular 
momentum are good quantum numbers, as helium-3 is a non-relativistic system.) Consequently, the 
order parameter breaks the global SO(3)s of spin as well as the global SO(3)l of angular momentum. 
This breaking can occur in many possible ways, giving rise to a plethora of phases in superfluid helium-3. 

Similarly, one would expect many different phases to occur in a 3-flavor color superconductor. 
However, in fact there are only two possibilities, one of which is likely to be realized in nature j41j . 
To see this, note the formal similarity between the order parameter and the one encountered in 
chiral symmetry breaking, in Sec. 12.3.31 While $ y transforms under SU(Nf) r x SU(Nf)e (in the 
presence of the U(1)a anomaly), transforms under [S77(3) c ] x SU(3)y- Consequently, the effective 
Lagrangian for is of the same form (fT3*|) as for The two possible patterns of symmetry breaking 
occurring in such an effective theory were already mentioned above. 
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If the coupling constant A2 > 0, the order parameter assumes the form 



(21) 



In contrast to the 2SC case, where blue quarks remained unpaired, now all quark colors and flavors 
participate in the pairing process. The order parameter (|21jl is similar to the one for chiral symmetry 
breaking, where = 5^ $. (In the ground state, $ V ac. = <j>o-) Similar to the chiral symmetry breaking 
pattern SU(N f ) r x SU(N f ) t -> SU(N f ) v , V = r + £, the condensate dUJ) breaks [SU(3) C ] x SU(3) V to 
the vectorial subgroup SU(3) c+ y. The condensate is still invariant under vector transformations in color 
and flavor space, or in other words, a transformation in color requires a simultaneous transformation 
in flavor to preserve the invariance of the condensate. Therefore, the discoverers [3?j of this 3-flavor 
color-superconducting state termed it the "Color-Flavor-Locked", or short, "CFL" state. The residual 
SU(3) c+ v symmetry is no longer a local symmetry, so that there are no massless gluons. The complete 
pattern of symmetry breaking, including the U(l)v symmetry of baryon number and the [{7(1 ) em ] 
symmetry of electromagnetism is [SU(3) C ] x SU(3) V x U(l) v x [U(l) cm ] -> SU(3) C+V x [{7(1)]. (This 
notation is slightly ambiguous: the [{7(l) em ] symmetry is generated by the quark charge operator 
Q = diag(2/3, — 1/3, — 1/3) which is traceless. Thus, the global part of the [{7(l) em ] symmetry is 
actually a subgroup of the global SU(3)v flavor symmetry. However, the local part of [{7(1 ) em ] is not. 
Therefore, here and in the following I choose to explicitly denote the [{7(1 ) em ] symmetry group.) In the 
CFL case, unlike the 2SC case, baryon number is broken, but a rotated electromagnetic [{7(1)] is again 
preserved. The symmetry breaking pattern leads to nine Goldstone bosons, eight of which are "eaten" 
by the gluons, i.e., all gluons acquire a Meissner mass. There is one Goldstone boson from the breaking 
of the U(l)y symmetry. 

In the chiral limit, the flavor symmetry of QCD is actually not just SU(3)y but SU(3) r x SU(3)g. 
Assuming that also the U(1)a symmetry of QCD is effectively restored at large quark densities, quark 
Cooper pair condensation of the form j2TJ) breaks [SU(3) C ] x SU(3) r x SU(3) e x U(l) v x U(1) A x [U(l) em ] 
to SU(3) c+ v x [{/(l)], i.e., not only the local color, but also the global chiral symmetry is broken. In 
addition to the eight massive gluons, there are also ten real Goldstone bosons, eight from the breaking 
of the SU(3)a chiral symmetry, and one each from the breaking of U(l)y and U(1) A - 

Closer inspection j^H] shows that the excitation spectrum in the CFL state bears a striking resem- 
blance to the one in the hadronic phase. Let us first focus on the fermionic sector. In the CFL phase, 
there are nine gapped fermionic quasiparticles (cf. also Sec. I5.2J1 . eight of which are degenerate in mass. 
These correspond to the baryon octet in the QCD vacuum. (For this argument we have to assume that 
there is no explicit SU(3)y flavor-symmetry breaking in the QCD vacuum.) The ninth quasiparticle is 
twice as heavy and does not have a counterpart in hadronic matter, but this is not a reason to worry, as 
such a particle would have a large decay width into lighter particles. In the bosonic sector, there are nine 
Goldstone bosons from the breaking of the axial SU(3) A x U(1) A symmetry, which correspond to the 
pseudoscalar nonet in the hadronic phase. Only the tenth Goldstone boson from the breaking of U(l)y 
does not have a counterpart in the QCD vacuum. Such a boson exists, however, in dense nuclear matter, 
where a superfluid A A condensate may form, which also breaks the {7(1) y baryon number symmetry. 
The preceding arguments have led to the conjecture of "continuity" between hadron and quark matter 
|48j . This conjecture states that, since there is no difference in symmetry between quark matter in 
the CFL state and (SU(3)y flavor-symmetric) hadronic matter, there need not be any phase boundary 
between these two phases at all. Of course, this requires that there is no other color-superconducting 
phase, for instance the 2SC state, which separates CFL matter from hadronic matter. 

For A2 < 0, the order parameter is given by 



(22) 
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where the 3-direction is arbitrary. This condensate breaks [51/(3) c ] x SU(3)v to [SU(2) C ] x SU(2)v- In 
this case, [SU(2) C ] is still a local symmetry. Consequently, like in the 2SC case only five gluons become 
massive, and blue up, down, and strange quarks do not participate in the formation of Cooper pairs. 
Additionally, there are also five Goldstone bosons from the breaking of the flavor symmetry. There is 
a more technical and a more physical argument, why the CFL state is most likely realized in nature. 
From the more technical point of view, one can show |H] that, to one-loop order, in QCD Ai = and 
A2 > 0. From the more physical point of view, the CFL state is energetically favored because all quark 
colors and flavors (instead of just a few) acquire a gap at the Fermi surface. The gain in condensation 
energy is thus expected to be larger than for a state with a condensate of the form (j22J)- 

Although a single quark flavor cannot form Cooper pairs with total spin J = 0, it can pair in the 
J = 1 channel. (An exhaustive discussion of possible pairing channels for a single quark flavor is given 
in Ref. jH|-) This channel corresponds to the symmetric [3]f representation of the SU{2)j spin group. 
If one still assumes pairing to occur in the color [3]° channel, the Cooper pair wave function is, as 
required, overall antisymmetric. The condensate is a 3-vector in space which points in the direction of 
the spin of the Cooper pair. It has the form |3B1 HOI EH EJ] 

% = a* n , (23) 

where a = x,y, z denotes the spatial component of the spin vector. Condensation breaks the local color 
[5£/(3) c ] symmetry and the global 50(3) j spin symmetry. This is similar to superfluid helium-3 where 
condensation breaks 50(3)s x 50(3)^ (see discussion above). While many different phases arise, let us 
just mention two which are quite similar to the ones discussed in the context of three and two quark 
flavors, the so-called "Color-Spin-Locked" or CSL phase, where the order parameter assumes the form 

n = 81 $ , (24) 

and the so-called polar phase, where 

^ a k = 5 k3 5 az ^. (25) 

In the CSL phase, the order parameter (|24jl is strikingly similar to the one in the CFL phase, cf. Eq. (|21|). 
All quark colors participate in the formation of Cooper pairs. Also the symmetry breaking pattern is 
similar, [5f/(3) c ] x 50(3) j x U{l)y x [£7(l) em ] — > 50(3) c+ j. The main difference is that now there 
is no rotated electromagnetism [E/(l)], cf. Sec. 15.41 Consequently, all eight gluons and the photon 
become massive |S2]- I n the polar phase, the order parameter resembles that of Eq. (|22j) . Neglecting 
electromagnetic interactions for the moment, the symmetry breaking pattern is [5C/(3) C ] x 50(3)j x 
U(l) v -» [SU(2) C ] x 50(2) j x U(l). Like in the 2SC phase, the residual [SU(2) C ] color symmetry is a 
local symmetry, and there are three massless and five massive gluons. The breaking of the rotational 
50(3) j symmetry to 50(2)j also leads to two real Goldstone bosons. Baryon number is not broken, but 
merely rotated. Including the [£/(l) em ] of electromagnetism there is a small subtlety which is explained 
in more detail in Sec. 15.41 If there is only a single flavor present, or if all flavors carry the same electric 
charge, a rotated electromagnetic [C/(l)] symmetry exists. If there are at least two flavors which differ 
in charge, the [t/(l) Gm ] symmetry is broken. Table ^ summarizes the results of this section for the 2SC, 
the CFL, the CSL, and the polar phase of color-superconducting quark matter. 

2.4.3 Color-superconducting phases in the nuclear matter phase diagram 

How does color-superconductivity affect the phase diagram of nuclear matter? Let us first assume 
that the temperature is sufficiently small to favor a color-superconducting over the normal-conducting 
state. As long as /i > m s , the CFL state is likely to be the ground state of quark matter. Since one 
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Table 1: Color-superconducting phases in dense quark matter. For the polar phase, there is 
an additional [C/(l)] symmetry, if all flavors in the system carry the same electric charge. 



phase 


condensate 


order parameter 


residual symmetry 


2SC 


${; = <„,.. < J:i <i'/. 


$fc = 5 k3 $ 


[SU(2) C ] x U(l) x [tf(l)] 


CFL 




= 5% $ 


S*7(3) c+ y X [[7(1)] 


CSL 






SO(3) c+J 


polar 




= 5*3 5°*$ 


[SU(2) C ] x LT(1) 



has (approximately) equal numbers of up, down, and strange quarks of colors red, green, and blue, 
the system is (approximately) neutral with respect to color and electric charge. However, when one 
extrapolates down to smaller quark chemical potentials, say of the order of [i ~ 500 MeV, the strange 
quark mass is no longer negligibly small and causes, for a given /i, a mismatch in the Fermi surfaces 
between non-strange and strange quarks jSH]- In general, a nonzero strange quark mass reduces the 
number of strange quarks as compared to the massless species. This, in turn, leads to nonzero electric 
and color charge in the system. Consequently, one is forced to introduce chemical potentials for electric 
and color charge, which have to be tuned to again ensure overall electric and color neutrality. The 
chemical potential for a quark species of color i and flavor / thus reads 

\i{ = [i — q S lie + t\ \i3 + t\ us , (26) 

where qf is the electric charge of flavor / (q u = 2/3, q d,s = —1/3), fi e is the electron chemical potential, 
tf and are the color charges associated with the third and eighth generator of [SU(3) C ], respectively 
(tf. = 1/2, tg = — 1/2, = 0, = 1/(2 v^3), tf = — 1/a/3), and /i3, /is are the associated color chemical 
potentials. (One could also introduce individual chemical potentials for red, green, and blue quarks, 
but these can be written as linear combinations of //, /13, and //g.) The mismatch in the Fermi surfaces 
of different quark species forming Cooper pairs is then 

2 2 

5k F {? = k F { - k F ] ~-(qt- tf) /i e + (tl - t]) /i 3 + (f? - # ^8 - ^ + ^ , (27) 

where only the first correction in mf//i{ <C 1 was taken into account. Equation (}2Tj) shows that the 
mismatch in the Fermi surface between different quark species is proportional to the electric and color 
chemical potentials, as well as their mass difference. 

The formation of Cooper pairs occurs at the Fermi surface. Typically, a Cooper pair consists of 
fermions with momenta which are in magnitude close to the Fermi momentum, but which have opposite 
directions, such that the total momentum of the Cooper pair is zero. However, when the mismatch 5k F {f 
increases, it becomes increasingly more difficult to form such pairs with zero total momentum. For the 
species with the smaller Fermi surface one may take a fermion right at its Fermi surfaces, but in order to 
match the momentum, one has to go deeper into the Fermi sea of the other species. Pictorally speaking, 
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forming a Cooper pair becomes energetically disfavored once the cost of "diving" into the Fermi sea 
to find such a matching fermion is higher than the gain in condensation energy by forming a Cooper 
pair. Whether this condition is fulfilled depends on the magnitude of the color-superconducting gap at 
the Fermi surface, (j) Q , compared to the mismatch in the Fermi surfaces, 6kp{j . As long as O ^ 
the Cooper-paired state remains the ground state of the system. However, when Skpff becomes of the 
order of O , or even considerably exceeds it, the Cooper-paired state becomes energetically disfavored 
as compared to normal-conducting state |54j . 

It was recently realized, however, that instead of a transition to the normal-conducting state many 
other possibilities can be envisioned. For instance, imagine being in the CFL state and for the moment 
neglect /i e , /i3, and /is in Eq. (J27j) . Then, the CFL state will become energetically disfavored when 
m 2 s /2n exceeds O jSHI- Nevertheless, quark matter will not simply become normal-conducting, because 
there is nothing to prevent the up and down quarks to form a 2SC state. Of course, one cannot simply 
discard fi e , /i^, and /ig from the consideration. Taking these chemical potentials into account to ensure 
overall neutrality with respect to color and electric charges, the 2SC state may become unstable with 
respect to the formation of a gapless superconductor jSS], a crystalline color superconductor jSB], or 
some other state with an even more exotic pairing scenario |oT| H)5] . 

However, also a more conventional pairing scenario is conceivable [52]: the dominant terms in the 
mismatch (J27)) are the ones ~ fi e and the mass difference (the color chemical potentials fi% and /i§ are 
parametrically of order 0o//U <C 0o)- Consequently, instead of realizing one of the more exotic pairing 
scenarios, it could be energetically favorable to simply pair quarks with the same charge and the same 
mass, i.e., of the same flavor. As discussed above, these Cooper pairs must have spin one. Although 
spin-one gaps are orders of magnitude smaller than spin-zero gaps [5TH I5T} I5T?] , the gain in condensation 
energy IS.E co ^ is parametrically larger than for some of the aforementioned exotic pairing scenarios, for 
instance AE con a. ~ ^ 2 4>1 f° r spin-one pairing vs. AE coni i. ~ /uVcK^loff/^o) 4 for the crystalline color 
superconductor j^. (Here, 0loff is the value of the gap in the LOFF phase, while 0o is the gap in 
a superconductor with equal Fermi surfaces for the particle species forming Cooper pairs.) Whether, 
and if yes, which of these pairing scenarios are realized in nature, can only be decided by a quantitative 
comparison of the pressure in the various cases. This has not been done so far and to draw definite 
conclusions about the structure of the phase diagram of nuclear matter at small temperatures and 
chemical potentials of the order of \x ~ 500 MeV appears to be premature at this point. 

On the other hand, it is far simpler to decide what happens to a particular color-superconducting 
phase when one increases the temperature at a given chemical potential. Like in any other super- 
conducting system, thermal motion will break up quark Cooper pairs. In BCS theory, the transition 
between superconducting and normal-conducting phases is usually of second order and occurs at a 
temperature T C BCS proportional to the size of the superconducting gap parameter O) 

T C BCS = -0o - 0.56700, (28) 

7T 

where 7 ~ 0.577 is the Euler-Mascheroni constant. At least in the mean-field approximation, in all color- 
superconducting phases studied so far, T c either rigorously obeys this relation or differs only by a factor 
of order one from it, for details see Sec. 15.31 and Ref. [51 J . The value of the color-superconducting gap 
parameter O is therefore of great importance in order to locate the transition line between the normal 
and the color-superconducting quark matter phases in the nuclear matter phase diagram. In Sec. 03 it 
will be discussed how to compute this gap parameter. Here it suffices to know that an extrapolation 
of the result of solving a gap equation in weak coupling QCD down to moderate densities suggests gap 
parameters of the order of 10 MeV for pairing in the spin-zero channel. NJL model calculations suggest 
somewhat larger values around 100 MeV. With Eq. (|28|) . this would lead to transition temperatures 
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of the order of 6 to 60 MeV. In the spin-one channel, the gaps and critical temperatures are typically 
smaller by two to three orders of magnitude |oTH I51| loTi]. 



3 Lattice QCD 



3.1 Basic Concepts 

The most fundamental approach to compute thermodynamic properties of strongly interacting matter 
and, in particular, its equation of state, are lattice QCD calculations [BU]- In these calculations, one 
directly computes the grand partition function (JTJ) on a discretized space-time lattice, V x 1/T — 
(a a N a ) 3 a T N T , where a a = L/N a is the lattice spacing in spatial direction, a T = 1/ (N T T) is the lattice 
spacing in Euclidean time (i.e. temperature) direction, and N a and N T are the number of lattice points 
in spatial and temporal direction, respectively. Any space-time point in V x 1/T is then parametrized 
as x^ = (r, x) = (a T I, k), with < I < N T , and < i,j, k < N a . Commonly, one uses 

symmetric lattices, where a a = a T = a. A space-time point on the lattice, a lattice site, is then uniquely 
determined by the 4- vector n^ = (Z,n), n = (i,j,k). Quantities with the dimension of energy are 
measured in units of the inverse lattice spacing a -1 , and different lattices are simply characterized by 
their extension N% x N T . The smallest length scale on a lattice is the lattice spacing a, corresponding 
to a maximum momentum scale Auv ~ a -1 - This scale serves as ultraviolet cut-off which regulates the 
ultraviolet divergences commonly appearing in quantum field theories. The largest length scale on the 
lattice is the lattice extension aN a . It determines the minimum momentum scale Ajr ~ (aN^)^ 1 . 

The next step is to define the QCD action S = j x C, with C given by Eq. (J2J), on the discretized 
space-time lattice. As a first guess, one could replace all derivatives with finite differences between 
lattice points. For reasons explained below, this naive prescription is, however, not particularly suitable, 
neither for the gauge field (gluon) nor the matter (quark) part of the action. To find an alternative, 
note that the only condition a discretized version of the QCD action has to fulfill is to reproduce the 
continuum action in the limit a — >• 0. The choice of a discretized QCD action is therefore not unique. 
This apparent shortcoming can, however, be turned into an advantage by choosing a form of the action 
which reduces or even completely eliminates discretization errors (so-called improved or perfect actions, 
respectively). 

To find a suitable discretized version of the gauge field part of the action, one first observes 61j that, 
on a finite-size lattice, the gauge fixing term C gauge in Eq. (j2J) is no longer necessary, as the integration 
over gauge fields becomes convergent. Nevertheless, a naive discretization of the gauge field part of the 
action is still not gauge- invariant, and will remain so even when taking the continuum limit a —>■ 0. It 
is therefore advantageous to formulate the gauge field part of the action in a gauge-invariant form. A 
suitable choice was proposed by Wilson [HTj . 



The link variable describes the parallel transport of the gauge field between two neighbouring lattice 
sites n x and n x + /a a , where fi x = 5^ is the 4-dimensional lattice unit vector pointing in /z-direction. 




(29) 



The sum over n runs over all lattice sites n M and the link variable U n>ll is defined as 




= an' 



(30) 
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Visualizing the product of the four link variables on the right-hand side of Eq. (|29|). one realizes 
that this product transports the gauge field around an elementary lattice plaquette; it is therefore also 
called the plaquette operator. The trace of the plaquette operator, and thus also the Wilson action (|2*9]). 
is gauge-invariant. The Polyakov loop (J7J) is related to the link variables via 

1 N T 

L(x) = — Trn^n.o, n" = (Z,n), xEan. (31) 
iVc 1=1 

Expanding the Wilson action for small lattice spacing a, one obtains the continuum limit 

-P S A^j x {~\ W) + 0(a 2 ) , (32) 

where (3 = 2N c /g 2 = 6/g 2 is the so-called bare coupling. The correction terms to the continuum result 
are of order 0(a 2 ). The construction principle behind an improved action is to add further terms 
to Sa in Eq. (|2Ti|) in order to eliminate corrections of order 0(a 2 ). The improved action then reproduces 
the continuum limit up to corrections of order 0(a 4 ). Repeating this procedure, one can systematically 
eliminate discretization errors up to a given power of a 2 . Extending this procedure in order to eliminate 
all corrections leads to so-called perfect actions 

The naive discretization of the fermionic part of the QCD action is not particularly suitable because 
of the so-called doubling problem for massless fermions on the lattice [HQ]. Fermion doubler states 
originate from the periodicity of the fermion dispersion relation within the Brillouin zone. One obtains 
one extra doubler state per space-time dimension, such that there are in total 2 4 = 16 instead of a single 
fermion species. One way out is to break chiral symmetry explicitly by introducing a mass term. This 
leads to the so-called Wilson fermion prescription [HT]. Wilson fermions eliminate the doubler states 
completely, but they have the disadvantage that one can in principle no longer study the restoration of 
chiral symmetry at the QCD transition. Another possibility is to distribute components of the fermion 
Dirac spinor over several lattice sites. These so-called staggered or Kogut-Susskind fermions [HI] do not 
completely solve the fermion-doubling problem: the number of doubler states is merely reduced to four. 
However, the solution to this problem is to interpret the doubler fermions as different flavor states. 
Hence, the standard staggered fermion action is interpreted as describing QCD with Nf = 4 flavors. 
The advantage of the staggered fermion prescription is that it preserves a subgroup of the original chiral 
symmetry. The chiral condensate is thus an order parameter for chiral symmetry restoration at the 
QCD transition. Other attempts have been made to solve the fermion-doubling problem, while at the 
same time improving (or even preserving) the chiral symmetry of the lattice action. To name a few, 
there are the so-called overlap jHH] > domain- wall [HH] , fixed-point jHZ] or chirally improved [HE] fermions. 

Let us take a closer look into the staggered fermion prescription, where the fermionic part of the 
QCD action reads 

S£ s = J2^nM«^ m , (33) 

n,m 

with the inverse staggered fermion propagator 

+ 2 ( S n+6,m U n,o ^ - S nm+6 U ] mfl e - *) + 5 n , m m . (34) 

Here, the fermion mass (in units of the inverse lattice spacing) is denoted as fh = am. This notation 
prevents confusion of the fermion mass with the lattice site vector m = m' i . The chemical potential 
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(in units of the inverse lattice spacing) is jl = a/x. As shown in Ref. [HO], the correct prescription to 
introduce the chemical potential in the discretized fermion action is as indicated in Eq. ([34)1. i.e., in 
exponential form on a temporal link. 

The fermionic part of the QCD action is bilinear in the Grassmann fields if) and if). The fermion 
fields can thus be integrated out exactly. The result for the QCD partition function (Q) on a discretized, 
4-dimensional N% x iV T lattice is 

Z{N a ,N T ,fi,fh,p) = I HdU n>li [debM KS (m,fi,U)] Nf/A e~ pSA . (35) 

The integration is over all link variables U n ^. The power Nf/4 of the fermion determinant takes 
into account that, in the continuum limit, the standard staggered fermion prescription leads to Nf = 4 
fermion species. In order to obtain results with less than Nf = 4 flavors, one has to take the appropriate 
root in Eq. (|35|). In this way, one can also obtain results for fermions with different masses. For instance, 
in order to compute the partition function for two light (say, up and down) and one heavy (say, strange) 
quark flavor (also called the "2 + 1" flavor scenario) one replaces the fermion determinant in Eq. (JHHj) 
by the product [detM KS (m q , fi q , 

For vanishing quark chemical potential, p, = 0, the fermion determinant in Eq. ()35|) is real and 
positive, and standard Monte Carlo methods [00] can be applied to evaluate the integral over the 
link variables U n ^. However, for nonzero quark chemical potential, the fermion determinant becomes 
complex. It is clear that the partition function itself cannot have an imaginary part, thus the imaginary 
part of the fermion determinant has to cancel when integrating over U n ^ m . However, for a particular 
configuration of the gauge field, or equivalently, the link variables U n ^ on the space-time lattice, the 
real part of the fermion determinant is no longer strictly positive. This so-called sign problem prevents 
the application of standard Monte Carlo techniques to evaluate the partition function. For this reason, 
most lattice QCD calculations have been performed at zero quark chemical potential, with data reaching 
an impressive level of quality. Results for the QCD phase transition and the equation of state, i.e., the 
pressure as a function of temperature, are presented in Sees. 13.21 and 13.31 respectively. Only recently, 
attempts have been made to compute the partition function also for nonzero values of the quark chemical 
potential. This will be discussed in Sec. 13.61 

Finally, let us note that, in order to extract continuum physics from lattice calculations, one has to 
extrapolate the results to the case of vanishing lattice spacing, a — > 0. In order to change the value 
of a, one has to change the value of bare coupling f3 = 6/g 2 - Since QCD is an asymptotic theory, the 
strong coupling constant at the momentum scale a" 1 vanishes as a goes to zero, g(a) — > for a — > 0. 
This, in turn, implies that (3(a) — > oo as a — > 0. Asymptotically, the relation between a and (3 is given 
by the leading-order renormalization group result 

(6b y bl/i2b2 ° ] / M 11-2^/3 102 - 38^/3 

where A^ is the lattice scale parameter that can be unambiguously related to the scale parameter in 
other regularization schemes, for instance A^ in the MS scheme. 

In principle, this allows one to convert the value of a physical quantity, say the pressure, which 
on the lattice is computed in units of a -4 , i.e., as p = a 4 p, into physical units, i.e., MeV 4 . However, 
as Eq. ()36|) is strictly valid only for asymptotically small values of a, in practice one uses a different 
prescription. Consider a physical quantity, for instance a hadronic mass m#, which is well-known in 
the continuum. Compute this quantity on the lattice, where its value is given in units of the inverse 
lattice spacing, fhn = arriH- Then, any other quantity with the dimension of energy can be determined 
in units of m#, say the temperature, which is T/rriH = (aN T )~ l j '{jhu I '0) = l/(m>HN T ). 
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When decreasing the lattice spacing a (by increasing the value of the bare coupling f3), the temper- 
ature T = l/(aN T ) increases, if one keeps the number of lattice points in temporal direction N T fixed. 
(Simultaneously, for a fixed number of lattice points in spatial direction N a , the volume V = (aN a ) 3 
decreases.) Therefore, in order to determine the temperature dependence of a quantity, one simply has 
to compute it on a lattice with a fixed number of temporal points N T , but for different values of a, 
resp. f3. Thus, one often finds lattice data presented as a function of (3 rather than as a function of T. 
Both presentations are equivalent, but note that the temperature increases with the bare coupling f3. 
One should therefore never confuse the bare coupling f3 with the quantity (3 = (fcgT) -1 from thermody- 
namics and statistical mechanics, which decreases with temperature. (The way (3 appears in Eq. (jHSJ) 
certainly does not help to avoid this mistake.) 

When extrapolating lattice results to the continuum limit a — > 0, one does not simultaneously 
want to increase the temperature T or decrease the volume V of the system. Rather, one has to 
ensure that these quantities are kept fixed. In other words, the continuum limit a — > is obtained 
by simultaneously increasing the number of lattice points in space and time direction, N a , N T — > oo, 
such that aN a = V 1/3 and aN T = l/T are constant. This is obviously quite costly numerically. There 
is, however, also another problem of numerical nature with this limit. Consider, for instance, a lattice 
computation of the pressure, which yields values for the quantity p = a 4 p = (p/T 4 )N~ 4 . A given value 
for the physical temperature corresponds to some value for the physical pressure, such that p/T 4 assumes 
a certain value. Consequently, as N T — > oo the numerical value for p on the lattice rapidly decreases as 
N~ 4 when a —>■ 0. Since lattice QCD calculations are subject to statistical errors, it therefore becomes 
increasingly more difficult to extract the physically relevant quantity from the statistical noise. It is 
thus important to use improved actions (see discussion above), which reduce the discretization errors 
and allow one to perform calculations for moderate values of iV T where p is still significantly larger than 
the statistical noise. 

Finally, not only is one interested in the continuum limit for a finite volume V, but one would like 
to extrapolate to the thermodynamic limit V = (aN a ) 3 — > oo as well. At a given nonzero temperature, 
however, l/T = aN T remains finite (in fact, it decreases as T increases). Therefore, simulations at 
nonzero temperature, which aim towards the thermodynamic limit, require N a ^> iV T , which represents 
another numerically expensive condition. Nowadays, typical "hot" lattices have space-like extensions 
N„ ~ 16 — 32 while the time-like extension is N T ~ 4 — 8. The only situation where one also has to 
have a large extension of the lattice in the time direction is the zero-temperature case, T = 0. "Cold" 
lattices typically have N a = N T ~ 16 — 32. 

How close to the thermodynamic limit are present-day lattice QCD calculations? Suppose one is do- 
ing a simulation at the critical bare coupling, i.e., where the QCD transition occurs (see Sec. l3.2|) . Phys- 
ical values for the transition temperature are of the order of T c ~ 150 MeV, cf. Sec. 12.3.31 Consequently, 
the lattice spacing at the critical bare coupling j3 c is a c = 1/(T C N T ). For temporal lattice extensions of 
the order of N r ~ 4 — 8, this corresponds to values a c ~ 0.15 — 0.3 fm. For typical spatial lattice exten- 
sions N a ~ 16-32 on a "hot" lattice, the physical volume is then V ~ (2.5-10) 3 fm 3 ~ (15-1000) fm 3 . 
Is such a system sufficiently close to the thermodynamic limit? The answer is not necessarily "no", 
as this depends on how large the system is in comparison to the size of its constituents. The latter 
can be estimated via their Compton wavelength Ac = mT 1 . For nucleons, the Compton wavelength 
is Ac ~ 0.2 fm, so that many nucleons would comfortably fit into the system. (Of course, this is an 
optimistic estimate: taking the nuclear charge radius r ~ 0.8 fm instead of the Compton wavelength 
drastically worsens the situation.) For a pion, Ac ~ 1-4 fm, such that the lattice volume for these light 
particles appears to be on the verge of being too small (unless the pion becomes much heavier at the 
phase transition, cf. Sec. l4.6|) . In any case, not more than a few pions would fit into the physical volume, 
which certainly casts doubts on whether one is able to reach the thermodynamic limit with present-day 
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Figure 3: Deconfinement and chiral symmetry restoration in QCD with Nf = 2 dynamical 
quark flavors, (a) The expectation value of the Polyakov loop and the corresponding suscep- 
tibility as functions of the bare coupling, (b) The chiral condensate and the corresponding 
susceptibility as functions of the bare coupling. From Ref. |To] . 



lattice sizes. 

In the following, results from lattice calculations at zero and nonzero quark chemical potential will 
be reviewed. For /i = 0, a wealth of data is available; for the purpose of this introductory review, I only 
focus on the QCD phase transition, the equation of state, the heavy quark free energy, and mesonic 
spectral functions. The case /i ^ has only recently received a fair amount of attention. The main 
activity is still to find solutions of (or ways around) the sign problem of the fermion determinant. For 
more details, see Refs. [T5t fTB] . 



3.2 The QCD phase transition 

As already discussed in Sec. 12.31 lattice QCD calculations have numerically established the existence 
of the quark-hadron transition. Figure El (a) shows the expectation value of the Polyakov loop, (L(x)), 
with L(x) as defined in Eq. (fBTj). as a function of the bare coupling (3 (i.e., as explained in Sec. 13.11 
as a function of temperature) for Nf = 2 quark flavors. For the pure gauge theory, i.e., for quark 
masses m —>■ oo, the Polyakov loop is an order parameter for the deconfinement transition: it changes 
its value from zero in the confined phase below T c to a nonzero value in the deconfined phase above 
T c , cf. discussion in Sec. 12.3.21 However, the presence of dynamical quarks in the calculation of Fig. |3] 
breaks the Z(3) symmetry of the pure gauge theory explicitly Thus, the transition is no longer of first 
order, but crossover. This is also observed in the data. 

In Fig. 01 (b) the chiral condensate {tpip) is shown as a function of the bare coupling (3. For vanishing 
quark masses, the chiral condensate serves as an order parameter for chiral symmetry breaking: it is 
nonzero below and vanishes above T c , cf. discussion in Sec. 12.3.31 Since the calculations of Fig. |3]have 
been done for a nonzero quark mass, chiral symmetry is explicitly broken. Consequently, the chiral 
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Table 2: The critical temperature T c for QCD with different quark flavors (extrapolated to 
the chiral limit) and the pure [577 (3) c ] gauge theory (Nf = 0). For the "2+1" case, T c is 
close to the 2-flavor case. From Refs. [TBI I16j. 
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1 ( 1 ± 4 ) Me V 


Wilson fermions 
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(173 ± 8) MeV 


Kogut-Susskind fermions 
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(154 ± 8) MeV 


Kogut-Susskind fermions 





(271 ± 2) MeV 


pure gauge theory 



transition is not of second order, as expected for Nf = 2 flavors, but crossover, which is also seen in 
Fig. 0(b). 

Also shown in Figs. H3 (a) and (b) are the susceptibilities corresponding to the Polyakov loop and 
the chiral condensate. These are defined as 

X L = N* (<L 2 > - (L) 2 ) , X m = £r $ 4>) • (37) 

These quantities have a maximum at the value of f3 where the Polyakov loop and the chiral condensate 
change most rapidly. This value is the critical bare coupling /3 C , which corresponds to the critical 
temperature T c for the QCD transition. In this way, one can define a critical temperature, even if the 
transition is not of first or second order, but only crossover. The interesting observation one can make 
from Fig.[3]is that f3 c assumes the same value for the deconfinement transition as for the chiral symmetry 
restoration transition. A possible explanation for this strong correlation between deconfinement and 
chiral transition is provided by the Polyakov-loop model of Ref. [S3], see also Sec. 14.51 Current results 
for the phase transition temperature in the pure [S77(3)] c gauge theory, as well as in QCD with different 
flavors, extrapolated to vanishing quark mass (for details, see Ref. [71]) are summarized in Tabled 

Lattice QCD calculations have also begun to explore the quark-mass diagram discussed in Sec. 12.3.41 
in order to decide the question about the order of the QCD transition. The present knowledge is 
summarized in Fig. 0] The open triangles are results from Ref. [72] and correspond to the line of 
second-order transitions separating the first-order from the crossover region in Fig. [21 The other data 
points confirm that the transition is of first order below the second-order line and crossover above. It is 
somewhat difficult to locate the physical point on this diagram. Naively, one would think that it suffices 
to determine the lattice spacing a in physical units, after which one finds the physical point in lattice 
units via multiplying the physical quark masses by this value of a, fn^J 5 ' = am^f s \ This deceptively 
simple method does not work in practice, because the physical mass in lattice units rh^ h / s - also receives 
contributions from renormalization, which violate this simple relationship. Present estimates seem to 
indicate, however, that the physical point is deep in the crossover region [76 . 
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Figure 4: The quark-mass diagram as computed in lattice QCD. Below the line of open 
triangles, the transition is of first order, above it is crossover. The data points labelled 
"Columbia" are from Ref. jZHj, the ones labelled "JLQCD" are from Ref. [71] . The other 
points are from Refs. j?3 ES] . 



3.3 Equation of state 

The equation of state is determined by the pressure p(T, u) as a function of temperature T and chemical 
potential u. Hence, according to Eq. (£Q) one has to compute (T/V) InZ. From the pressure, other 
thermodynamic quantities can be derived via differentiation, cf. Eq. (JHJ), and the fundamental relation 
of thermodynamics, e = Ts + un — p. For any quantum field theory in the continuum as well as on the 
lattice, the calculation of the absolute value of (T/V)\nZ is plagued by ultraviolet divergences [TO] . 
These arise from vacuum fluctuations and have to be subtracted in order to obtain a finite value for 
p(T,fj). The simplest way to achieve this is to subtract the value of (T/V)lnZ in the vacuum, i.e., at 
T = n = 0, 

p (T,u) = ^lnZ-(^\nz) . (38) 

V \ V J T=/i=0 

In this way, the value of the pressure in the vacuum is normalized to zero, p(0, 0) = 0. 

A direct computation of the pressure using this formula is still cumbersome, because it requires the 
calculation of the absolute values of \iaZ(T, V, //) and In Z(0, V,0). which then have to be subtracted 
from each other. On the lattice, it is much simpler to compute average values of quantities. Therefore, 
one uses the following method to determine the pressure. First, note that 
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-InZ - - InZ 

V \V 



T=n=0 



_ 4 /In Z(N a ,N T , (3, fh,p) In Z(N a , N a , ft,rh, £i)\ 

= K { m J • (39) 



The assumption underlying this identity is that one can approximate the vacuum subtraction by 
the value of (T/V)lnZ computed on a "cold" lattice with N T = N a , but at the same value of the 
bare coupling (5 (i.e., with the same lattice spacing a) as for the "hot" lattice (where N T N a ), 
(T/V) lnZ|y =At=0 = (aiV CT )~ 4 In Z(N a , N a , /?, rh, jj). Now introduce the expectation value of the (dimen- 
sionless) Wilson action density (§a) = a 4 y(SA) = (N^N t )^ 1 (Sa), 

(3a) = -— 7£ ln Z{N„N T ,p,m,jk) , (40) 
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Figure 5: Left panel: Pressure (divided by T 4 ) as a function of temperature for the pure 
gauge theory and for QCD with two and three light flavors, as well as with two light and a 
four times heavier quark flavor (curve labelled "2 + 1"). Arrows denote the corresponding 
Stefan-Boltzmann values. For the "2+1" case, the arrow is slightly below the three-flavor 
case, due to the nonzero strange quark mass. Right panel: Pressure normalized to its 
corresponding Stefan-Boltzmann value as a function of temperature normalized to the corre- 
sponding transition temperature T c for the four cases shown in the left panel. From Ref. [To] . 



and its zero-temperature value (§a)o, which is computed on a "cold" lattice, i.e., by setting N T = N a 
on the right-hand side of Eq. ()4(jp. The quantity p/T A in Eq. (JH9*|) can now be obtained through an 
integration of Eq. (pfUj) with respect to the bare coupling (3, 

£ - f* = ~K £ d(3' ((s A ) - ( S ~4> ) . (41) 

In order to determine the second term on the left-hand side, pi/T 4 , one would like to choose a rather 
low value for the temperature 7\. For temperatures T\ <C tuh, where m# is the lightest hadronic 
particle, the pressure is exponentially small, pi/Tf ~ exp(— m#/Xi). This argument holds to very good 
approximation in the pure gauge theory, since the lightest glueball state has a mass of order 1 GeV. 
It does not hold in full QCD in the chiral limit, where there are Nj — 1 massless Goldstone particles, 
cf. discussion in Sec. 12.3.31 For lattice QCD calculations, however, chiral symmetry is always broken 
by a nonvanishing dynamical fermion mass, thus m# is always positive. To very good approximation 
one may therefore set pi/Tf = 0. Once the pressure is known, other thermodynamic quantities can be 
determined from thermodynamic identities, for more details see Ref. 1 1 ll . 

As discussed in Sec. 13.11 before one can draw definite conclusions about the thermodynamic proper- 
ties of hot quark-gluon matter, one has to extrapolate the lattice data to the continuum limit (and hope 
that present-day lattices are sufficiently large to be close to the thermodynamic limit). This has been 
done in Fig. which shows the pressure (normalized to T 4 ) as a function of temperature (in physical 
units) for the case of two light flavors, three light flavors, and the "2 + 1" case, i.e., two light plus one 
heavy flavor [77], in comparison to the pressure for the pure [S77(3) c ] gauge theory [78| . 

One observes that the pressure is rather small at low temperatures. This is to be expected, as the 
contributions from hadronic resonances (or, in the pure gauge theory, from glueballs) to the pressure 
are exponentially suppressed, Ph/T 4 ~ exp(— tuh/T) for a hadron (glueball) of mass m#. However, 
at the critical temperature T c for the QCD transition (cf. Table EJ), the pressure increases rapidly, 
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approaching the so-called Stefan-Boltzmann limit psb for a system of quarks and gluons as T — > oo. 
The Stefan-Boltzmann limit is the pressure for an ideal (i.e., non- interacting) ultrarelativistic gas of 
particles. For an ultrarelativistic gas at // = 0, the temperature is the only scale with the dimension 
of energy, consequently Psb/T 4 = const.. The value of this so-called Stefan-Boltzmann constant only 
depends on the number of degrees of freedom in the system. For an [SU(N C )] C gauge theory with Nf 
massless quark flavors one obtains 



Psb 

y4 



2(N 2 c -l) + 2N c N f 7 - 



7T 2 



(42) 

90 V ; 



Here, the first term in brackets is the contribution from the gauge fields, while the second corresponds 
to that from the matter fields. The factors of 2 in these terms arise from the spin degrees of freedom of 
massless gauge fields and fermions. The factor iV 2 — 1 counts the number of gauge fields which are in the 
adjoint representation of the gauge group. The factor N c Nf counts the number of colors and flavors of 
the fermions which are in the fundamental representation of the gauge group. The factor 7/4 accounts 
for the difference between Bose-Einstein and Fermi-Dirac statistics and for the fact that at fi = there 
are as many antifermions as fermions. Finally, the factor 7r 2 /90 is the value of the (dimensionless) 
Bose-Einstein integral (1/67T 2 ) J °°da;x 3 (e x — 1) _1 occurring in the calculation of the pressure. 

In all cases, the function p/T 4 approaches the Stefan-Boltzmann value for a gas of quarks and gluons, 
which indicates that there is indeed a transition from hadronic degrees of freedom to quark and gluon 
degrees of freedom, i.e., from hadronic matter to the QGP. From the behavior of the Polyakov loop 
and the chiral condensate discussed in Sec. 13.21 in the QGP color charges become deconfined and chiral 
symmetry is restored. However, the approach of the pressure to the corresponding Stefan-Boltzmann 
value is rather slow; even at temperatures ~ 3T C , deviations are typically of the order of 20%. This 
indicates that at such temperatures the QGP cannot really be considered as a non-interacting gas of 
massless quarks and gluons. 

In order to understand the deviations from the Stefan-Boltzmann values, one has to resort to an- 
alytic calculations of the pressure, taking into account interactions between quarks and gluons. In an 
analytic approach, deviations from the ideal-gas behavior are well under control and can be physi- 
cally interpreted. For instance, in a perturbative calculation of the QCD pressure, deviations from the 
Stefan-Boltzmann limit are due to corrections proportional to powers of the strong coupling constant, 
see Sec. 14.11 Another possible explanation for the deviation of the pressure from Psb/T 4 is that quarks 
and gluons are actually quasiparticles, i.e., they are not massless, but assume a thermal mass due to 
interactions with the hot environment, see Sees. 14.31 and 14.41 

An important step to understand the deviations from ideal-gas behavior might be the observation 
that, when normalizing the pressure to the corresponding Stefan-Boltzmann value and the temperature 
to the critical temperature, the curves p/psb as a function of T/T c exhibit a universal behavior for the 
pure gauge theory and for QCD with various dynamical quark flavors, see right panel of Fig. 03 A possible 
explanation for this behavior is provided by the Polyakov-loop model of Ref. [7Dj, see Sec. 14.51 where 
the dynamics of chiral symmetry restoration is exclusively driven by the dynamics of the deconfinement 
transition. 



3.4 Heavy quark free energy 

The behavior of the heavy quark free energy as a function of temperature is another indication for 
deconfinement in the QGP. The heavy quark free energy Fqq(R, T) is the free energy of a heavy quark 
and an antiquark, separated by a spatial distance R, at a temperature T j21]. It is related to the 
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Figure 6: Left panel: The heavy quark free energy (in units of temperature) as a function 
of the distance (in units of inverse temperature) for various values of the bare coupling (3 in 
the confined phase. Right panel: The same for various values of (3 in the deconfined phase. 
In this phase, the Polyakov-loop correlation function of Eq. ()43j) is normalized to (|L(0)| 2 ). 
Calculations are done for pure [S77(3) c ] gauge theory on a 32 3 x 4 lattice. The critical bare 
coupling for this lattice is (3 C = 4.0729(3). From Ref. [75j . 



Polyakov-loop correlation function via 



cxp 



T 



(L(O)Lt(x)), i2=|x|, (43) 



where L{x) is the Polyakov-loop operator defined in Eq. (|3ip. At T = 0, the heavy quark free energy is 
identical to the heavy quark potential, Vqq, which is expected to have a form motivated by the string 
model, 

V QQ( R ) = + o- R + const - ■ ( 44 ) 

The second term ensures confinement of color charge due to the linear increase of Vqq with distance. 
The constant a is the string tension. The first term is an attractive Coulomb-like contribution arising 
from fluctuations of the string. 

Lattice QCD data [79 confirm these expectations, cf. Fig. El Below T c (left panel of Fig. EJ, at 
small distances the heavy quark free energy is dominated by an attractive Coulomb-like part, while 
at large distances it linearly rises with the distance, indicating confinement. The linear rise becomes 
less pronounced with increasing bare coupling (3, indicating that the string tension decreases with 
temperature. 

At temperatures above T c , color charges are deconfined, i.e., the linearly rising part of the potential 
in Eq. (J44j) has to vanish, leaving only a Coulomb-like part. The latter is, however, screened due to the 
presence of a hot medium. This is confirmed by lattice QCD data above T c (right panel of Fig- EI) • It 
turns out [7^] that a fit to the numerically computed potential can be achieved by the formula 

F QQ( R > T ) = < T ) -^{T)R Ur>) 

T " (RT) d ( T ) ' 1 J 

where fi(T) is the temperature-dependent screening mass (or inverse screening length). This function 
is shown in Fig. 
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Figure 7: The temperature dependence of the screening mass fi(T) (in units of temperature) 
as obtained from a fit of Eq. to the lattice data of Fig. The solid line is a fit inspired 
by the perturbative expectation fi pcT t.(T) = 2m D (T). From Ref. [79] . 



While the qualitative picture of deconfinement and screening of color charges is certainly applicable, 
the deconfined gluon-plasma phase cannot be described perturbatively at temperatures in the range 
from T c to a few times T c . This is indicated by the fact that the fit function d(T) in Eq. ()45|) is 
temperature-dependent [79] and always smaller than ~ 1.5 in the range of temperatures considered 
here, while from perturbation theory one expects <i per t. = 2. Furthermore, the screening mass n(T) 
deviates from the perturbative value n per t.(T) = 2m D (T), where m D {T) = gT is the Debye mass in 
a hot gluonic medium. The solid line in Fig. [7J represents a fit inspired by /i per t. (T) to the two data 
points corresponding to the highest temperatures, for details see Ref. [79J. One observes that while 
the qualitative behavior of the data follows the perturbative expectation at large temperature, near 
T c the data strongly deviate from the perturbative result. They even suggest that the screening mass 
goes to zero when T — > T c . This is an indication for critical behavior and is naturally explained by the 
Polyakov-loop model of Ref. [70], see Sec. 14.51 

More recent developments in the study of the heavy quark free energy include a calculation in full 
QCD with dynamical quark flavors |71j . Below T c , the string breaks when creation of dynamical quark- 
antiquark pairs becomes energetically favorable. Consequently, the heavy quark free energy saturates 
at larger distances instead of increasing linearly. In another recent paper [SU] the color-singlet and 
color-octet contributions to the heavy quark free energy were studied separately within the pure gauge 
theory, using a novel prescription to renormalize the expectation value of the Polyakov loop. It was 
found that the singlet and octet contributions only deviate at smaller distances. As expected, the 
color-octet channel is repulsive, while the color-singlet channel is attractive. 

3.5 Mesonic spectral functions 

The correlation function of a mesonic state Gh(t, x) is defined as 

G h (t, x) = $(0) T H V>(0) $(r, x) T H ij(r, x)> , (46) 



32 



0.6 



_ «oooo^x xxxxxx 


Jt i — I — i 

a ^ 

Xv 

x:; < >< ><X>O0<X><><X><X>0<X><><XKK 




m = 0.02 " 
L = 8 



5.24 5.26 5.28 5.3 5.32 5.34 



Figure 8: Mesonic masses as a function of the bare coupling, computed on an 8 3 x 4 lattice 
for QCD with Nf = 2 dynamical flavors. The vertical line indicates the critical bare coupling 
f3 c for the phase transition. From Ref. [To] . 



i.e., it is the overlap of a mesonic state with quantum numbers determined by the 4x4 Dirac matrix 
Th at the origin with the same mesonic state at (r, x). Fourier-transforming Eq. ()46|) with respect 
to the spatial variable, one obtains the mixed correlation function Gh(t, p), which has the spectral 
representation 

roc do; COsh UJ (T — ^p) 

G H (t, p) = / — a H (u, p) X - 2r ^ J . (47) 

J o 2vr S i n h (^J 

Here, <7fl-(u;,p) is the spectral density in the quantum number channel under consideration. Suppose 
the spectral density is dominated by a single, stable, mesonic state with mass m#. In this case, 
<Jh(u, 0) = 7i\ 2 5(uj—m H )/uj, where A 2 is a constant with the dimension [MeV 4 ]. Then, the susceptibility 

V r l / T r , 

XH = - / dr / d 3 xG ff (r,x) (48) 
1 Jo Jv 

assumes the value xh = (V/T)X 2 mjj 2 , i.e., it is proportional to the inverse mass (squared) of the meson. 
In a lattice QCD calculation, one can thus infer the mass of a mesonic state in a given quantum number 
channel from the corresponding susceptibility. 

The masses for the pion, the a meson, and the a® meson computed in this manner are shown in 
Fig.|Hlas a function of the bare coupling (i.e., the temperature). The results indicate restoration of chiral 
symmetry, i.e., the mass of the pseudoscalar meson (pion) becomes degenerate with those of the scalar 
mesons (cr and a®) at large temperature. The fact that the pion and the a meson become degenerate 
in mass at smaller temperatures than the pion and the a>o meson indicates that the SU(2) r x SU(2)i 
symmetry is restored prior to the U(1)a symmetry when increasing the temperature. 

Instead of the susceptibility, one could also try to compute the complete spectral density an(uj,p) 
of a mesonic state from lattice QCD data. An important motivation for such a calculation is the fact 
that the spectral density in the vector channel, ay(uj, p), is directly proportional to the rate of dilepton 
emission [81J, which is an experimentally observable quantity 
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Figure 9: Left panel: The spectral density in the vector channel, divided by u 2 , for various 
temperatures. Right panel: The dilepton rate computed in lattice QCD for two different 
temperatures as a function of energy (in units of temperature). Also shown is the Born 
approximation and the result from HTL-perturbation theory. From Ref. |16j . 



To determine oy(u;, p), one would have to perform an inverse Laplace transformation of Eq. (|47j). This 
requires complete knowledge of the correlation function GV(t, p) on the left-hand side. On the lattice, 
however, this function is only known at a few discrete points in r-direction. Moreover, its value at these 
points is subject to statistical fluctuations. Consequently, a computation of ov(u;, p) via inversion of 
Eq. (jT7|) with lattice data for Gy(r, p) is impossible. 

Nevertheless, a solution of this problem is provided by the so-called "Maximum Entropy Method" 
(MEM). The basic idea is to construct that particular spectral density cr#(u;, p) under the integral in 
Eq. (J47|) . which is the most probable one to yield a given correlation function Gh{t, p) on the left-hand 
side of that equation; for details see Ref. jHH|. Figure M shows the spectral function in the vector meson 
channel computed with this method [S3] (left panel) and the corresponding dilepton rate (|49|) (right 
panel). One observes that the peak in the spectral density broadens and shifts towards larger energies 
as the temperature increases. Consequently, the dilepton rate is depleted for small dilepton energies. 
This behavior is in stark contrast to the dilepton emission rate computed in the Born approximation 
and in the so-called "hard thermal loop" (HTL-) resummation scheme, which are also shown in Fig. El 

Finally, the low-energy behavior of the spectral density determines the value of transport coefficients 
in a hot medium [S3]. I do not elaborate further on this point, as it concerns the non-equilibrium 
properties of the QGP, which are beyond the scope of the present review. 

3.6 Nonzero Chemical Potential 

As discussed in Sec. 13.11 for nonzero values of the quark chemical potential, /i ^ 0, a straightforward 
evaluation of the QCD partition function on the lattice is not possible due to the sign problem of 
the fermion determinant. However, for sufficiently small /x progress has recently been made by apply- 
ing methods which explicitly avoid the sign problem. Most notably among these are multiparameter 
reweighting i32J, Taylor expansion around \x = 85 , and analytic continuation from imaginary values 
of /a, where the fermion determinant is real- valued and positive, to real values of /i jSSj. For the sake 
of brevity, here I only discuss the multiparameter-reweighting method; for a detailed comparison of all 



34 



approaches see the review [TEj . 

The multiparameter- reweighting method is based on the so-called Glasgow method [HI]. The idea 
of the Glasgow method is to treat the fermion determinant at nonzero fx in the partition function (|35j) 
as an observable rather than as a part of the integration measure. The integration measure itself is 
computed with a fermion determinant at fx = 0, which is real-valued and positive and thus causes no 
problems when applying standard Monte Carlo methods to sample gauge field configurations, 



(50) 



•> n,M det M (m, 0, C/) \ det M (m, 0, U) / A=0 

Here, the expectation value (O)^ =0 of an operator O is defined with respect to an ensemble of gauge 
fields and fermions at zero quark chemical potential. 

This method is limited to small values of fx. In order to understand this, one has to remember 
the principle behind a Monte Carlo computation of the functional integral (|5()j) [60J. A Monte Carlo 
computation assumes that in order to obtain a reasonable approximate value of the functional integral in 
Eq. (joUj) it suffices to sum only over (a few 10 2 to 10 4 of) the "most probable" gauge field configurations, 
rather than performing the integrals over the link variables U n ^ explicitly. The "most probable" gauge 
field configurations are obviously those which minimize the action Sa- However, the most probable 
configurations at fx ^ are not the same as the ones at fx = 0. Thus, approximating the partition 
function at fx ^ by configurations obtained for fx = becomes increasingly worse as fx increases. In 
other words, the "overlap" between the ensemble of most probable configurations at fx = and the 
ensemble that consists of the configurations which are actually most probable at fx ^ diminishes. 

It has been recently realized jH2j that a way to increase this overlap is to also include the exponential 
of the action into the operator which is averaged over the ensemble, 

Z(N a ,N T ,f3,m,fl) EE ( _ PqSa j-^ ^ n rr\ ) > ( 51 ) 



detM(m, 0, U) 



A=o,/3 



i.e., the ensemble one averages over is generated at fx = and a value f3$ for the bare coupling. In this 
way, one not only reweights the ensemble in the parameter fx, as in the Glasgow approach (}50|) . but also 
in the bare coupling (3 (hence the name "multiparameter reweighting"). 

How does one choose the second reweighting parameter (31 This depends on which physical question 
one asks. Suppose one wants to compute the QCD phase transition line for nonzero values of fx. One 
first generates an ensemble at p, = and f3$ = (3 C - This ensemble is "maximally" critical in the sense 
that it is generated at the phase transition point (f3, fx) = (/3 C , 0) (which corresponds to the point 
(T,fi) = (T c , 0) in the continuum). For each nonzero value of fx one then determines (3 such that one 
remains on the phase transition line. 

The criterion for "remaining on the phase transition line" is the position of the Lee- Yang zeroes 
/?*, /?!,... of the partition function Z in the complex /5-plane jHHI- For a given set of parameters 
N a , N T ,m, fx, there are many Lee- Yang zeroes, i.e., roots of the equation Z{N a , N T , (3*,fh, fx) = 0. 
(In fact, the total number of Lee- Yang zeroes, M, increases linearly with the volume of the system, 
M ~ N%.) In the case of a first-order phase transition, one root, say /?*, has a vanishing imaginary part, 
i.e., it lies on the positive real (3 axis. Then the value of f3 which corresponds to the phase transition 
line in the (f3, fx) plane is (3 = Re/3*. 

Note that in a finite system, such as the space-time lattice considered in lattice QCD, all Lee- Yang 
zeroes have nonzero imaginary parts. Then one has to study lattices of different sizes and extrapolate 
to the infinite- volume limit. In Ref. [S2], this is done via linear extrapolation in the variable 1/V, 
(3*{V) = /3*(oo) + a/V. In the case of a crossover transition, the imaginary parts of the extrapolated 
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Figure 10: The QCD phase transition as computed in lattice QCD |32j . To the left of the 
critical point the transition is crossover, to the right it is of first order. At the critical point, 
the transition is of second order and in the universality class of the Ising model. Note that 
\i B = 3/i. 



Lee- Yang zeroes never vanish. In this case, the value of (3 corresponding to the phase transition line is 
determined by the real part of the Lee- Yang zero with the smallest imaginary part. 

The phase transition line calculated in this way is shown in Fig. El It agrees with the expectations 
discussed in Sees. ECU 13.21 there is a line of first-order phase transitions, ending at the point (T, fi) cr = 
(160 ± 3.5, 242 ± 12) MeV, at which the transition is of second order. To the left of this point, the 
transition is crossover. One should mention that the lattice QCD calculation underlying Fig. El was 
done on fairly small lattice sizes, with probably unrealistically large quark masses. As discussed in 
Sec. 12.3.41 for smaller quark masses the endpoint should move towards the temperature axis. For three 
massless flavors, it should reach the temperature axis, since in this case the transition is of first order. 
For realistic quark masses, however, as discussed in Sec. 13.21 the transition is crossover at fi = 0, and 
the line of first-order transitions should always end at some nonzero value of \i. 

The position of the phase transition line determined by multiparameter reweighting is in good 
agreement with that computed by the other approaches mentioned previously, namely the Taylor ex- 
pansion method, and method of analytic continuation from imaginary values of [i ^B]- Recent devel- 
opments [HH| are the application of the multiparameter-reweighting method to compute the equation 
of state at nonzero quark chemical potential. Figure El shows the results for the pressure difference 
Ap = p(T, n) — p(T, 0), normalized to T 4 , as a function of T for various values of [i. There is a strong 
increase of Ap around the phase transition temperature. This increase is larger for larger values of \x. 

Multiparameter reweighting, as well as the other aforementioned methods, is restricted to values 
of the quark chemical potential, which are not too large as compared to the temperature. In order 
to compute the partition function of QCD for large quark chemical potential at small or even zero 
temperature, and possibly study the color-superconducting phases of quark-gluon matter, one has to 
resort to other methods. A promising approach is the so-called meron-cluster algorithm which has been 
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Figure 11: Ap as a function of T/T c for various values of \ib- From bottom to top, //# 
100, 210, 220, 310, 530 MeV. From Ref. [BH]. 



shown to solve the sign problem of the fermion determinant for the Hubbard and the Potts model [90J . 
For QCD, as of yet no solution has been found. 

Another possibility is to study a model for QCD which does not have the sign problem. Such a 
model is, for instance, the NJL model which has been investigated on the lattice at T = and /i / 
in Ref. PJT] . Although this model has no color gauge symmetry which could be spontaneously broken, 
and thus strictly speaking cannot exhibit color superconductivity, quarks can still form Cooper pairs 
and the system may become superffuid. Figure ^] shows the chiral condensate, the baryon density, 
and the superffuid diquark condensate as a function of the baryon chemical potential, computed on the 
lattice and then extrapolated to the thermodynamic limit. One observes that, as the baryon density 
increases, the chiral condensate vanishes and the diquark condensate increases, signalling the onset of 
superfluid behavior. These results are in agreement with analytic calculations [HH] for the NJL model 
in the mean-field approximation, which stimulated recent interest in color superconductivity. 



4 Analytic Approaches 



4.1 Perturbation Theory 

The QCD partition function (JI} can be expanded in a power series in the strong coupling constant g. 
In the following, I present the general idea behind this approach, neglecting contributions from gauge 
fixing and from Fadeev-Popov ghosts. Of course, these have to be properly accounted for in order to 
obtain the correct answer; for more details, see Ref. ^U]- I also focus on the case /x = 0. The first 
step is to split the QCD Lagrangian (j2J) into two terms, the non-interacting part, Cq = C g= o, and the 
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Figure 12: The chiral condensate, the baryon density, and the diquark condensate as a 
function of baryon chemical potential /ig. The solid line is the chiral condensate computed 
analytically within the Hartree approximation. From Ref. |91j . 



interaction part, Cj 
S 
Co 



C — Cq. Then, the QCD action S 

S + Sj= f (C + Ci) , 
Jx 



f x L can be written as 
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Here Gq 1 



i'y^dfj, — m is the free inverse quark propagator and Aq ab = (□(? /iI/ — d^d u )8 a b is the 
free inverse gluon propagator (which will eventually receive another contribution from the gauge fixing 
terms neglected here) . The next step is to introduce source terms for fermions and gauge fields (which 
eventually have to be set to zero). One can then replace all fields in Si in terms of functional derivatives 
with respect to the sources, and thus extract e 1 from the functional integral. The functional integration 
over the exponential of the non-interacting part and the source terms is a Gaussian integral and can be 
performed exactly. The result is 
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(55) 



rj=T]=J=0 



where Zq is the partition function for a system of non-interacting quarks and gluons. Obviously, the 
pressure po = (T/V)\nZo is identical to the Stefan-Boltzmann pressure defined through Eq. (j42j) . 
p = p SB . The full pressure in QCD also receives contributions from the remaining two terms in 
Eq. (J55|) . After introducing Feynman rules for propagators and vertices ^U], these terms have a graphical 
representation as an infinite series of diagrams with no external legs. The diagrams can be sorted 
according to powers in the strong coupling constant g associated with the vertices. Thus, one obtains 
a perturbative series in powers of g. 

Inspecting the topology of these diagrams, one would naively conclude that this perturbative series 
is an expansion in powers of g 2 . In fact, it turns out that this is only true at zero temperature [92j. 
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At nonzero temperature, the expansion is in powers of g, rather than g 2 , due to the different infrared 
behavior of a field theory containing massless modes (such as gauge fields) at nonzero temperature. 
Roughly speaking, the difference arises from the infrared behavior of single-particle phase space, which, 
at zero temperature, is ~ dkk 3 , while at nonzero temperature it is ~ Tdkk 2 . The missing power of 
k at nonzero temperature leads to a completely different infrared behavior as compared to the zero- 
temperature case. 

At zero temperature the theory is well-behaved in the infrared and the terms of the perturbative 
series are probably computable to all orders in g 2 j"f"[]. Freedman and McLerran computed the series 
up to terms of order g 4 , for details see Ref. [02] • I do not present a more detailed discussion of their 
results at this point, as quark matter at zero temperature is a color superconductor, cf. Sees. 12.4.31 and 
EJ Color superconductivity is a nonperturbative phenomenon, which cannot be described in a purely 
perturbative calculation of the pressure. 

At nonzero temperature, the infrared behavior of the theory leads to terms proportional to odd 
powers of g in the perturbative expansion (|55|) of the partition function. Technically, they arise from 
a resummation of an infinite subset of diagrams describing the screening of long-range electric fields. 
Moreover, there are infinitely many diagrams at order 0(g 6 ), and the perturbative expansion breaks 
down [TOl - This is sometimes called the Linde problem of QCD, after its discoverer [II]. Nevertheless, 
what is perturbatively computable has been evaluated. These are all terms up to 0(g 5 ), and the terms 
of order 0(g 6 \ng). How to obtain the latter is discussed in greater detail in Sec. 14.21 While the 
terms which are genuinely of order 0(g 6 ) cannot be computed perturbatively, they can in principle be 
evaluated via a lattice calculation. 

At zero chemical potential, the pressure assumes the form 

p = T 4 c + c 2 g 2 + c 3 g 3 + (c' 4 \ng + c 4 ) g 4 + c 5 g 5 + c 6 g e . (56) 

The coefficient Co is equal to the Stefan-Boltzmann constant (|42j) . The coefficient c 2 arises from the 
lowest-order perturbative correction to the pressure of an ideal gas. It consists of two-loop diagrams, 
and was first computed by Shuryak ["""""] 



^-Isr^i*') • (57) 

The computation of the coefficient C3 requires a nonperturbative resummation of plasmon ring diagrams 
in the infrared limit. This was first done correctly by Kapusta [9*4*] . with the result 

N 2 -l / 1 \ 3 / 2 

<* = izHr( N c+n N f) ■ ( 58 ) 
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The coefficient c' 4 has been computed by Toimela ["I*"] . 

< = ^ N '{ N ' + k N ') ' (59) 
The coefficient C4 is due to three-loop diagrams and has been computed by Arnold and Zhai 
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Here, ((x) is Riemann's zeta function, 7^ is the Euler-Mascheroni constant, and p is the renormalization 
scale in the MS scheme, for details see Ref . jSH] • The coefficient C5 arises from corrections to the three- 
loop diagrams due to Debye screening of electric gluons. It has been computed by Zhai and Kastening 
|9~T| and by Braaten and Nieto 
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The coefficient C6 contains terms ~ Ing and constant terms. The former contribute to order 0(g 6 lng) 
and can be evaluated perturbatively, while the latter are genuinely of order 0(g 6 ) and can only be 
computed e.g. via a lattice calculation, for more details, see Sec. 14.21 At finite temperature and nonzero 
chemical potential /i, the contributions of order 0(1), 0(g 2 ), and 0(g A \n.g) to the pressure have been 
computed by Toimela PU] . 

To be more explicit, consider pure [SU(3) C ] gauge theory, i.e., N c = 3, Nf = 0. The pressure up to 
terms of a given order in g is shown in the left panel of Fig. ED The strong coupling constant g is taken 
to be running and evaluated at the scale p,. After applying the principle of fastest apparent convergence 
to minimize the two-loop corrections to the running of g, this scale is chosen as p, ~ 6.742 T; for more 
details sec [100J. The scale p also enters under the logarithms in Eqs. (jHOj) . (jfiT|) . In principle, the 
complete result for the pressure, being a physically observable quantity, must be independent of the 
renormalization scale p. The way this works out is that p under some logarithm, such as occurring 
in Eqs. (|6U|). (|6ip. is cancelled by a similar logarithm from the running of the coupling constant in a 
lower-order contribution. Nevertheless, while terms ~ In p must cancel, there still exist physical terms 
~ Ing, and here g has to be evaluated at the scale p. The cancellation of the /i-dependence holds for 
the complete result for p, but this does not happen if one terminates the perturbative expansion at 
some given order. This is the reason why, for instance, the 0(g 2 ) contribution to the pressure in Fig. [TBI 
is not flat. Here the curvature arises from the logarithmic running of the strong coupling constant with 
the scale p ~ 6.742 T. 

The perturbative series (J56j) converges badly. The second-order term ~ cig 2 gives a negative con- 
tribution to the Stefan-Boltzmann pressure, which is less than 10% of psb at T ~ 10 3 A-^ and at most 
40% of psb at T ~ ^-jts (~ ^c)- However, the next contribution ~ c^g 3 is positive and so large that 
the pressure overshoots psb up to the largest values of T shown in Fig. [T^l The terms of order g 4 are 
again small, but also positive, such that, to order 0(g 4 ), the pressure is larger than p SB . The terms of 
order 0(g 5 ) are negative and so large in magnitude, that the pressure even vanishes at T ~ Aj^. Thus, 
naive perturbation theory is clearly not applicable for temperatures of order T c . 
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Figure 13: Left panel: Pressure (divided by T 4 ) as a function of temperature for the pure 
[S77(3) c ] gauge theory. Several perturbative contributions (up to a given power in g) are 
shown, as well as lattice QCD data from Ref. j7S]. Right panel: Sensitivity of the pressure 
to the value of the constant 5 in the term ~ g 6 . From Ref. 100 . 



In the sections following Sec. 14.2) several ways to improve the situation will be explained. All of 
them are based on the observation that the odd powers of g in the perturbative expansion 1)56)1 are 
responsible for the bad convergence properties, i.e., the latter are caused by the infrared properties of 
QCD. Note that there have also been attempts to improve the convergence properties of perturbation 
theory by using mathematical devices such as Pade approximates 101J and Borel resummation 102J . 
Here, I do not discuss these methods in more detail, because the physical problem of improving the 
description of the infrared sector of QCD cannot be solved in this way. For the sake of completeness, 
one should also mention Ref. 103 , where a phenomeno logical solution to the problem of convergence 
of the perturbative series was presented. 

4.2 Dimensional Reduction 

Consider a quantum field theory at nonzero temperature in the limit T — > oo. In this limit, the 
Euclidean time interval in the partition function (0) shrinks to zero, 1/T — > 0. Consequently, the original 
3+1-dimensional theory reduces to a theory in three spatial dimensions. This is called dimensional 
reduction |104j . What are the degrees of freedom in the dimensionally reduced theory? Recall that 
the compactification of the Euclidean time interval [0, 1/T] at nonzero temperature leads to discrete 
energies for the field modes, the so-called Matsubara frequencies JTQ . Bosonic degrees of freedom 
have periodic boundary conditions in Euclidean time, and thus their Matsubara frequencies are even 
multiples of 7rT, = 2mrT, n = 0, ±1,±2, .... On the other hand, fermionic degrees of freedom 
are antiperiodic in Euclidean time, and consequently their Matsubara frequencies are odd multiples of 
7rT, = (2n + l)vrT, n = 0, ±1, ±2, .... As T — > oo, all modes with nonzero Matsubara frequency 
become infinitely heavy, and are thus removed from the spectrum of physical excitations. These are 
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all fermionic modes, and all non-static bosonic modes. Consequently, dimensional reduction leads to a 
theory of static bosonic fields in three spatial dimensions. 

The dimensionally reduced theory can be viewed as an effective theory at energy scales much less 
than the temperature. Consider, for example, QCD in weak coupling, g < 1, where there is a distinct 
separation of energy scales, g 2 T <C gT <C T. The dimensionally reduced theory is then the effective 
theory for modes at energies of order O(gT) which one obtains from the underlying theory, i.e., QCD, by 
integrating out modes at energy scales of order 0(T). One can then take this idea one step further and 
integrate out modes at energies of order O(gT) and obtain an effective theory at an energy scale of order 
0(g 2 T). In Ref. |l()5j it was suggested to apply this principle of constructing a series of effective theories 
to compute the pressure in QCD. This task was recently carried out to order 0(g 6 ) in a succession of 
papers [100J. In the following, I outline the idea and discuss the results, which are also shown in 
Fig. EI| Note that the idea of constructing an effective theory valid on a certain energy scale has also 
been applied to non-Abelian transport theories 106J. As transport theory concerns non-equilibrium 
situations, a discussion of these aspects are beyond the scope of the present review. 

After the first step of integrating out modes at energy scales of order 0(T), the pressure in QCD 
takes the form (at \i = 0) 



where Pt{T) is the pressure of the modes at energy scales of order 0(T) and the remaining term is 
the contribution from modes at energy scales of order O(gT). The argument of the logarithm is the 
partition function of the effective theory for these modes. Since the energy scale gT is that of the 
Debye mass, mo = gT, which determines the screening length of static color-electric fields, quantities 
appearing in this partition function will be labelled with a subscript U E" . The action Se of the effective 
theory is that of a three-dimensional non-Abelian gauge theory (i.e., consisting of the color- magnetic 
fields of the original theory) coupled to a Higgs field in the adjoint representation of the gauge group 
(corresponding to the static color-electric fields of the original theory) [100 , 



Here, T { j = {i/g E )[D h Dj] = F%T a , A = ft - ig E A, and A^ = A a ^T a . There are five unknown 
quantities on the right-hand side of Eq. (|fi2~J) : the pressure pr, the mass of the adjoint Higgs field 
Aq and the coupling constants gg, X E , X E . Their values have to be determined by "matching" the 
effective theory to the original theory at some matching energy scale. At this point, however, one can 
already determine their scaling behavior from counting dimensions. The pressure pt is that of a 3+1- 
dimensional theory of modes with momenta of order 0(T). Thus, at nonzero temperature (see above) 
Pt ~ T J dkk 2 ~ T 4 . Since the action Se has to be dimensionless, one can deduce the dimensionality 
of the fields Ai from the kinetic term ~ T 2 y As lengths have dimension T~ l , the fields have to scale 
as Ai ~ T 1 / 2 . The adjoint Higgs field must have the same dimension, Aq ~ T 1//2 . From this one 
deduces that the mass term scales asm£~mo = gT, the coupling constant g E scales ~ gT 1 / 2 , and the 
four-point couplings behave as A^ ~ g A T. The dots in Eq. (|64|) denote higher- dimensional operators. 
One can show by power counting that they are not relevant if one is interested in a calculation of the 
pressure to order 0(g 6 ) |100j . 

The next step consists of integrating out modes at energy scales of order O(gT). Since physics at 
this scale is determined by static color-electric fields, or in other words, by the adjoint Higgs field Aq, 
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one has to integrate out this field, 



T 
V 



In [| VA1 VA% exp {-S E )j = p E {T) + - In [/ VA\ exp (S M ) 



(65) 



The term p E is the pressure of modes with energy of order O(gT). The argument of the logarithm on 
the right-hand side defines the partition function of an effective theory at energy scales of order 0(g 2 T). 
Since the energy scale g 2 T is associated with the scale of the magnetic screening mass m^/ ~ g 2 T in 
non-Abelian gauge theories at nonzero temperature, quantities appearing in this partition function are 
labelled with a subscript "M" . The action Sm entering the partition function of the effective theory 
at an energy scale 0(g 2 T) is simply that for a three-dimensional non-Abelian field theory for color 
magnetic fields, 



where Tij = (i/gM)[Di, Dj], Di = di — iguAi, and Ai = A°;T a . The two constants p E , gu on the 
right-hand side of Eq. (|65|) have to be determined by matching the effective theory at the energy scale 
g 2 T to that at the energy scale gT . However, their scaling behavior can already be determined by 
power counting. The pressure p E is again ~ T / dkk 2 , but now the integral runs only over modes with 
momenta of order O(gT), thus p E ~ (gT) 3 T ~ m%T. The dimensionality of the fields Ai is the same as 
in the previous effective theory, thus gM ~ gT 1 ^ 2 ~ g E . The dots in Eq. (jBTf) denote higher-dimensional 
operators which are again irrelevant if one is interested in a computation of the pressure to order 0(g 6 ). 
The final step is to compute the pressure of modes with energies of order 0(g 2 T), 



From power counting one deduces that pu ~ T Jdkk 2 ~ (g 2 T) 3 T ~ g 6 T 4 ~ gtjT, since the integral 
runs over modes with momenta of order 0(g 2 T). Due to the Linde problem, this contribution cannot 
be obtained perturbatively. What one can evaluate 100J is the contribution of order 0(g 6 In g) to pm, 
since this arises from ultraviolet divergences ~ hi{jx/mM) and not from the nonperturbative infrared 
sector which yields a genuine 0(g 6 ) contribution to the pressure. The latter has to be evaluated e.g. 
via a lattice calculation. 

The final answer for the pressure in QCD is then p(T) = Pt{T) + p E {T) + pm(T). If one wants 
to determine the pressure to order 0(g 6 \ng), one has to compute all terms appearing to this order in 
Pt, p E , and pu- This can be done perturbatively. The contributions to pt constitute a power series in 
g 2 , and not in g. They are needed explicitly only to order 0(g A ), since the full 0(g 6 ) contribution to 
the pressure is nonperturbative in any case. One then evaluates all four-loop diagrams in the effective 
theory at scales gT, in order to determine p E up to order 0(g 6 \iag) p.00 . As expected from power 
counting (see above), the lowest-order terms in p E are ~ m 3 E T ~ g 3 T A . Finally one adds everything 
to the 0(g 6 In g) term from pu- One obtains a well-defined expression for the pressure up to order 
0(g 6 lng). The term which is genuinely of order 0(g 6 ) remains unknown. The result for the pressure 
in QCD is then given by Eq. with the 0(g 6 ) contribution (iV c = 3, JV> = 0) 
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where 5 is an unknown constant. 
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In the right panel of Fig. EH the result for the pressure as a function of temperature is shown for 
various values of 5. Comparing to lattice QCD data for the pure [S77(3) c ] gauge theory, the optimum 
value appears to be ~ 0.7, since then the perturbative calculation nicely matches onto the results from 
the lattice computation, see left panel of Fig. One also observes that for the optimum value, the 
pressure up to 0(g 6 ) is rather close to the result to order 0(g 2 ), unless the temperature is very close 
to T c . This provides a certain amount of confidence that this perturbative evaluation of the pressure is 
reasonable. 

Finally, note that the above framework of constructing a sequence of effective theories via dimen- 
sional reduction was recently extended to include quark degrees of freedom at nonzero chemical potential 
|lU7j . At nonzero temperature, the quark Matsubara frequencies u f n = (2n + 1)ttT are always of or- 
der 0(T), consequently quark degrees of freedom have to integrated out in the first step in the 
construction of the sequence of effective theories. 

4.3 Quasiparticle Models 

In Sec. 14.21 the pressure of QCD was computed by evaluating the partition functions of various effective 
theories. This considerably improved the somewhat unsatisfactory situation of a purely perturbative 
evaluation of the pressure up to terms of order 0(g 5 ) as discussed in Sec. 14.11 Another way to improve 
the situation is based on the following observation. The results of Sec. 13 .41 suggest that nonperturbative 
effects still influence the physics at temperatures in the range from T c to a few times T c . Consequently, 
a perturbative expansion of the pressure around the perturbative vacuum in terms of massless quarks 
and gluons seems inappropriate. What is obviously missing in a perturbative description of the QCD 
partition function are nonperturbative effects which, when decreasing the temperature from T > T c , 
are responsible for the phase transition at T c . This was realized a while ago and attempts were made to 
incorporate them into the properties of the physical degrees of freedom. In the following, I discuss two 
such attempts, the so-called "cut-off model" and a model which treats quarks and gluons as massive 
quasiparticles. 

4.3.1 The cut-off model 

Let us again restrict the discussion to the pure [S77(3) c ] gauge theory. The cut-off model is motivated 
by the fact that QCD is an asymptotically free theory |Xj, i.e., only gluons with large momenta can 
be considered to be perturbative, while those with small momenta are subject to confinement. Quite 
similar to the effective-theory approach discussed in Sec. 14.21 one then introduces a cut-off momentum 
A to separate these two regions (1081 Il09j. Gluons with momenta larger than A are treated perturba- 
tively, while gluons with momenta smaller than A are assumed to remain bound inside colorless objects 
(glueballs in the case of pure [S77(3) c ] gauge theory). The dispersion relation for perturbative gluons 
then changes from the one for free massless particles, c^(k) = |k| = k to 



Gluons with momenta k < A are bound inside glueballs. The glueball mass scale M is of the order of 
1 GeV. The contribution of glueballs to the thermodynamic functions is then exponentially suppressed 
~ exp(— M/T), and can thus be neglected for the temperature range of interest. 



The leading-order contribution to the pressure arises from non-interacting gluons with momenta 
k > A, 



w(k) = Q(k- A)k . 



(70) 



p- t (T) = -2(iV c 2 -l)T 
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One can also compute perturbative corrections to this leading-order result. To this end, one has to 
evaluate the standard diagrams of the perturbative expansion of the pressure as discussed in subsec- 
tion |01 but with additional theta-functions like in Eq. (|70J) to restrict the phase space of the internal 
gluon lines. In Ref. |lU9j this has been done up to order 0(g 2 ). Due to the restricted phase space in 
the loop integrals, the perturbative corrections become relatively small compared to the zeroth-order 
contribution. In this sense, the perturbative series for particles with the dispersion relation ()70|) is 
better behaved than the original perturbative series when the cut-off A = 0. 

Besides the cut-off A, the cut-off model has another parameter, the MIT bag constant B which 
describes the energy difference between the perturbative and the nonperturbative vacuum. Fitting the 
parameters of the cut-off model to lattice QCD data for the pure [S77(3) c ] gauge theory, quite reasonable 
agreement could be obtained. I do not explicitly show results from Ref. |lU9j . because lattice data at 
that time were not yet extrapolated to the continuum limit. Consequently, the values for A and B 
obtained previously will quantitatively change once continuum-extrapolated data are used for the fit. 

Nevertheless, on a qualitative level, the values for A necessary to fit the data were on the order 
a typical glueball mass A ~ M ~ 1 GeV. This value for A nicely confirms the consistency of the 
assumption underlying the cut-off model, namely that gluons with momenta k < A ~ M are bound 
into glueballs of mass M. Moreover, if one interprets the cut-off model as an effective theory in the sense 
of Sec. I4.2[ the physics cannot depend on the precise value of A. In other words, if one properly matches 
the effective theory for gluons with momenta k > A to the effective theory for glueballs (i.e., gluons 
with momenta k < A), the cut-off would drop out. This matching procedure has not been done for 
the cut-off model. However, one may expect that a proper matching calculation would just confirm the 
result A ~ M obtained from the fit to lattice QCD data. One may thus simply replace the unphysical 
parameter A by the physical value of the glueball mass M to obtain a model which is independent of 
the arbitrary (and thus unphysical) cut-off scale A. The cut-off model has not been applied to lattice 
data for full QCD, because then the assumption that all colorless objects are heavy and are negligible 
when computing thermodynamic functions breaks down (the pion mass is of the order of T c ). 

The gluon dispersion relation (|70|) can be interpreted in the way that gluons with momenta below 
the cut-off momentum A have infinite mass, while those with momenta above A have zero mass. It 
is hard to believe that the true dispersion relation of gluons as computed in Yang-Mills theory would 
sustain the oversimplified and rather radical assumptions of the cut-off model. A more conservative 
model to improve our understanding of the thermodynamic properties of the QGP is explained in the 
following section. 

4.3.2 Models with massive quasiparticles 

In a hot and dense medium, particles attain a self-energy H(u, k), which (due to the breaking of Lorentz 
invariance) depends separately on energy uj and 3-momentum k, as well as on the properties of the 
medium (i.e., its temperature T and chemical potential fi). If the imaginary part of the self-energy 
on the dispersion branch of the physical excitations is not too large compared to the real part, these 
excitations are called quasiparticles The simplest situation is when the self-energy is independent 

of energy and momentum, i.e., constant and real, corresponding to a mass term, which depends on T 
and fi, but not on the energy and the momentum of the particle. In this case, the dispersion relation 
for quasiparticles of mass m reads 

w(k) = Vk 2 + m? . (72) 

The self-energies of quarks and gluons in QCD are certainly not constant (see also Sec. 14.4)1 . Never- 
theless, one can still simply assume that they are constant and explore the consequences. These so-called 
"massive quasiparticle models" have been investigated in great detail in the literature as a means to 
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describe and interpret lattice QCD data on thermodynamic functions of QCD II 13| II 14^ fTT3| HL 16j . 
The advantage of these models is that it is straightforward to extend them to nonzero quark chemical 
potential |117j . 

Here, as in the previous section let us only focus on a quasiparticle model for pure [SU(3) C ] gauge 
theory. For the generalization to QCD with dynamical quarks, see Refs. |115t lllffi I117j . For the pure 
gauge theory, there are only massive gluon degrees of freedom and the pressure reads 



m = _ DT j_ ln 



exp 



"00 ' 
T 



B , (73) 



with uj(k) given by Eq. (f72|) . As in the cut-off model, an MIT bag constant is introduced to account 
for the difference between the perturbative and the nonperturbative vacuum. 

Two comments regarding the pressure (|73|) are in order. First, while massless gluons have two 
transverse polarizations, massive gluons have an additional longitudinal polarization degree of freedom. 
Therefore, one would naively argue that the constant D in Eq. (|73|) should assume the value D = 
3 (AT? — 1), instead of D — 2(iV^ — 1), as for massless gluons. This is, however, not quite correct. 
First, the pressure ()73j) has to approach the correct Stefan-Boltzmann limit when T —>■ oo. For three 
polarization degrees of freedom, the Stefan-Boltzmann pressure is a factor 3/2 larger than the correct 
value. Another reason why this is incorrect is the following. In Sec. 14.41 we shall see that gluons indeed 
acquire a longitudinal degree of freedom in a hot or dense medium, but that the respective longitudinal 
spectral density vanishes rapidly for energies and momenta larger than gT. Thus, in a calculation of 
the pressure, which is dominated by modes with momenta of order T, one should in principle not count 
the longitudinal degrees of freedom. We shall therefore set D = 2(N 2 — 1) in the following. (Note that 
Ref. |115j also investigated a scenario where D is a function of temperature.) 

The second comment concerns the possibility to fit lattice QCD data with Eq. ([73)1 . It turns out 
that for constant mass and bag parameters, the quality of the fit is not satisfactory. Consequently, 
one needs to generalize the model (fT3*|) to allow for a temperature-dependent gluon mass, m — > m(T). 
A convenient parametrization is motivated by the dispersion relation for transverse gluons at large 



momenta, which takes the form u t (k) = Jk 2 + m 2 ^, see Eq. (I100|) . Consequently, 

m(T) = ^, ^)= llln[F(r g irc/A ^ )r F( X ,y) = K( X ) X y. (74, 

For a fixed value of T c /Aj^ — 1.03 ± 0.19 for pure [S77(3) c ] gauge theory |118j . one only needs to know 
the function K(T/T C ) in order to determine the mass, m(T), and thus the kinetic term in the pressure 
(fT3~|) . In Ref. |115j the function K(T/T C ) is simply fit to reproduce lattice QCD data. A surprisingly 
good fit is obtained with the functional form K(x) = 18/[18.4 exp(— 0.5x 2 ) + 1], see Fig. El Once 
the gluon mass is a function of temperature, thermodynamical consistency requires the bag parameter 
to depend on the temperature as well, B — > B(T). The functional form of B(T) can be uniquely 
determined from m(T); for details see Ref. |115j . 



4.4 HTL/HDL-Resummed Perturbation Theory 

Apparently, the idea that quarks and gluons are quasiparticles works rather well to describe the ther- 
modynamic functions of the QGP. Therefore, it seems appropriate to put this concept onto a more 
formal basis. In fact, the quasiparticle excitations in the QGP are well known in the weak-coupling 
limit, g<l, and for temperatures and/or chemical potentials much larger than the quasiparticles' en- 
ergies and momenta. At nonzero temperature, these quasiparticles form the basis of the so-called "hard 
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Fig. 4. SU(3), N,=0 EOS + Lattice QCD data 
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Figure 14: The temperature dependence of 3p/T 4 and e/T 4 in pure [S77(3) c ] gauge theory. 
Symbols are lattice QCD data from [7H]. The dashed lines are a fit within the massive gluon 
model. The dash-dotted lines represent the contribution of the kinetic term in Eq. (|73|) . The 
horizontal line is the Stefan-Boltzmann limit. (Since for an ultrarelativistic ideal gas e = 3p, 
this limit is the same for the functions 3p/T 4 and e/T 4 .) From Ref. |115j . 



thermal loop" (HTL-) resummation scheme [T3*l f!19j. At zero temperature, but large quark chemical 
potential, there is an equivalent approach, the so-called "hard dense loop" (HDL-) resummation scheme 
[TBI H^Uj - From the quasiparticle excitation spectrum, one can also construct the equation of state. All 
this will be discussed in detail in the following. 

4.4.1 The excitation spectrum in a hot and dense medium 

How does one determine the spectrum of physical excitations in a hot and/or dense medium? I shall 
illustrate this explicitly for the case of gluons. The case of quarks can be considered analogously, I 
only briefly report the results at the end of this section. The outline of the procedure is the following. 
First, one computes the gluon self-energy, n^(u;,p). From the self-energy, one then determines the 
full gluon propagator, A^(cj,p). From the gluon propagator, one then deduces the spectral density, 
which provides all information about the excitation spectrum in a hot and/or dense medium. Quite 
surprisingly, it turns out that one can follow this procedure in complete generality, without actually 
specifying EE^ until the very end. Although the derivation is somewhat formal, it is nevertheless a rather 
instructive exercise and will therefore be discussed in more detail |121j . Note that, in a medium, Lorentz 
symmetry is explicitly broken, and all quantities depend separately on energy u and momentum p. 
Nevertheless, to abbreviate the notation, I shall frequently use the 4-vector P M = (u, p) to characterize 
this dependence. Note also that, at nonzero temperature, one usually computes in Euclidean space 
time, i.e., all energies are discrete Matsubara frequencies on the imaginary energy axis. However, in 
order to determine the physical excitation spectrum, one has to analytically continue to real energies, 
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iu n — > u + if). This irj prescription produces the retarded Greens functions When writing real 

energies u in the following, the irj prescription will be suppressed. 

The self-energy can be decomposed according to its tensor structure. First of all, if the color 
[SU (3) c ] gauge symmetry is not broken, one may assume that the gluon self-energy is diagonal in adjoint 
color, 11^ = 5 ab IP" '. (In a color superconductor, this is in general no longer the case, for examples 
see Sec. 15.40 It thus suffices to consider IP". The next step is to decompose IP" in terms of tensors, 
multiplied by scalar functions of to and p [T3J I122| . Let us define 

= -pir- (75) 

as the projector onto the subspace parallel to P M . Then one chooses a vector orthogonal to P M , for 
instance 

N^{^,^f) = (<r-vnu, (76) 

with = (0, p). Now one defines the tensors 

With the help of these tensors one can write the gluon self-energy as 

IP" = IP A^" + II b B^" + IT CP" + n e W u . (78) 

The scalar functions n a,b,c ' e can be obtained by suitable projections of IP" onto the respective tensor 
structures. 

Using the explicit form of A^ 4 , one convinces oneself that the tensor A M " projects onto the spatially 
transverse subspace orthogonal to P^, 

A 00 = A 0i = , A ij = - (6 ij - ffP) . (79) 

This means that the self-energy function IP determines the excitation spectrum of the spatially trans- 
verse gluon fields 

A ± «(P) = A»Al(P) . (80) 

As A M " projects onto a two-dimensional subspace, there are two degrees of freedom associated with 
A±j 1 . In the vacuum, these are the only physical degrees of freedom, since gluons are massless. 
The tensor IP" projects onto the spatially longitudinal subspace orthogonal to P M , 

B 00 = -^ , B- = -^ , B« = -^ff. (81) 

Consequently, the polarization function IP determines the excitation spectrum of the longitudinal gluon 
degree of freedom, 

A°(P) 

A%(P) = jfi } , (82) 

which becomes a physical degree of freedom in a medium. 

The spatially transverse and spatially longitudinal gluon fields A^_ a and A a N are the only physical 
degrees of freedom. There is an unphysical degree of freedom associated with the projection onto P M , 

pn j\a ( p) 

Al(P) = p t ' ■ (83) 
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The scalar function IT 3 is the self-energy of this unphysical degree of freedom. It will be seen that gauge 
fixing ultimately removes this degree of freedom from the theory. A nonvanishing IF indicates that 
the spatially longitudinal, physical gluon degree of freedom mixes with the unphysical degree of 
freedom A?,. Before extracting the physical excitation spectrum, one has to remove this mixing term, 
as will be discussed below. 

Now use the tensor decomposition ([75)1 to determine the full gluon propagator Since the free 
inverse gluon propagator 



A" 1 ^ = 5 ab (P 2 g^ - P^P") ee 5 ab P 2 (A^ + EH 



34) 



is diagonal in adjoint colors, so is the full inverse gluon propagator, A 



ab 



5abA 



One obtains 



a- 1 ^ ee a uw + if" ee (p 2 + if) a"" + (p 2 + n b ) b"" + n c cr + n c w u . 



35) 



In an effective action, this inverse propagator is the coefficient of the term quadratic in the gauge fields. 
In momentum space (and choosing a normalization such that A a AP) retains dimensions of energy), 



S-2 



1 V 



N?-l 



-o t E E Al^-P) [P 2 + IF(P)j AT A\ V {P) - A a N (-P) [P 2 + U h (P)\ N 2 A%{P) 

Z 1 p a =l 1 

- A\{-P) n c (P) N 2 P 2 A a N (P) - A a N (-P) n c (P) N 2 P 2 A\{P) 

- A^{-P)U C {P)P 2 A^{P) | , 



36) 



where J2p = (Y/T) Jp- The physical excitation spectrum can be most easily extracted from a diago- 
nalized inverse gluon propagator, i.e., the term which mixes the physical field component A% with the 
unphysical component A^ has to be eliminated. This can be done as follows. Remember that in the 
partition function one functionally integrates the exponential of the action (a part of which is S2 in 
Eq. (|86Jl ) over all gauge field components, including Af,. For the purpose of diagonalizing the inverse 
gluon propagator one may simply redefine the integration variable 



W(P)N 2 

A;(P) - Al (P) = A\{P) + A-L A%(P) . 



37) 



This redefinition does not change the physics (as it only involves the unphysical component of the gauge 
field) and diagonalizes the action S2 in the components of the gauge field, 



So 



1 V 

2 f 



N?-l 



P a=l 



E E \ A ±^- p ) P 2 + R\P) A"" A a ±u (P) - A a N {-P) P 2 + n b (P) N 2 A a N (P) 



where 



Afd-P) n c (P) P 2 AUP) 



n b (p) ee n b (P) 



[W{P)] Z N 2 P 2 

n c (p) 



39) 



From Eq. (|88|l one can read off the inverse gluon propagator in diagonal form. However, in order to be 
able to invert it, it is necessary to fix the gauge. (To see this, consider the case where IF = 0. Then, the 
inverse propagator has a zero eigenvalue and is not invertible). For the sake of simplicity, let us choose 
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Figure 15: The one-loop diagrams contributing to the gluon self-energy (possible ghost 
contributions are not shown). 



covariant gauge, where the gauge fixing term in the action only involves the unphysical components of 
the gauge field, 

JV?-1 



S, 



gf 



1 V 
2XT 



E E A?,(-P) P 2 P 2 AMP) . 



P a=l 



Adding to S 2 , one reads off the (gauge- fixed) inverse gluon propagator 

1 



A-^(P) = P 2 + IT(P) 



A 



p 2 + n b (P) 



B 



-P 2 + IT(P) 
A 



which can be straightforwardly inverted, since A 11 ", B Mi/ , and E Ml/ are projectors, 

~ 2 A 



A^(P) = A t (P) AF» - MP)^ + p2 + x ne(p) 

where the transverse and longitudinal propagators are defined as 

11 P 2 1 



MP) 



p 2 + n a (P) p 2 - n t (P) 



p 2 p 2 + n b( P ) p 2 -n^(p) 



(90) 



(91) 



(92) 



(93) 



with the transverse and longitudinal polarization functions Ut = — Il a and = — (p 2 /P 2 )II b . The 
last term in Eq. ()92j) can be removed from the spectrum of physical excitations by the gauge choice 
A = 0. The physical excitations are described by the transverse and longitudinal propagators defined 
in Eq. (J9"3|) . From these one can derive the corresponding spectral densities in the usual way [TH] 



Pt(v,p) = -lmA t (u + i7],p) 
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(u, p) = - Im A e (u + ir], p) . 



7T 



(94) 



All that is left is to actually specify the form of the polarization functions Ut and IL\ In the weak- 
coupling limit, g 1, this can be done via a one-loop calculation of these functions. To this end, one 
has to compute the diagrams shown in Fig. using standard methods fOJEI]. Note, however, that 
the result depends in general on the choice of gauge and thus cannot determine the physical excitation 
spectrum, which is by definition independent of the choice of gauge. However, it was noticed many 
years ago |123j . that the high-temperature limit T 3> u,p of the polarization functions is actually 
independent of the choice of gauge, 



Re U t (u,p) 
Im U t (u,p) 



3 2 
- m n 
2 9 



UJ 2 / u 2 \ UJ , 

t + 1 - t hr ln 



oo + p 



LU — p 



UJ 



1 - 



UJ 



p 



p 



e{p-uj) 



(95) 
(96) 
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ReILj(a>, p) = — 3m g ( 1 — — In 



2p 



uj + p 



LU — p 



(97) 



3 n 00 



ImIL(w,p) = -rt^m 2 g -Q(p-uj) , (98) 
where the gluon mass parameter at a given T and /z is 

V, _1 

(99) 



2 2 

m g = g 



2N C + Nf 2 Nj_ 2 



18 ' 6tt 2 

This result is not restricted to the high-temperature limit, in fact it holds as long as either T or /i is 
much larger than u and p. Thus, it also describes gluonic quasiparticle excitations at T = and high 
density. It is only this gauge-invariant high-T (or high-/i) limit of the one-loop polarization functions, 
which is relevant for the HTL- (or HDL-) resummation scheme discussed in Sec. 14.4.21 

What is the physical meaning of the result (j93j) - ([EfH])? The condition T, \i 3> u, p implies that there 
is a separation of scales, just like in the construction of the effective theories in Sec. 14.21 The temperature 
(or the chemical potential) sets a "hard" energy scale, while the external energy and momentum are 
"soft". As will be seen in Sec. 14.4.21 this separation of scales forms the basis of the HTL- (or HDL-) 
resummation scheme. In this scheme, energies of order gT (or g/j,, in the HDL-resummation scheme) 
are "soft", while energies of order T (or /z, respectively) are "hard". In the loops in Fig. one 
integrates over the internal loop momentum, k, but the result must be finite, as there are no other 
ultraviolet divergences at nonzero T and/or /x than those already known from the vacuum. Therefore, 
the ultraviolet regularization must be provided by the distribution functions of quarks and gluons, 
which decrease exponentially with temperature. (At T = and /i ^ 0, the gluon distribution function 
vanishes, while the quark distribution function is a step function ~ 0(/x — k), which cuts off momenta 
k > /i). On the other hand, the phase space in the loop integral grows ~ / dkk 2 . One thus expects 
that the dominant contribution to the loop integral comes from "hard" momenta of order T (or, at 
T = and fi ^ 0, from momenta close to the Fermi surface, k ~ fi). For dimensional reasons, and 
including factors of g from the vertices, II ~ g 2 T 2 (or ~ g 2 f- 2 , at T = and \i ^ 0). This gives rise to the 
prefactor ~ m 2 in Eqs. (j93j) - ()98|). The "soft" external energy and momentum cannot significantly alter 
the kinematics in the loop, where the dominant contribution comes from "hard" momenta. In fact, it 
suffices to expand the integrands of the loop integrals to leading order in these external quantities. This 
gives rise to the particular dependence on u and p of the result (|95|) - (|98|). The essential approximation 
which leads to this result is the assumption that internal momenta are exclusively "hard". Therefore, 
the loops computed under this assumption are called "hard thermal loops' 1 (or "hard dense loops" , at 
T = and /i ^ 0). 

From the self-energies one can construct the propagators (J9l?j) and the spectral densities (JMj) . respec- 
tively. Following Ref. |124j it turns out that the propagators have poles above the light-cone, u > p, 
and a cut below, u < p. In the spectral densities, the poles become 5-functions. These determine the 
excitation branches u){p). Since a 5-function has no width, the quasiparticles corresponding to these 
excitation branches are stable, i.e., they have an infinite lifetime. (This changes if one computes beyond 
one-loop order.) The excitation branches are above the light-cone, i.e., they correspond to time-like, 
propagating gluons. In the left panel of Fig. they are shown for transverse and longitudinal gluon 
modes. For large momenta, the longitudinal mode has an exponentially vanishing residue |124j . In 
contrast, the transverse mode has a finite residue and a dispersion relation which approaches the form 

3 



w *(p) -* yp 2 + m too , ™ 2 oc = 2 m ? • ( 10 °) 

The cuts in the propagator become continuous distributions in the spectral density. They provide 
Landau-damping for space-like gluons. 
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Figure 16: The excitation branches for gluons (left panel) and quarks (right panel) in a hot 
and/or dense medium. For gluons, energy and momenta are shown in units of the gluon 
mass parameter, m g . For quarks, energy and momenta are shown in units of the fermion 
mass parameter, rrif. 



Performing a similar exercise for quarks 125J, one obtains the spectrum of fermionic quasiparticle 
excitations. It turns out that there are twice as many excitation branches as expected. There are in 
fact two solutions for positive energies, uj±(p), determined by the equation 



uj± (p) = ±p ± — - 
p 



<4t(p) TP , ( w±(p) +p s 



2p \v±(p)-p 
where the fermionic mass parameter (squared) is 

N 2 



(101) 



TO *=^wr + ^J- (102) 

These two solutions are shown in the right panel of Fig. (The two solutions for negative energies 
mirror the above solutions below the p axis.) The solution u)-(p) corresponds to a quasiparticle with 
the opposite chirality than the one associated with the solution u> + (p). This peculiar quasiparticle is 
commonly called the "plasmino". Note that, while the ordinary quasiparticle dispersion branch uj + (p) 
has a positive group velocity du + (p)/dp > 0, the plasmino branch has negative group velocity at small 
momenta. For large momenta, the residue of the plasmino branch becomes exponentially small, while 
the one for the ordinary quasiparticle remains finite and its dispersion relation approaches the form 



w+(p) - \ V 2 + m 2 foo , m 2 foo = 2 m 2 . (103) 



The quark propagator also features a cut below the light-cone, which gives rise to Landau damping. 
Finally, note that the fermionic quasiparticle spectrum in Fig. is shown for energies and momenta 
much smaller than either temperature and/or chemical potential. However, at T = and large fi, this 
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kinematic region is irrelevant, as it reflects the situation at the bottom of the Fermi sea. At T = 
and large //, the relevant fermionic excitations are those around the Fermi surface, where p ~ //. The 
excitation spectrum for this case will be discussed in more detail in Sec. 15.21 



4.4.2 HTL/HDL-resummation scheme 

The HTL-resummation scheme and its counterpart at T = and nonzero /i, the HDL-resummation 
scheme, are explained in great detail in textbooks on field theory at nonzero temperature and density 
(see, for instance, Ref. Therefore, for more details, I refer the interested reader to the literature 

and restrict myself here to a short discussion of the general idea behind these methods. 

As already discussed in Sec. 14.11 due to the infrared behavior of gauge theories at nonzero tempera- 
ture, naive perturbation theory breaks down. We have seen in Sec. !4.4.Tl that loop calculations involve 
propagators of the form 

P) ~ ~2 2~TW( V ' ( 104 ) 

The leading-order terms in the one-loop self-energy arise from "hard" particles with momenta k ~ T 
inside the loop. Together with factors of the coupling constant arising from the vertices, the self-energy 
is II ~ g 2 T 2 . Therefore, as long as either u or p are "hard", i.e., of order T, the self-energy II in the 
propagator (jl(J4j) can be neglected. However, for "soft" u and p of order gT, the self-energy is of the 
same order of magnitude as the first two terms in the denominator and cannot be neglected. 

This observation forms the basis of the so-called HTL-resummation scheme fSJ I119j in field theories 
at nonzero temperature. In simple words it states that whenever the energy and the momentum of 
a propagator in a given diagram is "soft", one has to use the "dressed" propagator (| 11)4)1 including 
the self-energy n, and if either energy or momentum is "hard", one may use the "bare" propagator 
Aq(u,p) = 1/(uj 2 — p 2 ) without the self-energy n. How does the name "HTL-resummation scheme" 
arise? The dressed propagator is the solution of the Dyson-Schwinger equation 

A = A - A n A . (105) 



Iterating this equation, one realizes that it stands for an infinite series of diagrams; consequently the 
solution ()104)) is a resummed propagator. As explained in Sec. 14.4.11 the quantity which is resummed 
is the self-energy n computed in the HTL approximation, i.e., one resums HTL's. 

It is now also easy to see why naive perturbation theory breaks down at nonzero temperature. 
Imagine a diagram with n vertices, such that naive perturbation theory would tell us that this diagram 
is of order 0(g n ). Now imagine that there is a loop in this diagram with propagators of the type (|1U4)) 
and that the dominant contribution to this loop arises not from the "hard" region of phase space, i.e., 
from momenta of order T, as in HTL's, but from the "soft" region, i.e., from momenta k ~ gT. The 
contribution of the propagator (jl(J4j) to the diagram is then ~ l/(g 2 T 2 ), instead of ~ 1/T 2 . This cancels 
two powers of the coupling constant in the naive perturbative counting scheme. The diagram is thus 
actually of order 0(g n ~ 2 ) (or even of lower order, if other propagators contribute additional powers 
of g~ 2 ). The occurrence of the additional energy scale T (or /z, at T = 0) compared to the vacuum 
invalidates the naive perturbative counting scheme. 

For gauge theories, the Ward identities require to extend the HTL-resummation scheme from propa- 
gators to vertices as well. For physical quantities which are determined by computing diagrams with at 
least one loop, depending on whether the dominant contribution in the loop arises from the hard or the 
soft region of phase space, one may be required to use only bare propagators, resummed propagators, 
or both resummed propagators and resummed vertices. The major success of the HTL-resummation 
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scheme was the proof that the leading-order result for the gluon-damping rate is independent of the 
gauge and positive [126J. 



4.4.3 The equation of state for quasiparticles in HTL/HDL approximation 

The objective of this section is to apply our knowledge from Sec. 14.4.11 about the quasiparticle excitation 
spectrum in QCD to determine the equation of state at high temperature and/or density. What one 
obviously needs is a thermodynamically consistent way to construct the pressure including information 
about the spectral density of the quasiparticles. Obviously, classical statistical mechanics, such as 
applied in Sec. 14.3.21 is of no use; we need a field-theoretical approach. The method of choice is 
the so-called "Cornwall- Jackiw-Tomboulis" (CJT) formalism |127j . The CJT formalism determines the 
effective action of a theory as a functional of the one- and and two-point functions. The stationary values 
of the effective action yield the expectation value of the field and the full propagator. The stationarity 
conditions are Dyson-Schwinger equations for these quantities. The CJT formalism is particularly useful 
for theories with spontaneously broken symmetries, see Sec. 14.61 For unbroken symmetries, the CJT 
formalism is equivalent to the so-called ^-functional approach [128J. 

Let us elaborate on this in somewhat greater detail. In the CJT formalism the effective action of 
a theory with bosonic fields <fi and corresponding propagators A, as well as fermionic fields ■?/>, ip with 
propagators G reads 

VMMAM = lMM~\^ lnA^-^Tr (iT 1 A - l) 

+ Tr lnG _1 + Tr (S^G-l) + r 2 [0, <0, <0, A, G] • (106) 



Here, is the classical action and all traces are taken in the functional sense. The quantities 

D^ 1 and S 1-1 are the inverse tree-level propagators for bosons and fermions, respectively, 

The functional T 2 is the sum of all two-particle irreducible (2PI) diagrams without external legs and 
with internal lines given by the propagators A and G. The stationarity conditions which determine the 
expectation value of the bosonic field, ip, as well as the full propagators for bosons, T>, and for fermions, 
Q, read 

5T[(f>, ^, A, G] 5T[<p, if,, A, G] 5r[<j>, ^ A, G] 



5<f>{X) ~ 5A(X,Y) ~ 5G(X,Y) 



(108) 



In principle, there are also stationarity conditions for the fermionic fields ip and ip, but their solution is 
always trivial, as fermionic fields are Grassmann-valued and thus cannot assume a nonzero expectation 
value. 

Inserting the explicit form of T from Eq. ()1U6|) into the last two equations of ()108|) . one obtains the 
Dyson-Schwinger equations for the full propagators for bosons, T>, and fermions, Q, 

D- 1 = IT 1 + II , ^ 1 = 5'- 1 + E, (109) 

where the bosonic and fermionic self-energies are 

n _ 2 5T 2 [^^,^A,G} s _ 5T 2 M^i^G] 
5 A 5G 
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Figure 17: The two-loop approximation to T 2 (possible ghost contributions are not shown). 



The right-hand sides of these equations have to be taken at the stationary point = ip, ip = ip = 
0, A = £>, G = Q. According to their definition (jllOj) . the self-energies are obtained from the set of 
2PI-diagrams T 2 by opening one internal line. 

For translationally invariant systems, (f)(X) = <p — const., A(X,Y) = A(X — Y), G(X,Y) = 
G(X — Y), and it is advantageous to work in energy- momentum space instead of in space-time. The 
effective action is, up to a sign and a factor of the four-dimensional volume of the system, V/T, equal 
to the effective potential, 

^[0,^-0,A,G] = -|r[0,^V,A,G] . (Ill) 

At the stationary point, the effective potential is, again up to a sign, equal to the thermodynamic 
pressure. Utilizing the stationarity conditions (jl()9j) and the definitions (jllOj) . 

p=-V[ip,V,g] = -U{<p) - ^Tr lnP^ + ^TrnP + Tr ln^T 1 - Tr £ Q - V 2 [<p, V, Q\ . (112) 

Here, U(cp) is the tree-level potential, V2 = — (T/V)r 2 , and Ti A = J K tr A(u n ,k), where tr runs over 
possible internal indices (Lorentz, Dirac, color, flavor, etc.) of A. 

Equation ()112j) is the desired result. All one has to do is to apply it to QCD. In QCD, there is 
no spontaneously broken symmetry (unless one considers color superconductivity, see Sec. EJ), and thus 
U((f) = 0. The bosons are the gluons and the fermions are the quarks. Of course, one also has to 
account for ghost degrees of freedom. These look like another fermion contribution in Eq. (|112|) . but 
the corresponding Matsubara sum in Tr has to run over even multiples of ttT [TU] . 

This sounds much simpler than it is to realize in practice. The difficulty obviously lies in solving the 
Dyson-Schwinger equations for the gluon and quark propagators. It is clear that, as V% consists of an 
infinite set of diagrams, one can never aspire to solve the problem exactly. However, a great advantage 
of the CJT formalism is that the solutions of the Dyson-Schwinger equations (jl(J9j) are self-consistent 
and conserving, even if one truncates this infinite set. Any truncation defines a meaningful many-body 
approximation scheme (for a well-known example, the Hartree approximation, see Sec. 14.6)1 . Let us 
therefore imagine we only take a finite subset of all diagrams in V2, for instance the one consisting of 
the two-loop diagrams shown in Fig. El This subset is particularly interesting, because when computing 
the self-energies according to Eq. ()110|) one obtains these self-energies to one-loop order. As discussed 
in Sec. 14.4.11 the high-temperature (high-density) limit of the one-loop self-energies defines the HTL 
(HDL) approximation, which is already known to provide a physically interesting, meaningful, and 
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gauge-invariant quasiparticle excitation spectrum. Another advantage of restricting V2 to a set of two- 
loop diagrams is that the entropy density s = dp/dT and the quark number density n = dp/dfi assume 
the particularly simple form [129] 

lining" 1 - Im£Re£ ,(113) 



s = — tr 
n = — 2tr 



d 4 K dn(uj) r T , ^ 1 T i n /■ d 4 if 0/fo;) 



ImlnP -1 - ImlTReD 



- 2tr 



(2tt) 4 <9T L J J (2vr) 4 <9T 

d 4 if 9/(o;) 



(2vr) 4 <9^ 



Imln^ 1 - Im£Re£ , (114) 



where = (u;,k) and n(u;) = (e ul / T — 1) _1 is the Bose-Einstein distribution function, while f(co) = 
[e( u ~M' T + 1] _1 is the Fermi-Dirac distribution function. While one has to compute two-loop diagrams 
to obtain the pressure ([112)1 . the entropy and quark number densities ()113I114|) are essentially one- loop 
quantities and thus much simpler to calculate. 

However, even when restricting V2 to the simple subset of Fig. El the solution of the Dyson- 
Schwinger equations is highly nontrivial (see also Sec. 14. 6|) . Let us suppose, however, that it can be 
achieved. Then there is still the issue whether the propagators thus obtained obey the Ward identities. 
In general, this is not the case [12QJ , and a cure of this problem would most likely require a self-consistent 
calculation of the three- and four-point functions on top of the two-point functions. An extension of the 
CJT formalism in this direction was proposed in Ref. |131j . but since that work considers only scalar 
field theories, it is not clear whether this approach also applies to gauge theories in a way which preserves 
the Ward identities. Another question is the gauge invariance of the solution. Fortunately, as shown in 
Ref. |13(Jj . the gauge dependence always enters at an order which is higher than the truncation order. 
Finally, there is the issue of renormalizability. For scalar theories, renormalizability was demonstrated 
in Ref. |132j . For gauge theories, it still remains an open question. 

In order to make progress, one has to simplify the solution of the problem. Instead of a completely 
self-consistent solution of the Dyson- Schwinger equations f!109|) . as a first educated guess one might 
simply use the self-energies in the HTL (or HDL) approximation to compute the propagators entering 
the pressure (1112)) . Since this approximation is independent of the choice of gauge, gauge invariance 
is not an issue anymore. There is also another advantage: since the HTL- (HDL-) approximated self- 
energies are just the high-T (high-//) limit of the one-loop self-energies, Eqs. ()113I114|) for the entropy 
and quark number densities apply. These expressions are completely ultraviolet safe and thus do not 
require normalization. They are one-loop quantities and thus simpler to compute than the pressure 
which contains two-loop diagrams. From the entropy and quark number densities one can always 
deduce the pressure via an integration with respect to T and \i. 

This approach has been followed in Ref. [133J , where a detailed discussion of the results and their 
interpretation can be found. In particular, since the self-consistent approach is based on a partial 
resummation of a subset of perturbative diagrams, i.e., a reorganization of the perturbative series, it 
is possible to compare with standard perturbation theory (Sec. 14.1)1 by re-expanding the results in 
powers of g; for more details, see Ref. [T71 1133) . Here, I restrict myself to a discussion of one of the 
main results, namely the pressure for pure [S77(3)] c gauge theory as computed in this approach in 
comparison to lattice QCD data, cf. Fig. HH1 One observes that the agreement at large temperatures is 
quite good, and that deviations occur only below T ~ 2.5 T c . The quality of reproducing lattice QCD 
data is comparable to that in the perturbative approach of Sec. 14.21 However, a distinct advantage of 
the self-consistent approximation scheme is that the physical interpretation of the results in terms of a 
gas of quasiparticles is much more appealing. 

Note that a very similar approach has been pursued in Ref. |134j . utilizing directly the expression 
(1112)1 for the pressure, with propagators in the HTL approximation and a suitably modified V2. Another 
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Figure 18: The pressure, normalized to the Stefan-Boltzmann limit, as a function of temper- 
ature, normalized to the critical temperature. Different resummation methods are compared 
to lattice QCD results for pure [SV^c] gauge theory. The band labelled "pert.th. to order 
g 3 " is the order 0(g 3 ) perturbative result computed in Sec. 14.11 and shown in Fig. ED The 



band labelled "ABPS" is the two-loop HTL-resummed pressure from Ref. [135J. The band 
labelled "BIR" is the result from Ref. |133| . Bands arise from varying the renormalization 
scale ft within certain limits. Lattice data from Refs. [7H1 I136| are shown as grey bands. 



related approach is so-called "HTL-resummed perturbation theory" |135j . Here, the pressure of the pure 
gauge theory is computed to two-loop order with HTL-resummed propagators. A variational procedure 
like that of Ref. |137j is applied to determine the gluon mass parameter. The pressure computed from 
this approach is closer to the Stefan-Boltzmann limit and fails to agree with lattice QCD data in the 
region of temperatures from T c to several times T c , cf. Fig.^3 The reason for the differences between this 
approach and the one of Ref. |133j lies in the way the perturbative series is reorganized when resumming 
certain subsets of diagrams. The differences can be explained by comparing to a dimensionally reduced 
version of HTL-resummed perturbation theory |138j . Finally, note that the CJT formalism |133j as well 
as HTL-resummed perturbation theory |139j were also applied to compute thermodynamic properties 
of the QGP at T = and nonzero //. 



4.5 Polyakov-Loop Model 

The pressure, normalized to the appropriate Stefan-Boltzmann limit, as a function of temperature, 
normalized to the appropriate critical temperature, shows a universal behavior, see right panel of 
Fig. El Since the normalized pressure for full QCD with dynamical fermions looks the same as for the 
pure gauge theory (i.e., without dynamical fermions), a natural conclusion would be to assume that 
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the dynamics of the gluons drives the QCD phase transition, not that of the fermions. Consequently, 
the order parameter for the transition in full QCD should be the same as in the pure gauge theory, i.e., 
the Polyakov loop, see Sec. 12.3.21 For an [SU(N C ) C ] gauge symmetry, the Polyakov loop (J7J) is invariant 
under [SU(N C ) C ] gauge transformations, up to an element of the center of the gauge group, Z(N C ), 
L(x) — > exp(2irin/N c ) L(x). The effective theory for the Polyakov loop consists of all possible terms 
invariant under Z(N C ) transformations [70J, 
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C cS = c | VL\ A + c 2 \L\ Z + c 3 [L 3 + (L*Y\ - c 4 [\L\ 2 ) + . . . , (115) 

where the dots denote higher-order terms (which will be neglected in the following). For N c = 2, one 
would have C3 = 0, because L 3 is not Z(2)-symmetric. For N c = 3, it is precisely this cubic term which 
drives the transition first order [2E], see also Sec. 12.3.21 

The ground state of the theory, Lq = (L(x)) (which is assumed to be real), is determined by the 
global minimum of the effective potential 

V eS (L) = -ca \L\ 2 - c 3 [l 3 + (L*) 3 ] + c 4 (|L| 2 ) 2 . (116) 

The Polyakov loop L is dimensionless, therefore all coupling constants in Eq. (Jllfcij) carry dimension 
[energy] 4 . Since the only dimensionful scale is set by the temperature, one may pull out a factor T 4 
from the right-hand side of Eq. ()116|) and after appropriately renaming the constants write 



V cS (L)=hT^-b 2 \L\ 2 -b 3 [L 3 + (L*y\ +(\L\') I , (117) 

where the new coupling constants are dimensionless. In order to have a potential that is bounded from 
below for large values of L, one has to assume 64 > 0. For the moment take 63 = (for N c = 2, 
63 = C3/(&4T 4 ) must always vanish). Then the sign of b 2 drives the transition. Below the transition, 
T < T c , there is confinement, which can only be achieved if the curvature of V c q at the origin is 
positive, b 2 < 0, such that the global minimum is at Lq = 0. Above the transition, T > T c , there is 
deconfinement, which is achieved by a negative curvature of V e s at the origin, b 2 > 0, giving rise to a 



nonvanishing Lq = ±y& 2 /2- The constant b 2 is therefore a function of temperature, b 2 = b 2 (T), while 
in the simplest version of the Polyakov- loop model, 63 and 64 are assumed to be constant. Through the 
temperature dependence of b 2 , L also becomes a function of temperature, L = L (T). At T = T c , 
one must have b 2 (T c ) = 0, i.e., L vanishes continuously at T c , and one has a second-order phase 
transition. A nonvanishing 63 (which is allowed for N c = 3) turns this second-order transition into a 
first-order transition. Lattice QCD data indicate that the transition for pure [SU(3) C ] gauge theory has 
a comparatively small latent heat, i.e., the transition is only weakly of first order. This suggests that 
the constant 63 is small. 

Now remember that the pressure is, up to sign, equal to the value of the effective potential at its 
minimum 

p _ V eS (L ) b A 2 

fi = fT~ = ^ L o ( fc 2 + h L ) . (118) 

For iV c = 2, where 63 = 0, one has p/T 4 = b^b 2 {T) 2 / '4, where the explicit value of L was used. The 
physical interpretation of this result is quite astonishing: the pressure in QCD, as calculated on the 
lattice, cf. Fig. is not determined by the kinetic energy of a gas of weakly interacting quasiparticles, 
as advocated in Sees. 14 .3| 14 A\ but is simply given by the potential for the Polyakov loop (at the 
respective minimum). All that one has to know is Lq(T) (or, equivalently, b 2 (T), under the assumption 
that this is the only coupling constant which depends on temperature), and the pressure in QCD 
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Figure 19: The functions 3p/T A (lower curve) and e/T 4 (upper curve) as a function of T/T c . 
Fit of the thermodynamic functions of the Polyakov-loop model to lattice data (shown as 
circles and boxes) for the pure [S77(3) c ] gauge theory. 



can be immediately computed from Eq. (jll8|) . Comparing Figs. El and this conjecture is at least 
qualitatively correct. Quantitatively, it is not as simple: not only is the absolute value of L as calculated 
on the lattice subject to renormalization, also the constants b 2 ,bs, and 64 are not known. In the 
meantime, one can at least convince oneself that there are no principal obstacles for this interpretation: 
one may insert the value of Lq (which is only a function of b 2 and 63) into Eq. (jll8p and simply 
fit 63, 64, as well as b 2 (T) to lattice QCD data. The result of this fit is 63 = 4/3, 64 ~ 0.1515, 
b 2 (x) = 2(1 - l.ll/x)(l + 0.265/x) 2 (l + 0.300/x) 3 - 0.974, where x = T/T c . The pressure and the 
energy density obtained from this fit are shown in Fig. IT91 

The Polyakov-loop model also allows to predict other physical observables, for instance, the ratio of 
screening masses related to the correlation function of the real and the imaginary part of the Polyakov 
loop [70]. Decompose the Polyakov loop into its real and imaginary part, L = R + il. In analogy to 
Eq. P3*|) . the correlation functions for the real and imaginary part are then defined as (-R(O) -R(x)) and 
(1(0) /(x)), respectively In the deconfined phase, one expects them to decrease exponentially with the 
distance, 

(R(0) R{x)) ~ exp (-m fl |x|) , (7(0) J(x)> ~ exp (-m,|x|) , (119) 

where wir, mi are the corresponing screening masses. In weak coupling, one can make a definite 
prediction for these masses. To this end, expand L in powers of g to leading nontrivial order. One 
obtains R ~ 1 — g 2 /(T 2 N c ) Tr„4.Q. Thus, the correlation function for the real part involves the exchange 
of two static Ao fields, which are Debye screened. One therefore expects (R(0) -R(x)) ~ exp(— 2m D |x|). 
Analogously, (1(0) /(x)) ~ exp(— 3m D \it\). Thus, the perturbative prediction for the screening masses 
is m R = 2m D , mj = 3m D , and mi/m R = 3/2. 

In general, because of the renormalization of the Polyakov loop it is not as simple to obtain an 
answer for the absolute values of m R and mj. However, in the ratio mijm R one expects unknown 
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factors to cancel out [7U]. The Polyakov-loop model makes a definite prediction for this ratio. First, 
write the potential ()116|) in terms of R and /. As usual, the curvature of the effective potential in the 
ground state provides the masses. Consequently, computing the curvature in R- and in /-direction, one 
obtains 

m|~ -b 2 -6b 3 L + GL 2 , mj ~ -b 2 + 6b 3 L + 2L 2 . (120) 

I have refrained from writing the constant of proportionality, which also provides the correct dimension 
for rriRj, because this constant will drop out anyway in the ratio mijm^. Let us compute the ratio 
close to T c . In the case of a first-order transition, at T = T c the effective potential V e g has two 
minima, one at the origin and one at Lq ^ 0, which are degenerate and separated by an energy barrier. 
The condition V^r(O) = V e s(L ) gives L = —62/^3, which together with the condition that L is the 
nontrivial minimum of V e s, Lq = (3b^Lo + b 2 )/2, allows to express all quantities in Eq. (jl20)l in terms 
of, say b 3 L . This yields mj/m^ = 3. This value differs from the perturbative expectation by a factor 
of two, and is a definite prediction of the Polyakov-loop model which can be tested, for instance, on the 
lattice. 

Another particularly appealing aspect of the Polyakov-loop model is that it correctly describes the 
behavior of the screening masses when approaching T c from above, cf. Fig. [7| Near T c , the temperature 
dependence of b 2 and L cause the screening masses m R , mj to decrease, in agreement with lattice QCD 
data. The fact that the perturbative prediction completely fails to reproduce this behavior was already 
mentioned in Sec. 13.41 

Finally, let us remark that the Polyakov-loop model may also naturally explain why the chiral 
symmetry restoration transition occurs at the same temperature as the deconfinement transition. There 
is nothing to prevent the effective theory (|115j) from having a term ~ c|L| 2 Tr <3>t<3>, where $ is the chiral 
condensate and c > a constant. Thus, when \L\ condenses (in the deconfined phase), this term 
suppresses condensation of $ on account of the positive coupling constant c, effectively restoring the 
chiral symmetry. Thus, chiral symmetry restoration and deconfinement occur at the same temperature. 

4.6 Linear Sigma Models with Hadronic Degrees of Freedom 

In Sec. 13. 21 we have seen that the QCD transition is most likely crossover for small values of the chemical 
potential fi. Therefore, there is no real distinction between hadronic degrees of freedom on the one side 
and quark and gluon degrees of freedom on the other side of the transition. Consequently, there is no 
compelling reason why the thermodynamic properties of QCD around the phase transition should be 
described in terms of quarks and gluons; a hadronic description should be equally adequate. (For a 
more complete discussion of this argument, see Ref. |14U| .) 

Let us assume that, somewhat contrary to the idea behind the Polyakov-loop model discussed in 
Sec. 14.51 the QCD transition is driven by chiral symmetry restoration instead of deconfinement. In this 
case, the order parameter for the transition is the chiral condensate, see Sec. 12.3.31 and the effective 
theory is given by the linear sigma model (jl3j) . The degrees of freedom in this effective theory are the 
fluctuations of the order parameter field $ around the ground state, (<&). Physically, these fluctuations 
correspond to the scalar and pseudoscalar mesons. Following the above line of arguments, this effective 
theory could therefore not only serve to understand the dynamics of chiral symmetry restoration, it could 
equally well be used to describe the thermodynamic properties of QCD around the phase transition. 

Of particular interest is the question how chiral symmetry restoration is exhibited in the meson 
mass spectrum. In the framework of lattice QCD, this was discussed in Sec. 13.51 Here, I explain how 
to answer this question using the effective theory (fT3*j) at nonzero temperature. The standard definition 
of a particle mass is via the pole of the propagator at zero momentum. Consequently, the goal is to 
determine the mesonic propagators. The method of choice is obviously the CJT formalism discussed 
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Figure 20: Double-bubble diagrams (left-hand side), where full lines denote scalar mesons 
and dashed lines denote pseudoscalar mesons. Cutting the bubbles produces the tadpole- 
diagram contributions to the self-energies shown on the right-hand side. The tadpole dia- 
grams constitute the so-called Hartree approximation. 



in Sec. 14.4.31 because this approach allows to compute the full propagator from the stationary points 
of the effective action T, cf. Eq. (jl()8j) . or in other words, from the Dyson-Schwinger equations (jl()9)) . 
Of course, the solution of these equations requires T 2 , i.e., the complete set of 2PI diagrams. For 
practical purposes, a solution of the Dyson-Schwinger equations is therefore impossible. However, as 
already discussed in Sec. 14.4.31 one may truncate T 2 at some given order. Such a truncation defines a 
many-body approximation, within which a solution of the Dyson-Schwinger equations becomes feasible. 

The most simple, nontrivial truncation of T 2 is to include only the set of double-bubble diagrams 
shown on the left in Fig. 1201 The self-energies (|11(J|) computed from these diagrams consist only of the 
tadpole diagrams shown on the right-hand side of Fig. I2H This is the well-known Hartree approximation 
In general, all particle species in a particular effective theory contribute via tadpole diagrams 
to the self-energy of a given particle species. Since the tadpole diagrams do not have any dependence 
on the external momentum, the self-energies are (temperature-dependent) constants in the Hartree 
approximation. They simply shift the meson masses as compared to their vacuum values. In principle, 
the Dyson-Schwinger equations (jl(J9|) for the propagators are coupled integral equations, but in the 
Hartree approximation, they reduce to a system of coupled fix-point equations for the meson masses. 
While numerically much simpler than solving integral equations, the solution can still be a formidable 
task, if the underlying chiral symmetry group is large. For an 0(4) chiral symmetry, this problem was 
solved in Refs. |Jjij[EI2]- The U(3) r x U(3) e case was discussed in Ref. [HI]. The cases U{2) r x U(2) e 
as well as U(4) r x U (4)^ were investigated in Ref. |144j . (The U(l)y of baryon number is always trivially 
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Figure 21: Scalar and pseudoscalar meson masses as a function of temperature in the Hartree 
approximation, for effective theories models with different chiral symmetry: 0(4), U(2) r x 
U(2)i, U(3) r x U(3)g, and U(4) r x £7(4)^. Chiral partners are shown in the same panel to 
demonstrate that their masses become degenerate at high temperature. From Ref. 



respected in these models, cf. the discussion in Sec. 12.3.31 Nevertheless, in order to simplify the notation 
I include it in characterizing the symmetry of a particular chiral effective theory.) 

In Fig.|^the masses for the scalar and pseudoscalar mesons are shown as a function of temperature, 
calculated within the framework of linear sigma models with 0(4), U(2) r x U(2)g, U(3) r x U(3)e, and 
[7(4) r x U(A)i chiral symmetry. The different subpanels show the masses of the respective chiral partners. 
In the non-strange sector, these are (a) the a meson and the pion, and (b) the ao and the r\ meson. In 
the strange sector, these are (c) the k meson and the kaon, and (d) the fo and the rj' meson. Finally, 
in the charmed meson sector one has (e) the scalar D and D sQ meson and the pseudoscalar D and 
D s meson, and (f) the r] c and Xco meson. One observes that in all cases the chiral transition occurs at 
temperatures of the order of 200 — 300 MeV. For the results shown in Fig. |^ realistic, nonzero values 
for the quark masses were assumed, and consequently the transition is crossover, cf. the discussion in 
Sec. I2.3.4I Figure |^ also allows to compare the results for models with different chiral symmetry. 
For instance, Fig. |^ (a) nicely demonstrates the effect of enlarging the symmetry group. In general, 
the higher the symmetry, the more particle species contribute via tadpole diagrams to the mass of a 
particular species, and consequently the larger is its mass at a given temperature. Furthermore, one 
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Table 3: The critical temperature T c for the chiral transition, computed for various chiral 
effective theories in the chiral limit with and without U(1)a anomaly [144], in comparison 
to results from lattice QCD, extrapolated to the chiral limit. For Nf = 2 flavors, the QCD 
results for Wilson and Kogut-Susskind fermions from Tabled have been averaged, assuming 
uncorrelated statistical errors. 



Chiral model 


T c with U(1)a anomaly 


T c without t/(l) a anomaly 


LQCD 


0(4) 


(159.5 ±0.2) MeV 


N/A 


(172 ± 9) MeV 


U(2) r x U(2) t 


(154.5 ±0.2) MeV 


(149.2 ± 0.2) MeV 


(172 ± 9) MeV 


U(3) r x U(3) t 


(165.5 ±0.2) MeV 


(147.5 ± 0.2) MeV 


(154 ±8) MeV 



can learn from all subpanels that the difference between the U(3) r x U(3)e and the {7(4) r x U(4)e 
model is at best marginal. The reason is that the additional charmed meson degrees of freedom in the 
latter are so heavy that the contribution of the respective tadpole diagrams is exponentially suppressed 
~ exp(— m/T). Consequently, they only minimally influence the behavior of the non-charmed mesons 
in the temperature range considered here. In turn, in this temperature range the charmed meson masses 
change little from their vacuum values, cf. Fig. |^ (e), (f). 

Lattice QCD calculation cannot be done in the chiral limit, because as the quark mass decreases, 
the computation of the fermion determinant becomes more and more difficult. The only way to make 
predictions about the chiral limit is to extrapolate data obtained for nonzero quark masses. This is 
computationally expensive, as it requires calculations at several different values of the quark mass. In 
contrast, in the framework of linear sigma models, taking the chiral limit actually simplifies the calcu- 
lation. A comparison between quantities computed in lattice QCD and extrapolated to the chiral limit 
with the corresponding values obtained in linear sigma models is therefore straightforward. Let us di- 
cuss two examples. The first is the critical temperature for the chiral transition. The way to identify the 
critical temperature is the following. For a first-order phase transition, look for the temperature where 
two minima of the effective potential become degenerate. For a second-order transition, the criterion 
is a vanishing second derivative of the potential in one direction in order parameter space (indicating 
a massless degree of freedom and critical behavior). It should be mentioned that the Hartree approx- 
imation does not always reproduce the correct order of the transition as predicted from universality 
arguments. (Aside from that, being a mean-field type approximation, it always fails to predict the 
correct critical exponents.) For instance, for the 0(4) model and the U(2) r x U(2)i model with U(1)a 
anomaly, the transition should be of second order, but the Hartree approximation predicts a first-order 
transition. For the U(2) r x U(2)i model without U(1)a anomaly, as well as for the U(Nf) r x U(Nf)i 
models with Nf > 3, the Hartree approximation predicts a first-order transition in agreement with 
universality arguments. In chiral effective theories, one can simply "switch off" the U(1)a anomaly by 
setting the constant c = in Eq. (|14|). This is not possible for lattice QCD: the U(1)a anomaly may or 
may not be present, depending on how strongly instantons are screened at T c . This question has been 
studied on the lattice |145j . with the result that the U(1)a anomaly is not completely absent at the 
critical temperature, but at least rapidly decreasing. 
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Figure 22: The quark-mass diagram computed in a U(3) r x U(3)e chiral model |146j . with 
or without U(1)a anomaly for different values of the a meson mass in vacuum. The physical 
point (m q , m s ) = (10, 150) MeV is shown as a diamond. 



The critical temperatures obtained in various chiral models, with and without the U(1)a anomaly, 
are shown in Table El in comparison to the results from lattice QCD, cf. Table |2j Two things are 
noteworthy. First, the values of the critical temperature obtained in chiral models are rather close 
to the ones computed in lattice QCD. This is surprising given the fact that chiral effective theories 
only contain a certain subset of the degrees of freedom of full QCD, do not describe confinement, and 
moreover are treated here in the framework of a very simple many-body approximation. Second, the 
ordering of the critical temperatures as a function of quark flavors is the same in lattice QCD and chiral 
models without U(1)a anomaly, but the opposite in chiral models with U(1)a anomaly. This also lends 
support to the above mentioned results |145j regarding the rapid decrease of the U(1)a anomaly near 
T 

- 1 c- 

The second example is the investigation of the quark-mass diagram. Figure 1221 shows this diagram 
as determined within a U(3) r x U(3)e linear sigma model |146j . Of course, the linear sigma model 
does not have quark degrees of freedom. Nonzero quark masses correspond to a term of the form (fTB^) 
in the Lagrangian, where the matrix H is proportional to the quark mass matrix. Since the vacuum 
expectation value of $ is a diagonal matrix, the matrix H must also be diagonal. If one assumes SU(2)y 
isospin symmetry, H = h T + h 8 T 8 . Consequently, m g = a(h + hs/V2), m s = b(h — \/2h$). The 
fields ho, h 8 can be determined from the vacuum values for the pion and kaon masses, as well as the 
pion and kaon decay constants. Then, setting m q = 10 MeV, m s = 150 MeV fixes the values for the 
constants of proportionality a and b; for details see Ref. |146j . 

One observes in Fig. E21 that the position of the line of second-order phase transitions between the 
first-order region around the origin and the crossover region depends sensitively on the value of the a 
meson mass in vacuum. However, the physical point is always in the crossover region, if the U(1)a 
anomaly is present. Without the U(1)a anomaly, the transition could be of first order, if the a meson 
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is sufficiently heavy. 

Recently, there have been attempts to go beyond the Hartree approximation by including diagrams 
with more complicated topologies in T 2 [1321 1147] . In this case, a self-consistent calculation becomes 
technically rather difficult, because these topologies lead to momentum-dependent self-energies, so that 
the Dyson-Schwinger equations turn from simple fix-point equations for the masses into integral equa- 
tions for the full propagators. Moreover, momentum-dependent self-energies have in general nonzero 
imaginary parts, such that quasiparticles, which are stable in the Hartree approximation, develop a 
finite decay width. The spectral density carries the complete information of the spectral properties of 
the quasiparticles in the system. It is therefore natural to solve the Dyson-Schwinger equations for the 
full propagators as self-consistency equations for the respective spectral densities. Work in this direction 
is in progress |148j . 

5 Color Superconductivity 

5.1 Derivation of the Gap Equation 

As discussed in Sec. I2.4.3[ quark matter at high density and sufficiently low temperature is most probably 
a color superconductor (HEUIS]- While the color quantum numbers of a Cooper pair are determined by 
the attractive gluon interaction in the color-antitriplet channel, there are still many different ways to 
combine flavor and spin quantum numbers, giving rise to a plethora of possible color-superconducting 
phases. In general, the energetically most favorable phase will prevail at a given T and u. In order 
to decide which is the most favorable phase one has to compute the gain in condensation energy when 
forming Cooper pairs. To this end, one has to calculate the color-superconducting gap parameter. In 
general, the larger the gap, and the more degrees of freedom participate in forming Cooper pairs, the 
larger the gain condensation energy. The gap is computed from a so-called gap equation. In this section, 
I shall outline the derivation of this equation. 
Consider the QCD action 



is the action for non-interacting fermion fields in the presence of a chemical potential. For reasons which 
will become clear in the following, Sp and the quark-gluon coupling are rewritten in terms of ordinary 
quark fields, ijj, ip, and their charge- conjugate counterparts, tpc = ip T C, 4>c = Cip T , where C = £7 2 7o 
is the charge conjugation matrix. To this end, introduce the so-called Nambu-Gor'kov basis, with the 
2 • 4 iV c iV/-dimensional quark spinors 





where Sa is the gauge field action (including ghost and gauge fixing contributions) and 




(122) 




(123) 



For the fermion action one then obtains 




(124) 
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where 

So 1 = ( [Go °]-i)' [G^ 1 (X,Y) = -i(i 1 .d x ±^ -m)S^(X-Y) , (125) 

is the free inverse fermion propagator in the Nambu-Gor'kov basis. The quark-gluon coupling becomes 

i;(X)YT a ^(X)Al(X) = l^(X)T^(X)A;(X) , (126) 



where 



ft* 7 7" _ 7 ° T , (127) 



is the Nambu-Gor'kov matrix which couples the corresponding quark spinors to gluon fields. 
Now add a bilocal source term to the QCD action, 



S[)C]=S+[ V(X)JC(X,Y)V(Y) , (128) 
Jxy 



where, in the Nambu-Gor'kov basis, 



K= 1° . (!29) 



Here, a + and a~ are sources which couple to adjoint quark spinors and quark spinors, while ip + cou- 
ples to two quark spinors (remember that ipc ~ 4> T ) an d f~ couples to two adjoint quark spinors, 
respectively (ipc ~ r )■ Charge-conjugation invariance of the action relates the sources a + and o~~: 
(t~ = C [o~ + ] T C^ 1 . The action must also be real-valued, which leads to the condition ip~ = 7o[v 3+ ]^7o- 

The next step is to derive the effective action from the partition function of QCD in the presence of 
the external source /C, 

Z[K\ = J VA°VVVV expS[/C] . (130) 

The details of this derivation are beyond the scope of this review, but the interested reader can readily 
convince himself that, in the presence of bilocal sources, this problem is solved precisely by the CJT 
formalism |127j discussed in Sec. 14.4.31 Consequently, if one takes into account that there is no nonvan- 
ishing expectation value for a fermionic one-point function, and if one adapts Eq. ()106|) to the notation 
of this section, the effective action reads 

T[A,A,S] = I[A]--Ti lnA^-^Tr (A^A-lW-Tr mS'^-Tr fa" 1 S - l) +r 2 [A, S] , (131) 



2 2 V u J 2 2 

where A is the free and A the full gluon propagator, and S is the full fermion propagator in the Nambu- 
Gor'kov basis. The factors of 1/2 in front of the fermionic terms account for the double- counting of 
fermionic degrees of freedom in the Nambu-Gor'kov formalism. The full propagators for the physical 
situation (where /C = 0) are obtained from the stationarity conditions (|108j) . These conditions are 
Dyson-Schwinger equations of the type ()109|h 

A" 1 = A 1 + IT , S- 1 = + S . (132) 

(In a slight abuse of notation, we also denote the full gluon propagator at the stationary point of 
T[A, A, S] by A and, similarly, the full quark propagator by S.) 

In order to proceed one has to make an approximation for IV As in Sec. 14.4. 3^ the discussion will be 
restricted to the set of diagrams shown in Fig. El however, with the quark propagators in the quark loop 



66 




Figure 23: The quark self-energy. Writing out the individual components in the Nambu- 
Gor'kov basis one obtains the diagrams on the right-hand side. 



now given by the Nambu-Gor'kov propagator S. The gluon self-energy is computed as II = —2 ^I^/^A, 
i.e., by cutting a gluon line in the diagrams of Fig. El Thus, II is given by the diagrams shown in Fig.lTol 
(with Nambu-Gor'kov propagators S in the quark loop). The quark self-energy in the Nambu-Gor'kov 
basis is 

H^H^- < i33 » 

It is given by the diagram shown on the left-hand side in Fig. 1221 The diagonal components 
correspond to the ordinary self-energies for particles and charge-conjugate particles. In space-time, 
S + (X, Y) (£~(X, Y)) has a particle (charge-conjugate particle) entering at X and another particle 
(charge-conjugate particle) emerging at Y . On the other hand, the off-diagonal components $ ± have 
to interpreted as follows: a particle (charge-conjugate particle) enters $ + (X, Y) (<£>~(X, Y)) at X and 
a charge- conjugate particle (an ordinary particle) emerges at Y. This is typical for systems with a 
fermion-fermion condensate in the ground state |lf lj . The self-energies $ ± symbolize this condensate. 
Note that the two off-diagonal components are related in the same way as the bilocal sources ip in 
Eq. (j!29|) . i.e., $~ = 7o[$ + ]^7o- I n the following, also the term gap matrix will be used for $ + . 
In momentum space, the self-energy is given by 

E(JO = ~9 2 J K S(Q) VI A*(K - Q) . (134) 

This equation can be decomposed in terms of its Nambu-Gor'kov components. To this end, let us 
first determine the full Nambu-Gor'kov quark propagator by formally inverting the Dyson- Schwinger 
equation (|132j) for the quark propagator [149J, 




S=[Z+ n- > ( 135 ) 



where 

G ± = |[G^]- 1 + S ± -^ ([G^-' + ^y 1 $±| 1 , (136) 

[G^]- 1 + Y^y 1 $ ± G ± . (137) 

Here G + (G~) is the propagator for quasiparticles (charge-conjugate quasiparticles). Besides these quan- 
tities describing the ordinary propagation of quasiparticles, there are also off-diagonal, or "anomalous" 
propagators H ± in Eq. (|135j) . These anomalous propagators are typical for superconducting systems 
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and account for the possibility that in the presence of a Cooper-pair condensate, symbolized by 
$ ± , a fermion can always be absorbed in the condensate, while its charge-conjugate counterpart is 
emitted from the condensate and continues to propagate. 

In terms of its Nambu-Gor'kov components, the quark self-energy is 

E+(J0 = -g 2 J Q rT a G + (Q)YT b Af u (K-Q) , (138) 

= -g 2 J Q rT^G-(Q)YT^Af 1/ (K-Q) , (139) 

S+(JO = g 2 J Q YT^E + (Q)YT b Af u (K -Q) } (140) 

dT(if) = g 2 J Q ^T a E~(Q)YT^Af I/ (K~Q). (141) 

The integrals in these equations are shown diagrammatically on the right-hand side of Fig. ESI In these 
diagrams, a normal full propagator G + (G~) is denoted as a double line with an arrow pointing to the 
left (right). According to Eq. (|137|h the anomalous propagators S + and S~ consist of a combination of 

propagators {[Gq ] _1 + E 1 ^ , gap matrices $ ± , and full propagators G ± . This combination is explicitly 

drawn on the right-hand side of Fig. I2S1 with propagators ([Go] -1 + as single lines with arrows 

pointing left /right, and gap matrices $ ± as full/empty circles. Inserting these self-energies (and the 
corresponding one for the gluons) into the Dyson- Schwinger equations ()132j) one obtains a coupled 
set of integral equations which has to be solved self-consistently. In particular, the Dyson- Schwinger 
equations for the off-diagonal components ^ of the inverse propagator S^ 1 , i.e., Eqs. (I14()j) and (jl41|) . 
are the gap equations for the color-superconducting condensate. 

A completely self-consistent solution of the Dyson- Schwinger equations is technically too difficult 
to be feasible, and one has to make certain approximations. It turns out that in weak coupling a 
well-controlled, approximate solution is possible. This will be demonstrated in Sec. 15.31 Prior to that, 
I shall discuss the excitation spectrum in a superconductor, which follows from the poles of the quark 
propagator ([135)1 . 



5.2 Excitation spectrum 

In order to find the quasiparticle excitations, one has to determine the poles of G ± and 5*. For an 
arbitrary quark mass m, this is a formidable task, see Ref. |15Uj . For our purpose it is sufficient to 
consider the ultrarelativistic limit, m = 0, where the situation simplifies considerably |151j . Let us also 
focus exclusively on G + ; the poles of G~ and 5 ± can be determined accordingly (in fact, they have 
precisely the same poles) . In order to proceed, we need to specify the color, flavor, and Dirac structure 
of the inverse free propagators [Gq and of the self-energies and $ ± . The inverse propagators 
[Gq] -1 are diagonal in color and flavor space. To determine the Dirac structure, it is convenient to 
Fourier transform them into energy-momentum space and then to expand them in terms of projectors 
onto states of positive or negative energies, A£ = (1 + e^oj ■ k)/2, e = ±, 

[G±]-\K) = 7 • K ± a l0 = 7o £ [*o - (ek T /i)] K • (142) 

For the weak-coupling solution discussed in Sec. 15.31 it turns out that it suffices to compute the self- 
energies E ± neglecting effects from the breaking of the [>S£7(3) C ] color symmetry due to condensation 
of Cooper pairs. Then, the self-energies are diagonal in color and flavor space, and one only needs to 



know their Dirac structure. One may expand them in terms of Dirac matrices, 



^ ± -47° + s ± -7 + ... . (143) 

Other Dirac structures are simply abbreviated by dots. It turns out that, in weak coupling, the dominant 
contribution to the quark self-energy arises from the exchange of almost static (i.e., Landau-damped) 
magnetic gluons and, in momentum space, is 



( M 2 \ 
Z + (K) = Z~(K)~ lo f jh \n — + i-K\k \ J , (144) 

where g = g/(3V2n) and M 2 = (3-n/A)m 2 g . The imaginary part gives rise to a finite lifetime of 
quasiparticle excitations off the Fermi surface |149j , but in weak coupling this contributes to the color- 
superconducting gap only at sub-subleading order. (How to count orders in weak coupling will be 
discussed in detail in Sec. 15.3)) . There are other contributions at sub-subleading order which have not 
been reliably computed so far. Only the subleading order terms are under complete analytic control. 
They determine the prefactor of the color-superconducting gap to order 0(1). The sub-subleading 
contributions contribute to the prefactor at order 0(g). Therefore, in the following the imaginary 
part will be neglected and only the real part will be kept, which, for ln(M/A;o) ~ 1/g contributes at 
subleading order to the color-superconducting gap |153| . Defining the wave function renormalization 
factor 



M 



2\ - 1 



Z(k )= ll+^ln-^-l , (145) 

the effect of the quark self-energy is to shift the poles in the propagator, k Q — > k /Z(k ). Note that the 
logarithm renders a normal-conducting system a non- Fermi liquid |154j . 

The color, flavor, and Dirac structure of the off-diagonal self-energies $ ± is less trivial. As they 
symbolize the Cooper-pair condensate, they must have the quantum numbers of the particular channel 
where condensation occurs. For the purpose of illustration, I specify this structure for the four color- 
superconducting phases already discussed in Sec. 12.4.31 Let us furthermore only consider parity-even 
spin-zero and spin-one condensates. In this case, one may expand the gap matrix $ + (in momentum 
space) as follows [ST] . 

$+(iv) = £0 e (K).M k A^, (146) 

e=± 

where 6 , the so-called gap function, is a scalar function of 4-momentum and A^k is a matrix in color, 
flavor, and Dirac space, which is determined by the symmetries of the color-superconducting order 
parameter (cf. Sec. 12.4. 3J) . An important property is that it commutes with the energy projectors, 
M k .\ k 0. 

In Table |U the explicit expressions for A^k are listed for the four phases considered here. For these 
expressions, I have used the fact that a color antitriplet is totally antisymmetric and thus has a repre- 
sentation in terms of the antisymmetric Gell-Mann matrices A2, A5, and A7. These Gell-Mann matrices 
form an SO(3) subgroup of SU(3), and are thus identical to the generators of 5*0(3), (Ji, J 2 , J3) = J. 
These matrices were finally used in Tabled to parametrize that the gap matrix $ + is a color- antitriplet. 
Similarly, r 2 (the only Pauli matrix which is antisymmetric) symbolizes that the gap in the 2SC phase is 
a flavor singlet. The Dirac matrix 75 in the 2SC and CFL cases is necessary to obtain even parity. For 
the CFL case, the flavor-antitriplet nature is represented by another set of generators (I 1 , 1 2 , 1 3 ) = I of 
SO (3). Color-flavor locking is obtained by a scalar product of I with J, cf. also Tabled For the two 
spin-one phases, the CSL phase and the polar phase, the order parameter is a 3-vector, see Sec. I2.4.3| 
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Table 4: The structure of the matrices .M k and L k in various color-superconducting phases, 
(r 2 ) f9 = ie f9 , (J k ) lj = -ie ijk , (I h ) f9 = -ie fah , and 7 ± (k) = 7 - 7 • kk. For the matrix L k 
in the CSL phase, the second term is a matrix in color space, formed by the dyadic product 
of the two vectors k + 7j_(k) and k — 7 ± (k), and a matrix in Dirac space, formed by the 
product of the two 7j_(k) matrices. The last two columns show the two eigenvalues A r of 
L k and their degeneracy d r (counting color-flavor degrees of freedom in the 2SC and CFL 
phases, and color-Dirac degrees of freedom in the two spin-one phases). 



phase 






Ai (di) 


A 2 (d 2 ) 


2SC 


75 T 2 J3 


(Js) 2 (r 2 ) 2 


1(4) 


0(2) 


CFL 


75l J 


(I.J) 2 


4(8) 


1(1) 


CSL 


J - [k + 7 _L(k)] 


l + [k + 7± (k)] [t-7i.(k)] 


4(4) 


1(8) 


polar 


J3 [jfes + T^k)] 


Us) 2 


1(8) 


0(4) 



but the gap matrix is a scalar. This requires it to be proportional to a scalar product of the order 
parameter with another 3- vector. There are only two other 3- vectors (in momentum space), the direc- 
tion of momentum of the quark in the Cooper pair, k, and the vector 7. Consequently, $ + has to be 
proportional to a linear combination of these two 3-vectors. It is convenient jlU] to use the projection 
of 7 onto the subspace orthogonal to k, 7 ± (k), because then M.^ commutes with the energy projectors. 
Finally, color-spin locking requires a scalar product of J with this linear combination of 3-vectors. In 
the polar phase, one may independently choose a direction for the gap in color space and in space-time. 
Conveniently, one chooses the 3- (anti-blue) direction in color space and the ^-direction in space-time. 

In order to proceed, however, one does not require the explicit form of .Mk hi the various phases. The 
existence of the decomposition ()146|) and the commutation property of M.^ with the energy projectors 
is sufficient to derive the quasiparticle spectrum as a function of the absolute magnitude of the gap 
function, \<f) e (K)\. To see this, compute $~ = 7o[$ + ]^7o an d, together with $ + and [Gq]^ 1 + £~, the 
quantity 

$- ([GoT 1 + 2-)^ $ + ([Gor + ir) = £ |0WL k A k e , (147) 

e=± 

where 

£k = 7o-A4-Mk7o (148) 

is another central quantity for the quasiparticle excitation spectrum. Expressions for the matrix Lk in 
the various phases are also listed in Table 0] Note that also L k commutes with the energy projectors, 
[L k ,A k ] = 0. Since L k is hermitian, it has real eigenvalues, A r , and can be expanded in terms of a 
complete set of orthogonal projectors, V^, 

Lk = £A r K. (149) 
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In the four phases considered here, there are only two distinct eigenvalues and therefore two distinct 
projectors. The eigenvalues are also listed in Table HJ and the projectors can be expressed in terms of 
L k via 

Lk_A2a (150) 



1,2 



Ai, 



A 



2,1 



Obviously, also these projectors commute with the energy projectors, [P k ' , Ajj] = 0. It is now straight 
forward to compute the full propagator by inverting the term in curly brackets in Eq. (H 
projectors V r A£ form a complete, orthogonal set in color, flavor, and Dirac space, 



since the 



G+(K) 



1 



[h/Z(k )] 2 - 



where 



-k,r v 



) = {ek - fi) 2 + \ r 



ie |2 



1/2 



(151) 



(152) 



Obviously, the poles of the full propagator are located at k = ±Z(k ) el_ r ((j) e ). Because of the k 
dependence of the gap function <j) e (K), this is a condition which has to be solved self-consistently. 

In order to get an impression what the excitation spectrum of the quasiparticles in a superconductor 
looks like, let us for the moment approximate <f) e (K) = <p — const., and also set Z(ko) = 1 (corrections 
are of order 0(g)). Let us also neglect the fact that there are two different sets of excitation branches 
depending on the value of A r . In Fig. |^ the excitation spectrum is shown for non-interacting massless 
fermions as well as for quasiparticles in a superconductor. The dispersion branches for the quasiparticle 
excitations corresponding to negative energies, e = — , i.e., the quasi-antiparticles and quasi-antiparticle- 
holes, hardly differs from the non-interacting antiparticle or antiparticle-hole branches. As we shall see 
in Sec. 15.31 in weak coupling, <C /i, such that to very good approximation e kr ~ k + \i. On the other 
hand, the dispersion branches for the quasiparticle excitations corresponding to positive energies, e = +, 
differ considerably from the non-interacting particle or hole branches. The most notable feature is an 
energy gap at the Fermi surface, k = fi, between the quasiparticle and quasiparticle-hole branches. This 
indicates that, in a superconductor, it costs an energy 2 to excite quasiparticle-quasiparticle-hole pairs 
at the Fermi surface. In contrast, in a non- interacting system it costs no energy to excite particle- hole 
pairs at the Fermi surface. The superconducting state is thus energetically favored compared to the 
normal-conducting (non-interacting) state. As a rule of thumb, the more fermionic excitations branches 
develop a gap (the more fermions form Cooper pairs), and the larger the associated gap, the lower is 
the ground state energy, and the more energetically favored is the particular superconducting state. 

In the CFL and CSL case, there are two different excitation branches with two different gaps, 
e M,i(^) = V^i^ = 20 and e^ j2 (0) = y/^2<fi = 0- Consequently, it costs twice the amount of energy 
to excite quasiparticle excitations from the first branch than from the second. In the 2SC and polar 
phases, there also two different excitation branches, but the one corresponding to unpaired fermionic 
excitations is gapless, et 2 — 0- In the next section, the gap equation ()140|) for the gap function 4> e (K) 
will be solved, which allows us to determine the magnitude of the gap at the Fermi surface. 



5.3 Solution of the Gap Equation 

In order to solve the gap equation (jl4()j) . insert Eqs. ()146j) and ()151j) into Eq. ()137|) for 



2+(AO = -£ 7 oMc7o7W 



[k /Z(ko)] 2 - e£ r (0 



(153) 
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Figure 24: The excitation spectrum for (a) non-interacting massless particles and (b) quasi- 
particles in a superconductor. The value of the gap function <fr e (K) is assumed to be con- 
stant, (f) = 0.5 jji,. The excitation energies for non-interacting particles are e k (0) — k — fi 
for particles (solid), e k (0) = // — k for holes (dashed), e k (0) = — (A; + /i) for antiparti- 
cles (dotted), and e k (0) = fc + /J for antiparticle-holes (dash-dotted). The excitation en- 
ergies for quasiparticles are e k (0) 

quasiparticle-holes (dashed), e k (0) 



\J{k — /i) 2 + 2 (solid), e k (0) = \J{k — /i) 2 + 2 for 
(> 2 for quasi-antiparticles (dotted), and 



/(A; + /i) 2 

\J {k + /i) 2 + 2 for quasi-antiparticle-holes (dash-dotted) 115 lj . 



(By the way, this result demonstrates the claim made earlier that G + and H + have the same poles.) 
Now insert this equation into the gap equation (jl4U|) . multiply from the right with .M k A k , and trace 
over color, flavor, and Dirac space. The result is an equation for the scalar gap function (p e (K), 



<P e '(Q) 



W Z [q /Z( qo )] 2 - 



-q,sv 



Q)UT (k,q) 



where 



ab,s 



(k,q) 



Tr 






Tr 







(154) 



(155) 



The form of this gap equation is the same for all color-superconducting phases considered here. The 
difference lies in the excitation spectrum and the structure of the term T^f e (k, q). 

At this point it is convenient to explain the power counting scheme in weak coupling, g <C 1. The 
right-hand side of Eq. ()154j) has a prefactor g 2 . Consequently, in order to satisfy the equality, after 
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performing the integral there have to be terms ~ <p/g 2 , which together with the prefactor produce a 
term ~ 0(0), i.e., which is of the same order as the left-hand side. These are the so-called terms of 
leading order in the gap equation. Then there are terms of so-called subleading order. These enter the 
right-hand side of the gap equation at order 0{g<p). The terms of so-called sub- subleading order are 
~ 0(g 2 (p). It turns out that only the terms of leading and subleading order can be reliably calculated 
in weak coupling. 

In order to proceed, one has to make further approximations. As shown in Ref. |155j . to leading 
and subleading order one does not need the fully self-consistent gluon propagator; it suffices to employ 
the gluon propagator in HDL approximation, cf. Sec. 14.4.11 The HDL gluon propagator is diagonal in 
adjoint colors, A°f v = 5 ab A^. For the sake of definiteness, I choose pure Coulomb gauge, where 

A 00 (P) = Ae(P) , A 0t (P) = , A^P) = A t (P) (<% - ptpj) , (156) 

with P = K — Q and the longitudinal and transverse propagators A^ introduced in Eq. (|93p. In 
fact, it is not even necessary to use the full form of the HDL propagator. In weak coupling, power 
counting along the lines of argument given above reveals |4lH H.56] that the dominant, leading-order 
contribution to the gap equation comes from almost static, Landau-damped magnetic gluons. Their 
propagator may be approximated by 

A^ DM (^) ^ 6 / M4 2 6(M - p) , (157) 
p b + ilf 4 CiT 

where M 2 = (3ir/4:)m 2 . To subleading order, there are two contributions, from non-static magnetic 
gluons, 

~±Q(p-M) , (158) 

pi 



and from static electric gluons, 



p l + 6m z g 



All other terms in the HDL gluon propagator contribute to sub-subleading order. 

In pure Coulomb gauge, one only needs the OO-component, T^f 6 (k, q), and the transverse projection 

of the ij-components of the tensor T^ u f e (k, q), 

T ^'(k, q) = - - p iPj ) VifiK q) , (160) 

where the extra minus sign is convention. In Ref. [21] it was shown that one can write these components 
in terms of a power series in p 2 / (kq), where p = |k — q|, with coefficients r^m which depend on k, q and 
the product ee' , and an overall normalization factor a s which is the same for the 00- and the transverse 
component. The power series over m start at m = — 1 and, for spin-zero and spin-one gaps, terminate 
at m = 2. The normalization factors satisfy the constraint J2 S a s — 1- 
Performing the Matsubara sum in Eq. ()154j) one then obtains 

a 2 r"*..^. -v^"'(^ tanh 



X 



{-2A^(p)4 m + [2Af M (p) + A^ M (ei s + el r ,p)+A 




Several approximations have been made to obtain this result. First, the integration over q has been 
restricted to a narrow interval of length 2 S around the Fermi surface, 5 ~ M. It turns out that this 
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approximation is good to subleading order; due to the momentum dependence of the gap function, the 
value of the cut-off affects the color-superconducting gap parameter only at sub-subleading order [40J. 
Second, when evaluating the Matsubara sum via contour integration, to subleading order the value of 
the poles q = ±Z(q ) e q s may be approximated by g ~ ±e qiS everywhere except in the residue of the 
contour integral. The single factor of Z(e qs ) under the integral arises from the residue. In the argument 
of this factor, one has also made the approximation go — =t e q S > since the logarithm in go m Eq. (jl45j) 
gives at most a subleading contribution to the integral [HU I15H] . Third, the gap function was assumed 
to be an even function of energy, (p e (k ) = (p e (—k ) J^UJ . Fourth, the gap function was assumed to be 
isotropic in momentum space, e (k) = <p e (k) [ST]. Finally, the imaginary part of <p e was neglected |40j. 

How do the leading, ~ 0(0), and subleading, ~ O(gtp), terms arise? To this end, one has to power- 
count the different contributions to the integral in Eq. (|161|) . One uses the fact (which will be confirmed 
below) that in weak coupling, cf) ~ jj, exp(— 1/g). Neglecting all subtleties regarding different excitation 
branches, and setting e q ' s = e q = (e'q — /i) 2 + 2 , for e' = + the integral over q gives rise to a term 

r S *JL = 2 f _i!£_ = 2 in ( H^HZ) K 2 ln „ ln \—JL— 1 „ i , (162) 

where I have substituted the variable £ = q — \i and used the fact that <5~M~gyu;^>0~ //exp(— 1/g). 
The logarithm appearing in the estimate ()162|) is called the "BCS-logarithm", because it also appears 
in standard BCS theory For e' = — , the BCS-logarithm does not occur, as e q ~ q + fi, such that 

the integral is parametrically only of order 8/ fi ~ Mj // ~ g (provided that the gap function <p~ falls off 
sufficiently fast that one may restrict the integral to a narrow range around the Fermi surface). 

For the leading-order contribution to the gap equation, we need another term which is also ~ 1/g, 
such that this term and the BCS-logarithm combine to give a contribution ~ 1/g 2 which cancels the 
prefactor g 2 in front of the integral in Eq. (jl61j) . To estimate the order of magnitude of the remaining 
terms, one notes that the coefficients are parametrically at most of order 0(1) [Slj, such that they 
can be neglected for the purpose of power counting. The term ~ 1/g which we are looking for arises 
from the term m = in the sum over m in conjunction with the Landau-damped magnetic gluon 
propagator. Abbreviating u± = e qs ± e^ r , one estimates 

/ dppA):™(uj ± ,p) = / dp P 2 = - In 

Jk-q Jk~q p b + M 4 U)± 6 

where the approximation k ~ q was used. (Only when k — q ~ 0, the logarithm may become large, see 
argument below. Otherwise, if k — q ~ M, the logarithm is only of order 0(1), not 0(1/ g).) If either 
e = — or e' = — , or both e = e' = — , the logarithm is parametrically of order 0(1), and not ~ 1/g. 
Consequently, the only case of interest is if both e — e' — +. In this case, u± ~ 0, and the logarithm is 
large, ln(/x/0) ~ 1/g. 

One readily convinces oneself that the p integral over the other terms in Eq. (J161)) gives at most 
a contribution of order 0(1). In combination with the BCS-logarithm, this leads to a subleading 
contribution in the gap equation. From the quasi- antiparticle poles, e' — — , one does not obtain a large 
BCS-logarithm, but a term ~ g. With the prefactor and a factor of <fi from the gap function under 
the integral, their contribution to the gap equation is of order g 3 <p, i.e., even beyond sub-subleading 
order. In the following, one may therefore safely neglect the contribution from quasi-antiparticles when 
computing the gap for the quasiparticle excitations, e = +. 

In order to proceed, one performs the p integrals, which can be done exactly [STJ, but only need 
to be known to order 0(1). One furthermore analyzes the coefficients ^ 2 m an d realizes jHI] that, to 



M 6 + MW ± 
(k - g) 6 + M 4 c4 



in U ' (163) 
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Table 5: The normalization factors a s , the coefficients an d the constant d from Eq. ()165|) 
in various color-superconducting phases. In the polar phase, is the angle between the 
direction of the color-superconducting order parameter and the momentum of the quarks in 
the Cooper pair. 
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subleading order, one may approximate k ~ q ~ fx in the expressions for these coefficients. Thus, 
they become pure numbers of order 0(1). Moreover, to subleading order the coefficients rjj^ = and 
need not be considered further. I list the coefficients for m = 0,1,2 in Table El together with the 
normalization factors a s for the four color-superconducting phases considered here. The final result for 
the gap equation can be written in the concise form ^T] (let us omit the superscript "+" for the sake 
of simplicity) 

0(e k , n k) = f \l d(g - n) E a s Z(e qjS ) tanh (^) \ In ( — ^ ) , (164) 



which is exact to subleading order. In Eq. ()164j) we have introduced 

/ 2 \ 5/2 6 r i 

b = bex V (-d), 6 = 256vr 4 -— , d = -- 4 + rf 2 + 2(rji + rjl) . (165) 

\N f g 2 J rf Q L 

The Nf- dependence of b arises from the corresponding dependence of the gluon mass parameter m g , 
cf. Eq. ()99|) . The values for the constant d are also listed in Table For the spin- zero color- 
superconducting phases, d = 0, due to an accidental cancellation of the coefficients rf 2 and rf 2 . This does 
not happen in the spin-one phases and, consequently, d ^ 0. 

In order to solve Eq. (J164)) . one makes the following approximation which was first proposed by Son 
and is valid to subleading order, 

* ©(e q , s - e Kr ) In (J^-) + 6(e k , r - e q , s ) In (^j . (166) 

The remainder of the calculation is technical, but straightforward and given in detail in Ref. [51J. To 
summarize the steps, a suitable substitution of variables allows to rewrite the gap equation ()164|) . which 
is an integral equation, in terms of Airy's differential equation jSU I153j . The result for the gap function 
has the form 

<f>(x r ) =(j)oF(x r ) , (167) 
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where 0o is the value of the gap function at the Fermi surface, i.e., the color-superconducting gap 
parameter or "gap", and F(x r ) parametrizes the momentum dependence of the gap function. The 
variable x r is defined as 

Xr = g\n(- ^ ) . (168) 

\k - fJL + e^J 

At the Fermi surface, k = fi, one has 

ln[26/x/(\/A^0o)] — 7r /2 + 0(g) ~ 0(1). If one 
moves away from the Fermi surface, x r stays of order 0(1), as long as the momentum difference from 
the Fermi surface is \k — fi\ ~ O (</>). When \k — fi\ ~ M or larger, x r ~ 0(g). The precise form of 
the function F(x r ) is not very illuminating (it consists of a combination of Airy functions PxT] IT53]). 
and thus will not be discussed here. All one needs to know is that it has a narrow peak in an interval 
\k — fi\ ~ 0(4>) around the Fermi surface. At the Fermi surface, x r = x*, the function F(x r ) assumes 
the value F(x*) = 1 + 0(g 2 ). At a distance \k — fi\ ~ M from the Fermi surface, F(x r ) ~ 0(g). If one 
neglects the factor Z(e^ s ) in Eq. ()164j) . the differential equation satisfied by the gap function is that of 
the harmonic oscillator and, consequently, the solution of the gap equation becomes simpler and more 
amenable to interpretation: F(x r ) = sinx r [4U| I156j . 
The value of the gap function at the Fermi surface is 

0o = 2 bb' of i exp (-IjA (Af \ a 2 T 1/2 • (169) 

The constant b' Q = exp[— (it 2 + 4)/8] arises from the wave function renormalization factor Z(e^ s ) in 
Eq. ()164J) 5^1 I153| . The result (j!69J) differs from the standard BCS result in the power of the coupling 
constant g in the exponent. In weak-coupling BCS theory, O ~ exp(— l/g 2 ), while here 4>q ~ exp(— 1/g). 
The difference in the parametric dependence on g arises from the long-range nature of magnetic gluon 
exchange. In BCS theory, the attractive interaction is assumed to be short-range (point-like or at least 
exponentially screened). On the other hand, in QCD static magnetic gluon exchange is not screened 
[T3*j . Almost static magnetic gluons are dynamically screened, but the screening is rather weak. It gives 
rise to the large logarithm (|163|) in addition to the BCS logarithm (|162j) . This reduces the power of g 
in the exponent. 

In Table E] I list the value of the gap 0o m units of its value in the 2SC phase. For the spin-one gaps, 
the nonzero value of the constant d leads to a strong suppression ~ e~ d ~ 10~ 2 — 10~ 3 as compared 
to the spin-zero gaps. In the CFL and CSL phases, the second gapped excitation leads to a nontrivial 
factor (A^ 1 A2 2 ) < 1, which reduces the gap as compared to the 2SC and polar phases where there 
is only a single gapped excitation. 

The result ()169j) is rigorously valid in weak coupling, i.e., for asymptotically large quark chemical 
potentials, where the value of the strong coupling constant evaluated at the scale \i is small, g(ji) <^ 1. 
However, for phenomenology it is of considerable interest to determine the gap also at values of \i 
which might occur in nature, for instance in the core of compact stellar objects. To this end, one 
extrapolates the weak-coupling result ()169|) to large values of g(fi) ~ 0(1). Such an extrapolation has 
to be considered under the caveat that the sub-subleading terms are not really small for g(fi) ~ 1 
and could lead to large deviations of the actual value of the gap from the subleading result ([169)1 . 
Nevertheless, the computation of <po at g(fi) <C 1 is a well-posed problem with a definite result, and so 
is its extrapolation to large values of g(fj). In this sense, this approach should be considered to be more 
reliable than ad hoc calculations within NJL-type models which are very popular in the description of 
color-superconducting quark matter [THj . 

For the running of the strong coupling constant g with /x/A, where A is the QCD scale parameter, 
I take the standard 3-loop formula (US]- I assume that there are only Nj = 3 active quark flavors 
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Table 6: The value of the gap function at the Fermi surface, 0o i n units of its value in the 
2SC phase, and the critical temperature, in units of its value expected from BCS theory, 
Eq. (|28|). and in units of the critical temperature in the 2SC phase, T C 2SC . 
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involved in the running of the coupling constant, so that in order to obtain the correct value of g(fi) at 
the mass of the Z boson, one has to adjust the QCD scale parameter, A = 364 MeV. I also take Nf = 3 
in the factor b in Eq. (J165|) . Physically, this means that, independent of the number of quark flavors 
which form Cooper pairs, there are always three (massless) quark flavors which screen color charges. 
The result of extrapolating Eq. (j!69|) to realistic values of fi is shown in Fig. (a) for the various 
color-superconducting phases considered here. 

One observes that the 2SC phase has the largest gap, 0q SC — 10 MeV followed by the CFL phase. 
The spin-one phases have gaps which are about 2 — 3 orders of magnitude smaller, 0q =1 ~ 1CT 2 — 1CT 1 
MeV. The gap is approximately zero for chemical potentials below 500 MeV, rapidly rises to assume a 
maximum around fi ~ 600 MeV and then decreases. (For larger values of [i it will eventually increase 
again.) This behavior is due to the dependence of 0o on g. For large values of g (small values of 
/i), the power-law behavior g~ 5 from the prefactor b leads to a suppression, while at small values of g 
(large values of /i), the exponential suppression ~ exp(— 1/g) dominates. This leads to a maximum for 
intermediate values of /i. (For asymptotically large values of /x, the gap increases again, because the 
prefactor fi dominates the /i dependence of the remaining factors.) 

One can also solve the gap equation at nonzero temperature. One finds that the shape F(x r ) of the 
gap function hardly changes with T, but that the value of the gap decreases [301 • The gap equation (I140|) 
is equivalent to the one obtained in the mean-field approximation |151j , and therefore the temperature 
dependence of the gap follows the predictions from mean-field theory. In particular, the transition to the 
normal-conducting phase is of second order, irrespective of the symmetries of the order parameter. The 
critical temperature T c for this transition can be computed analytically, for details see Refs. |4UlloT| [153j. 
The result is [HI] 

T c =-<MArA^) 1/2 [l + 0(s)] . (170) 

7T 

This result is surprising for two reasons. First, in a phase with a single gapped excitation, like the 

1 /2 

2SC phase or the polar phase, where (A" 1 A2 2 ) = 1, the critical temperature in QCD, measured in 
units of the gap, is the same as in BCS theory [40J, at least to leading order in weak coupling. This 
is unexpected, since we have seen that the dependence of the gap itself on g is parametrically very 
different than in BCS theory. Second, in a phase with two different nontrivial excitation branches, like 
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oc s : 3-loop, N f =3, A=364 MeV 
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Figure 25: (a) The gap and (b) the critical temperature as a function of the quark chemical 
potential. Solid curves are for the 2SC phase, and dotted curves for the CFL phase. (In 
the case of the critical temperature, both curves coincide). The long-dashed curve is for the 
CSL phase, and the dashed and dash-dotted curves are for the polar phase with $ = n/2 
and $ = 0, respectively. 



the CFL and CSL phase, the factor (A" 1 A^ 2 ) 1 ^ 2 ^ 1 violates the expectation from BCS theory [3U| 1ST]. 
In Table Hi show T c in units of the critical temperature expected from BCS theory, T C BCS = (e 7/ 7r)</>o, 
to demonstrate this violation. In physical units, say the value of the critical temperature in the 2SC 
phase, T C 2SC = (e 7 /vr)0g sc , the factor (A? 1 \ a 2 2 ) l/2 cancels against its inverse in Eq. (|169j) . This leads 
to the conclusion that, in the mean-field type approach pursued here, the critical temperatures in the 
2SC and CFL phases are actually identical. The critical temperatures in the spin-one phases are just a 
factor e~ d smaller than in the 2SC and CFL phases [59|. 

The critical temperature (jl70|) is shown as a function of [i in Fig. (b). These curves also define 
the boundaries of the color-superconducting phases in the phase diagram of nuclear matter, cf. Fig. ^ 
The subleading result ()170|) for T c implies that one would have to cool quark matter below temperatures 
of order 5 MeV, before one enters a color-superconducting quark matter phase (the 2SC or CFL phase). 
This means that, unless sub-subleading corrections to the gap (|169|) (and thus to T c ) are large, color- 
superconductivity is irrelevant in the context of heavy-ion physics, but that it may play a large role 
for compact stellar objects which have a sufficiently dense core. While spin-zero color-superconducting 
matter may occur already quite early in the evolution of such a compact stellar object, i.e., while it is 
still comparatively hot, matter in a spin-one color-superconducting state only occurs after the core of 
the stellar object has cooled below a temperature of order 10 keV, i.e., in the later stage of its evolution 
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5.4 Gluon and photon properties 

In this section, I take a first step towards a self-consistent solution of the Dyson-Schwinger equations 
(j!32j) and compute gluon properties in a color superconductor. Within the two-loop approximation 
to T2, the gluon self-energy consists of the diagrams shown in Fig. ED At temperatures of relevance 
for color superconductivity T < T c ~ O ~ H ex p( — 1/flO <S /i, we may neglect the contributions 
from the gluon (and ghost) loops to the gluon self-energy: they are ~ g 2 T 2 , while the quark loop is 
~ g 2 ^ 2 ^> g 2 T 2 . Thus, the gluon self-energy in momentum space is 

IW) = Y j K ^[KS{K)TlS{K-P)} , (171) 

where the trace runs over color, flavor, Dirac, and Nambu-Gor'kov space. By introducing the Nambu- 
Gor'kov basis one has effectively doubled the degrees of freedom by introducing charge-conjugate quarks 
in addition to quarks The factor 1/2 in Eq. (|17H) prevents overcounting these degrees of freedom. 

Similarly to the gluon self-energy one can compute the photon self-energy 11^ replacing the quark- 
gluon vertices in Eq. (j!71|) by the corresponding ones for the coupling between quarks and 

photons, 

£ / 7 M Q 
g{ ~-fQ 



where Q = diag(2/3, — 1/3, — 1/3) is the quark electric charge matrix. Furthermore, as discussed in 
Sec. 12.4. 2^ in a color superconductor gluons can mix with the photon, leading to a "rotated" electromag- 
netic [f/(l)] symmetry in the 2SC and CFL phases. This fact manifests itself in a nonvanishing "mixed" 
gluon-photon self-energy IT^, which follows from Eq. (j!71|) by replacing just one of the quark- gluon 
vertices with the quark-photon vertex ([172)1 . In order to determine the gluon and photon properties in 
a color superconductor, one has to compute all these different self-energies. For the sake of convenience, 
in the following let us set the index 7 = 9 and consider Eq. ()171|1 for a, b = 1, 2, . . . , 9. I also introduce 
T 9 = (e/g)Q as the appropriate generator for [£/(l) em ]. 

Taking the trace over Nambu-Gor'kov space in Eq. (|171|) . one realizes that the quark loop consists 
of four contributions, two "regular" ones with normal propagators G for quarks and charge-conjugate 
quarks and two with anomalous propagators 5 , 

2 

Kb(P) = J J K Tr[^T a G + (K)YT b G + (K - P)+^TjG-(K)YT?G-(K - P) 

-^T a E-(K) YT? ~ + (K - P) - YTj E + (K) ^T b E~(K - P)] . (173) 



The two different topologies corresponding to these contributions are shown in Fig. | 

To further evaluate the trace one has to specify which color-superconducting phase one would like 
to consider. The form of the propagators G ± , 5 ± can then be determined following the method outlined 
in Sec. 15.21 see Eqs. ()151j) . ()153j) . After inserting these propagators into Eq. ()173|) . one performs the 
Matsubara sum. The resulting expressions for the self-energies are rather unwieldy and will not be 
shown here. For the 2SC phase they were first derived in Ref. jlS] and for the CFL phase in Ref. 158J. 
For the spin-one color-superconducting phases, this was done in Ref. |159j . To further evaluate these 
expressions, one has to compute the integral over d 3 k. For an arbitrary gluon 4-momentum P M , this 
has not yet been done. However, in the static, homogeneous limit, the self-energy n^(0) was computed 
in the aforementioned references in order to derive the Debye and Meissner masses in the respective 
color-superconducting phases. The results will be discussed in more detail in the following. The gluon 
self-energy in the 2SC phase was also evaluated for nonzero energies and momenta po, p, which are small 
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(a) 



(b) 



Figure 26: Contributions from (a) normal and (b) anomalous quark propagation to the 
self-energies of gluons and photons, and the mixed gluon-photon self-energy. The notation 
follows that of Fig. [21 



compared to the quark chemical potential. This calculation is rather technical, and I simply refer to 
Ref. |121[ I160| for the details. The main result was that the modification of the gluon self-energy in a 
color superconductor does not influence the value for the gap parameter at leading or subleading order 
in weak coupling. For the other color-superconducting phases, a similar calculation has yet to be done. 

In general, in a medium at nonzero temperature and/or density static, long- wavelength (color-) 
electric fields are screened. The screening length is determined by the (inverse) Debye mass. If the 
medium is normal-conducting, static, long- wavelength (color-) magnetic fields are not screened. This 
changes in a superconductor, where the Meissner effect expels (color-) magnetic fields. They can only 
penetrate a certain distance into the superconducting medium. For static, long-wavelength (color-) 
magnetic fields, the (inverse) penetration length is determined by the so-called Meissner mass. The 
acquisition of a Meissner mass by a gauge boson indicates that the corresponding gauge symmetry is 
broken via the Anderson-Higgs mechanism. The Debye and Meissner masses are defined as 



I present the values for the Debye masses in Table and for the Meissner masses in Table |H1 for various 
color-superconducting phases. 

In the 2SC and polar phases, the Meissner mass of the first three gluons vanishes. These gluons 
correspond to the unbroken [SU(2)] C subgroup, cf. Tabled] What is interesting is that they also 
have a vanishing Debye mass, indicating that the corresponding color-electric fields are unscreened. 
Implications of this result were discussed in Ref. |161j . The other five gluons acquire both a Debye as 
well as a Meissner mass. Electric and magnetic fields are always screened in these color-superconducting 
phases. Only the eighth gluon mixes with the photon. Another interesting aspect is, however, that this 
mixing only occurs in the magnetic sector, electric and color-electric fields remain unmixed. In order 
to obtain the eigenmodes of the gauge bosons, one has to diagonalize the mass matrices for electric and 
magnetic gluons. A zero eigenvalue in this mass matrix indicates the presence of an unbroken "rotated" 
[[/(l)] gauge symmetry. 

In the CFL phase, all gluons acquire a Debye as well as a Meissner mass, indicating that the [SU (3)] c 
color symmetry is completely broken. Photons are Debye- as well as Meissner-screened, and there is 
again mixing between the eighth gluon and the photon. In contrast to the 2SC and polar phases, 
however, this mixing extends also to the electric sector. In the CSL phase, all gluons and the photon 
obtain Debye and Meissner masses. There is no mixing between the gluons and the photon. This means 
that the mass matrix of the (former gauge) bosons is already diagonal, and it has no zero eigenvalue. 
Consequently, there is no unbroken residual symmetry, and no room for a rotation that could generate 




(174) 
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Table 7: The Debye masses for gluons, photons and from the mixed gluon-photon polar- 
ization tensor for various color-superconducting phases. The results are given in units of 
N f /j, 2 /(Qn 2 ), where N f = 2 in the 2SC phase, N f = 3 in the CFL phase, and N f = 1 
in the spin-one color-superconducting phases. The constants are ( = (21 — 81n2)/54, 
a = (3 + 4 In 2)/27, and (3 = (6 - 4 In 2)/9. 
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Table 8: The Meissner mass for gluons, photons and from the mixed gluon-photon polar- 
ization tensor for various color-superconducting phases. The results are given in units of 
iV7/i 2 /(67r 2 ), where N f = 2 in the 2SC phase, N f = 3 in the CFL phase, and N f = 1 
in the spin-one color-superconducting phases. The constants are ( = (21 — 81n2)/54, 
a = (3 + 4 In 2)/27, and (3 = (6 - 4 In 2)/9. 
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Table 9: The diagonal elements of the electric and magnetic gluon-photon mass matrix 
and the (square of the) cosine of the rotation angles for Debye and Meissner masses. The 
results are given in units of NfH 2 /(6n 2 ), where Nf = 2 in the 2SC phase, Nf = 3 in the 
CFL phase, and Nf = 1 in the spin-one color-superconducting phases. The constants are 
C= (21-81n2)/54 and (3 = (6-41n2)/9. 
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one. This is in agreement with the general arguments presented in Sec. 12.4.21 and summarized in 
Table [T] The particular pattern of gluon masses reflects the residual SO(3) c+ j symmetry in the CSL 
phase, cf. Table [D the gluons corresponding to the three antisymmetric generators of [SU (3) c ] (which 
are simultaneously generators of £0(3)) assume a different mass than the ones corresponding to the 
symmetric generators. 

The final step is to diagonalize the mass matrices m 2 D ab , m 2 M ab for electric and magnetic gluons. Since 
only the eighth gluon mixes with the photon, this diagonalization is achieved by a simple orthogonal 
rotation in the 2x2 block corresponding to the indices a = 8, 9. The resulting diagonal (squared) 
Debye masses m 2 D a and (squared) Meissner masses m 2 M a , as well as the (square of the) cosine of the 
rotation angles 9d, 8m are shown in Table El In the case of an unbroken [C/(l)] symmetry, cf. Tabled 
the "rotated" (magnetic) photon is massless. The "rotated" gluon remains massive, but its degeneracy 
with the other massive gluons is lifted. 

The case of the polar phase is special. If there is only one quark flavor, or all quark flavors have the 
same electric charge, the results shown in Tables[7|-|n]hold. In this case, the rotated photon is massless, 
and there is indeed an unbroken [t/(l)] symmetry. However, in the case of two or more quark flavors 
with different electric charges, the results change [52] ■ Let us assume that the chemical potentials of all 
quark flavors are identical. Then the (squared) gluon masses are the same, but in the mixed masses, 
the factor q has to be replaced by J2f Qfi while in the photon masses, the factor q 2 is replaced by J2f Q 2 - 
In this case, it is not hard to realize that a diagonalization of the Meissner mass matrix does not lead 
to a massless rotated photon. There is therefore no symmetry. Consequently, spin-one color- 

superconducting quark matter exhibits a Meissner effect, while color superconductors with spin- zero 
Cooper pairs do not. Due to the smallness of the ratio 4>o/m g for spin-one color superconductors, these 
are very likely of type I, i.e., the magnetic field is completely expelled. This is in contrast to the standard 
model of a neutron star, where the core is assumed to be a type-II superconductor and thus threaded 
by magnetic flux tubes. It was recently argued in Ref. |162j that the short precession period of some 
pulsars contradicts this assumption and requires the core of the pulsar to be a type-I superconductor. 
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The question then is whether the core could be made of spin-one color-superconducting quark matter 



6 Conclusions and Outlook 

In this review, I have presented the current knowledge of the equilibrium properties of strongly in- 
teracting matter at large temperatures and/or densities. In particular, I have qualitatively discussed 
the phase diagram. I have presented calculations of thermodynamic properties of strongly interacting 
matter, both via lattice QCD, as well as within analytic approaches. Finally, I have given an overview 
of color superconductivity in weak coupling. 

Our knowledge of the QCD phase transition and the QGP at zero quark chemical potential has 
tremendously increased over the last few years. Lattice QCD calculations are well under control for the 
pure [SU (3) c ] gauge theory, and the quality of the data is such that an extrapolation to the continuum 
limit as well as to the thermodynamic limit has become possible. Lattice calculations with dynamical 
fermions are more challenging and, consequently, the data are not of the same quality as for the pure 
gauge theory. The main problem is that, with present methods of putting fermions on the lattice, 
the pion comes out too heavy. Since pions dominate the equation of state in the hadronic phase, 
calculations of the pressure below the chiral restoration temperature do not yet reflect the correct 
physics. The challenge for the future is to improve the methods such that the pion mass on the lattice 
is close to the value in nature. Besides a reliable computation of the equation of state, this will also 
allow to decide the question about the order of the QCD phase transition in nature. 

Another important development in lattice QCD is to extend the investigation of thermodynamic 
properties to nonzero quark chemical potentials. For many years, the fermion sign problem has impeded 
progress in this direction. Recent attempts, like multiparameter reweighting, Taylor expansion around 
\x = 0, or analytic continuation from imaginary values of /i, have made an attempt to work around this 
problem. Much work remains to be done to improve these methods in order to correctly determine the 
location of the critical point in the (T, fi) plane. This is of great phenomenological importance: in order 
to find a signal for the first-order phase transition to the QGP in nuclear collisions, one has to tune the 
bombarding energy such that one probes the region of the phase diagram, which is to the right of the 
critical point. Nuclear collisions at very high energies most likely probe the region to the left, i.e., the 
crossover region of the quark-hadron transition. By definition, there is no qualitative difference between 
hadronic and QGP phase in this region, and a clear signal for the QGP will be hard to identify. 

Although the question about the location of the critical point is important and can be investigated 
with the above mentioned lattice methods, these methods circumvent the fermion sign problem rather 
than solving it. Moreover, they are only applicable for quark chemical potentials fi from zero up to 
values of order T. Therefore, one will ultimately have to find a true solution which also works at small 
temperatures and large chemical potential, so that the color-superconducting quark matter phase can 
be explored. 

Analytic approaches to compute the equation of state of strongly interacting matter at high tem- 
perature have advanced rapidly in recent years. The equation of state is now known to all orders which 
are perturbatively computable. Work is in progress to determine the remaining nonperturbative con- 
tribution of order 0(g 6 ). Resummation techniques have been applied to compute the thermodynamic 
properties of strongly interacting matter. At large temperatures T ^> T c , they suggest that the QGP 
is a gas of weakly interacting quasiparticles. However, when approaching the critical temperature from 
above, the approaches based on resummation techniques fail to describe lattice QCD data. At the 
moment, one seems to be forced to either abandon field-theoretical rigor in favor of simple quasiparticle 
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models with sufficiently many fit parameters to reproduce the data, or turn to an alternative descrip- 
tion, such as the Polyakov loop model, which is physically less intuitive. It remains to be shown how 
this model is related to the quasiparticle picture at large temperatures. 

Color superconductivity is a rapidly evolving field. It is fairly likely that color-superconducting 
quark matter can be found in the core of compact stellar objects. It remains to explore how this 
phase influences the properties of the star. Much work has still to be done, for instance to compute 
the transport properties of color-superconducting matter and the phase diagram under the conditions 
of electric and color neutrality. Although NJL-type models may give a qualitative picture of possible 
scenarios, they are unreliable when one wants to draw quantitative conclusions. The task is to improve 
existing weak-coupling calculations or apply nonperturbative techniques to obtain further knowledge 
about this interesting, exotic phase of strongly interacting matter. 
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